公众号
“医学论文与统计分析”公众号
扫码关注公众号

统计咨询
“公共数据库与孟德尔随机化”公众号
扫码关注公众号

意见反馈
邮箱:17357190071@163.com
微信:aq566665

生存数据如何做中介分析?中国学者用NHANES数据库发文一区top

Administrator
发布于 2025-07-31 / 95 阅读
0
0

甘油三酯-葡萄糖指数(TGI)作为一个可靠且易于获取的指标,可有效反映人体胰岛素抵抗程度。而幽门螺杆菌作为一种常见的致病菌,可通过多种途径影响人体健康。

有研究表明,幽门螺杆菌感染(幽门螺杆菌免疫球蛋白G(IgG)抗体阳性)可能会加剧人体胰岛素抵抗,从而促进心血管疾病(CVD)的发生和发展。

2024年12月18日,中国学者用NHANES数据库,在期刊《Cardiovascular Diabetology》(医学top一区,IF=8.5)发表题为:“Relationships among Helicobacter pylori seropositivity, the triglyceride-glucose index, and cardiovascular disease: a cohort study using the NHANES database”的研究论文,旨在探讨了幽门螺杆菌IgG抗体阳性、TGI和CVD风险之间的关联,以及TGI在这一关联中的中介作用。

研究结果表明,TGI在幽门螺杆菌增加心血管疾病(CVD)风险和死亡风险的过程中,均发挥了显著的中介作用。同时,研究还发现在一般人群中,TGI≥第75个百分位数且感染幽门螺杆菌者患CVD风险最高。而在CVD患者中,该特征的人群死亡风险最高。

在公共数据库与孟德尔随机化公众号回复 原文 即可获得文献PDF等资料。想用NHANES发文,看看这个可一键提取和分析数据的NHANES  Online平台!如感兴趣请联系郑老师团队,微信号:aq566665

研究团队基于美国国家健康与营养调查(NHANES)III(1988-1994)和NHANES数据库1999-2000年的数据,经过纳排,最终纳入了9,399名符合条件的年龄≥18岁的参与者,平均年龄为45岁,50.5%为女性。

  • 其中,4,488名(47.75%)为幽门螺杆菌IgG抗体阳性,3,934名(41.86%)被诊断为CVD。

图1 研究人群筛选流程图

研究团队根据TGI值(≥第75个百分位数为高TGI水平、<第75个百分位数为低TGI水平)以及是否感染幽门螺杆菌(幽门螺杆菌血清状态[阳性或阴性])对参与者进行分组。

主要暴露变量是TGI和幽门螺杆菌IgG抗体阳性。

主要结局指标为CVD风险及其全因死亡率(ACM)。

CVD包括冠心病(CHD)、心肌梗死、中风、充血性心力衰竭(CHF)、心绞痛和高血压。

图2 研究流程

主要研究结果

√TGI、幽门螺杆菌血清学状态与CVD死亡率之间的关联

研究团队首先进行了加权单变量和多变量logistic回归,探究TGI水平、幽门螺杆菌IgG抗体状态及其组合与CVD之间的关联。

在调整了混杂因素后,结果显示,在一般人群中,高TGI水平以及高TGI水平且感染幽门螺杆菌者患CVD风险显著升高。

表1 在一般人群中,TGI和幽门螺杆菌与CVD风险的关联

模型 1:未调整;

模型 2:根据年龄、性别和种族进行调整;

模型 3:根据年龄、性别、BMI、种族、教育程度、吸烟状况、HDL、T2DM进行调整

在多变量logistic回归的基础上,团队进一步探讨TGI、幽门螺杆菌IgG抗体状态与CVD之间潜在的非线性关系。RCS结果表明,在一般人群中,仅TGI和CVD风险之间存在显著的线性关系。

图3 TGI和幽门螺杆菌与CVD风险之间的RCS曲线

