零代码复现SCI!P trend+多模型控制混杂一站式搞定!

avatar
作者
猴君
阅读量:0

今天给大家带来的复现内容是CHARLS数据库的一篇横断面研究文章,文章统计方法部分包括了基线差异性分析、②线性趋势分析计算P trend、③多模型线性回归控制混杂因素

这些统计方法风暴统计统统一站式搞定,零代码操作,新手小白也可以轻松上手完成!


案例文献是一篇基于CHARLS公共数据库的一项横断面研究,旨在了解SUA和FPG与认知功能(3个层面)之间的关联。

442709df23c74b21d3c060147872fac1.pngSUA和FPG与认知功能的关联:一项横断面研究

1.1 基线差异性分析

首先,合并SUA与FPG分组【仅Low SUA、仅High FPG、无Low SUA和High FPG(Non)以及Low SUA和High FPG(Both)】,对人群的基本特征进行描述。

abcd7386e2c0d31b262b21eec501563f.png

1.2 多因素线性回归+趋势性P trend

table2展示的是自变量SUA及协变量与认知功能结局(3个层面:Global Cognition、Episodic Memory、Mental Status)的多因素回归结果。

table3展示的是自变量FPG及协变量与认知功能结局(3个层面:Global Cognition、Episodic Memory、Mental Status)的多因素回归结果。

这两个表格分别探讨了SUA与FPG对认知功能的单独影响。

2d5af72546a177b703424ddbf3d3078e.png

5f65d85435130439d3da0ce6036ff7d0.png

1.3 多模型线性回归

合并SUA与FPG分组【仅Low SUA、仅High FPG、无Low SUA和High FPG(Non)以及Low SUA和High FPG(Both)】,通过多模型策略逐步校正混杂因素,探讨两者对认知功能的联合影响。

82a0215a78ac27451613862cc305ac02.png

以上就是本篇CHARLS数据库文章的主要内容啦!下面我们将通过风暴统计平台进行结果复现!


风暴统计平台是由浙江中医药大学郑卫军教授基于R语言开发的统计分析平台,它具有以下优点:

  • 结果输出简洁,直接生成三线表,支持word版下载!

  • 平台是基于R语言进行开发,结果准确性有保障

  • 全部菜单式操作,统计小白也可以轻松上手!

a53ea99eaff25350b20cef44cb9dcd7d.png

2.1 复现数据介绍

我们复现用到的数据集是根据原文献变量及纳排标准整理的,样本量n=6875(原文献6509),样本量有所出入,但是变量和原文保持一致,包含有:

2f78b5545ed5885019e3a81f7d719cbe.png

便利起见,针对这篇文章的统计方法制作了专用链接,按顺序操作即可复现全部统计方法,不用在风暴统计平台换模块分析啦!

https://shiny.medsta.cn/xujm/

2.2 数据导入

进入网页后,可以选择导入自己的数据进行分析,也可以选择示例数据(也就是我们本篇文献的复现数据)进行测试。

f3c6c8d17d515d68675be0b3ad05cd8f.png

2.3 基线差异性分析

来到"分组差异性分析"板块,分别选入分组变、正态定量变量、偏态定量变量、分类变量。

72967a2f31f4010cb07f94f525eae0bb.png

平台就会直接给出基线差异性结果啦!平台还给出了差异性分析结果的统计量,信息更全面。由于原文中,二分类变量仅展示阳性人群的分布比例,因此表格长度会更简短些。

6b93cfc8ec6ce296964feb5608b3c220.png

2d583314fe5566f8afe96ae0b685925d.png

2.3 线性趋势分析计算P trend

趋势性分析是针对等级数据开展的,连续性变量需要提前转换为等级变量


选入我们想要开展分析的连续型变量,平台给出了多种分组方式,这里我们选择"四分位数",接着点击"开始分组"

e5c761b5ba7713f5345a0e480bee51a1.png

右侧预览界面就会新增两列数据!

以“quantile”为后缀的变量名,代表四分位数转换后的等级数据,将原始的连续型变量赋值为1,2,3,4。

以"M_value"为后缀的变量名,则是用每组的中位数重新赋值。这个变量就是计算P for trend的核心变量了!

0d1268be562bd721e5085f0a31c353f6.png

首先,选入定量结局变量(3个层面的认知功能,以global cognition为例),以"M_value"为后缀的组中值分类变量。

最最重要的是,取消勾选"分类变量"选项!使变量以定量形式进入回归分析!

338f860e6cd5f8757c8e0f89d704c04a.png

这时的P值就是我们的P trend值了!

81f7d6848a889add671c1653cb014901.png

如需校正协变量,可以点击"增加回归模型",依次选入协变量即可。

94f6aac3fe0190ddf0be0a48a5335611.png

2.4 多因素线性回归

接着来到"线性回归"模块进行logistic回归分析!

依次选入结局变量(以global cognition为例),根据文献选入协变量,点击多因素回归查看多因素结果。

87c7a8ae9bba9b42c7fe1d7fbcadb35c.png

平台依旧是直接给出三线表结果,原文献中省略了reference行以及P值,SE值,t值等,平台全部给出,大家可以根据自身需要进行取舍。

a8fe64b342d21a5e99b875df2dae1fbe.png

2.5 多模型线性回归

接着来到多模型线性回归策略模块,首先选入我们的结局变量,以及我们SUA和FPG联合变量,接着点击"增加回归模型",选择"自定义",按照原文献勾选我们的协变量。

32c6eedeedcc11cd3e3c8177d529e47e.png

接着右侧结果框就会出现我们的多个模型的回归结果啦!下方脚注还会标明我们校正了哪些协变量!

d1385f5182f46963577b82446f4277a2.png

与原文献对比,可以看出整体格式形式比较相似,还免去了绘制三线表填写结果的繁琐!

aa2ab85c6467cbd64ccd0cdeaeb8e1ca.png

以上就是我们本次零代码复现SCI文章的全部内容啦!有P trend 、多模型控制混杂等常见统计分析方法!欢迎试用下方链接,进行复现与测试哦!

https://shiny.medsta.cn/xujm/

如果您在使用过程中出现报错,可以参考下方推文,排查一下问题所在哦!

答疑 | 风暴统计更换服务器后,为什么总出现同样的报错?


详情请点击下方:

临床、护理研究1对1全程统计指导,支持发表SCI论文或顺利毕业

    广告一刻

    为您即时展示最新活动产品广告消息,让您随时掌握产品活动新动态!