【零代码复现SCI】多模型生存分析+累积发生曲线！

今天给大家带来的复现内容是CHARLS数据库的一篇横断面研究文章，文章统计方法部分包括了①基线差异性分析、②多模型生存分析+累计发生曲线、③非趋势性分析。

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

今日文章分为两部分

1. 复现文章介绍

2. .风暴统计平台快速复现

1.复现文章介绍

---

案例文献是一篇基于CHARLS公共数据库的一项横断面研究，旨在了解中国中老年人群内脏脂肪指数与中风发病风险之间的关联。

中国中老年人群内脏脂肪指数与中风发病风险的关系：一项大型全国队列研究的证据

1.1 基线差异性分析

首先，将研究的焦点暴露CVAI进行四分位数分组，并作为分组变量，对人群的基线特征进行简单描述。

1.2 绘制累积发生曲线

以CVAI四分位数作为分组变量，绘制KM曲线，并提供了log-rank检验P值。

1.3 多模型策略+趋势性分析计算P trend

采用多模型策略，逐步矫正混杂因素，研究CVAI四分位数变量与结局之间的关系。

1.4 非趋势性分析，绘制RCS曲线

将CVAI作为连续变量，绘制RCS曲线，探讨与结局间的关系，得出P for non-linear trend=0.319，两者之间不存在非线性关系。

1.4 亚组分析（多模型策略）

在不同血压分组的人群中，采用多模型策略，分别构建多个模型，研究CVAI四分位数变量与结局的关系。

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

2.风暴统计平台快速复现

---

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

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

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

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

2.1 复现数据介绍

我们复现用到的数据集是根据原文献变量及纳排标准整理的，样本量有所出入，但是变量和原文保持一致，包含有：

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

https://shiny.medsta.cn/limy/

2.2 数据导入

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

2.3 基线差异性分析

首先，在"数据整理"模块将CVAI转换为四分位数变量。

操作完成后，数据集就会新增一列四分位数变量。

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

根据原文献，将连续变量全部选入"正态变量"，分类变量选入"分类变量"板块。

选择完毕后，平台就会直接给出基线差异性结果啦！和原文献对比发现，提供的信息很全面且一致，并且额外提供了统计量！

2.3 绘制累积发生曲线

根据界面提示，分别选入结局变量、时间变量、分组变量。这里结局变量，建议以0和1进行赋值，0代表删失事件、1代表死亡事件。

接着右侧结果界面就会出现我们的累积发生曲线啦！

2.4 多模型策略+趋势性分析计算P trend

首先选入我们的结局变量、生存时间与焦点暴露因素

接着选入我们的model2、model3需要调整的协变量。

然后我们的多模型回归结果就出现啦！

2.5 RCS曲线（非线性趋势分析）

接着来到“非线性趋势RCS图”模块，因为本篇文章是生存分析，因此选择"COX回归"。

分别选入结局变量、生存时间变量、连续型自变量以及我们的混杂协变量。

接着右侧就会出现调整协变量后的RCS曲线啦！

以上就是我们本次零代码复现SCI文章的全部内容啦！有累积发生曲线、生存分析多模型控制混杂 、RCS曲线（校正混杂因素版）等常见统计分析方法！欢迎试用下方链接，进行复现与测试哦！

https://shiny.medsta.cn/limy/

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

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