同时,研究团队通过生存分析来评估CVD患者中TGI水平、幽门螺杆菌IgG抗体阳性与CVD患者死亡风险之间的关系。

  • KM曲线结果显示,在CVD患者中,高TGI水平且感染幽门螺杆菌者死亡风险最高。Cox比例风险模型结果与其一致,进一步证明了该结果。

图4 不同TGI水平和幽门螺杆菌IgG抗体状态下,CVD患者中ACM的Kaplan-Meier曲线

表2 在CVD患者中,TGI以及幽门螺杆菌与ACM的关联

在此基础上,研究团队再次绘制了RCS曲线,探讨CVD患者TGI水平、幽门螺杆菌IgG抗体阳性与死亡风险之间可能存在的非线性关系。

结果显示,在CVD患者中,也是仅TGI与死亡风险之间存在显著的线性关系。

图5 TGI水平和幽门螺杆菌与全因死亡风险之间的RCS曲线

√亚组分析

为了进一步验证研究结果,团队还根据糖尿病前期、糖尿病状态、代谢综合征和肥胖对CVD患者进行了分层,并开展了亚组分析。

分析结果与上述一致,表明高TGI水平且感染幽门螺杆菌的患者面临最高的全因死亡风险

√中介效应分析

最后,中介分析的结果揭示了,TGI在幽门螺杆菌增加CVD风险和死亡风险的过程中,均发挥了显著的中介作用。

这意味着幽门螺杆菌可能通过影响胰岛素抵抗等机制,间接增加患CVD风险和CVD患者的死亡风险。

图6 中介分析结果

综上所述,幽门螺杆菌可通过TGI水平间接提高心血管疾病风险和死亡风险。此外,将幽门螺杆菌与TGI水平相结合,可以提高对个体患心血管疾病风险以及心血管疾病患者的死亡风险的预测能力。

统计知识点汇总

1. 数据来源

美国国家健康与营养调查(NHANES)III和NHANES数据库:

NHANES数据库是一项旨在评估美国成人和儿童的健康和营养状况的研究计划。该数据库涵盖了多种健康指标的数据收集和分析。其中,NHANES III,即第三次国家健康和营养调查,特指在1988年至1994年期间进行的这一轮调查。

2.研究设计思路:

P(Population)参与者:从NHANES III(1988-1994)和NHANES数据库1999-2000年的数据中,纳入的9,399名接受幽门螺杆菌血清学检测并提供空腹血液样本的年龄≥18岁的参与者;

E(exposure)暴露因素:甘油三酯-葡萄糖指数(TGI)和感染幽门螺杆菌

O(outcome)结局:心血管疾病(CVD)风险及其全因死亡风险(ACM);

S(Study design)研究类型:横断面和队列研究相结合。

3. 统计方法

  • 加权单变量和多变量Logistic回归

  • 变量Logistic回归:指在模型中只包含一个自变量的情况。这种类型的模型用来评估单一因子对结果概率的影响大小和方向。

  • 变量Logistic回归:包含两个或两个以上的自变量。这种模型允许研究者同时考察多个因素对结果变量的影响。

在本文中,研究团队先分别对TGI水平、幽门螺杆菌IgG抗体状态进行加权单变量logistic回归分析,以初步了解它们各自与CVD的关系。

加权处理是为了考虑样本的代表性,使得结果更能反映总体情况。

在此基础上,团队进一步考虑了其他可能影响CVD的协变量(如年龄、性别、种族、BMI、糖尿病状态等),并进行了加权多变量logistic回归分析。这样做可以更加准确地评估TGI和幽门螺杆菌IgG抗体状态在多个因素共同作用下的影响。

单变量与多变量Logistic回归的关系

变量和多变量Logistic回归之间的关系是包含与被包含的关系。单因素Logistic回归是多变量Logistic回归的一个特例,后者是更为一般的形式。

在实际应用中,研究者通常首先进行单变量Logistic回归分析,以识别可能的重要预测因素。然后,在多变量Logistic回归中包含这些因素,以控制混杂变量,准确估计每个因素的效应。

  • 生存分析

生存分析是一种用于分析事件发生时间的统计方法,其关注的重点是从观测开始到目标事件发生的时间,旨在了解和预测这些事件发生的时间及其影响因素。

常见的生存分析方法包括:Kaplan-Meier曲线、Log-rank检验和Cox比例风险模型。

本文主要用了这个三个方法,评估CVD患者中TGI水平、幽门螺杆菌IgG抗体阳性与CVD全因死亡率之间的关系。

限制性立方样条(RCS)曲线

RCS曲线是一种平滑的曲线拟合方法,可以用于评估自变量与因变量之间的非线性关系。

在本文的研究中,研究团队绘制了两次RCS曲线。

  • 一般人群中,暴露与CVD风险之间的RCS曲线;

  • CVD患者中,暴露与全因死亡风险之间的RCS曲线。

  • 亚组分析

亚组分析是将研究人群按照某种特征划分为不同的亚组,然后分别对每个亚组进行分析,以探究该特征是否会对主要关联结果产生影响。

在本文中,研究团队在每个亚组内,运用加权多变量Cox比例风险模型分析,以评估TGI水平、幽门螺杆菌IgG抗体状态在该亚组中与CVD及全因死亡风险的关系。

  • 中介效应分析

介效应分析旨在探究自变量X如何通过中介变量M,对因变量Y产生影响的作用机制。通过该方法,我们可以更加明确X对Y的影响是否以及是如何通过中介变量M间接发生,而不仅仅是直接的因果关系。

在本文的研究中,研究团队构建了一个中介效应模型,探究了TGI是否是感染幽门螺杆菌增加CVD风险及全因死亡风险的中介变量及其中介效应大小,该研究步骤通过'mediation' 包实现。

√生存数据的中介分析

本文中用到的就是中介分析中最常用的R包mediation。该R包使用简单,只需分别运行中介模型和结局模型,然后使用mediate函数获得中介分析结果,可分析多种类型的中介和结局,但是只能实现单中介模型。

在公共数据库与孟德尔随机化公众号后台回复“中介分析”即可获得该资源。

研究团队采用了特定的公式来计算中介效应的大小,考虑了直接效应和间接效应(即通过TGI的中介作用)对结果变量的影响。

中介分析的过程包括以下四个步骤:

(1)验证X与Y之间的相关性(模型 Y = βTot X),其中βTot表示总效应;

(2)建立了X与M之间的关联(模型 M = β1 X),其中β1表示间接效应;

(3)确定了经X调整后的Y部分(模型 Y = β2 m + βDir X),其中β2代表间接效应,βDir代表直接效应;

(4)计算中介效应的百分比:中介效应(%)=(β1 × β2 / βTot)× 100%。

这篇文章的分析思路全面但清晰易懂,对于想用NAHNES数据库发文的同学而言,是一篇值得参考的好文章。想了解更多有关NHANES数据库发文思路的同学,郑老师的NHANES一对一统计服务课程值得了解一下!现在报名还会送一年的NHANES Online平台使用权,像本研究中涉及的NHANES指标TGI,平台可零代码一键提取和分析数据。


评论

收银台

订单信息

应付金额 积分

模块介绍

请加我们助教二维码或加入Zstats使用群
Zstats交流群

Zstats交流群

助教二维码

联系助教

选择支付方式

请输入助教告诉您的积分券

如果不填写积分券,将直接使用当前余额支付

正在创建订单...

请稍候,正在为您生成支付订单

支付二维码

请使用扫描二维码完成支付

等待支付中...

二维码获取失败

支付二维码获取失败,请点击重新获取

正在处理余额支付...

请稍候,正在为您完成支付

正在处理充值并支付...

正在使用积分券兑换,然后完成支付 正在使用当前余额完成支付

支付成功!

您的订单已支付完成,页面将在 秒后自动关闭

支付失败

支付过程中出现错误,请重新选择支付方式

平台说明