生存分析在SPSS上的实践

文摘   2024-10-04 15:58   安徽  


PSYCH统计实验室


前言

对于R使用不是很熟练的同学,SPSS上提供了完整的生存分析操作程序,可以辅助我们更快得到想要的答案。


一、进行K-M检验


找到K-M检验键入,选择变量进入分析框:

其中时间为样本的生存时间,状态为样本在结束随访时的状态,并对其进行定义:

由于死亡案例标注为1,而我们在分析过程中重点关注因假设而死亡的病例,所以单值设置为1。

这里的因子就是在整个随访过程中的分组方式,可以检验任意的假设分组。

组间比较的方式可以将三种都选中,其中Log-Rank与布雷斯洛的区别在上期推送中有相关说明。

在保存中,可以选中以上内容,来作为后续结果的参考。

在选项中,我们可以选择以上需要展示的内容。

第一部分的结果:

首先是对sex进行分组时组间是否具有统计学差异(在自己的实验中,也可能为干预方式分组等等)

可以很明确地看到,这三种统计学方式的检验结果都不具有显著差异。

此图为生存分析函数形式,其看法与R中做出的结果看法相同。因为两组线条有相交,所以两组之间不具有显著差异。

二、 因素分析


分析其他影响生存时间分布的因素,也同样按照此前的单因素与多因素的方式分别进行。

首先,找到COX回归的按钮。

打开后输入时间,状态,并定义事件,再输入协变量(单因素cox回归可以先输入一个,再将单因素具有统计学意义的筛选入多因素cox回归中。)

图中需要输出生存分析图。

保存生存函数,Xbeta比例系数(>1,危险因素;<1,不是危险因素)

选项中重点选择左边模型统计中的95%置信区间以及估算值相关性。

结果解读:

需要报告-2对数似然。

显著性方面<0.05,代表SO2越低越容易死亡。 其HR值为2.372。单因素cox分析得出此变量可以进入多因素cox回归。

三、 多因素cox回归:cox比例风险模型


执行步骤与之前单因素cox回归相同,但在变量选择中,输入多个协变量。


对分类变量进行分类协变量标记

并在图中选择分类变量进行单独线条绘制。

结果:

最终得到cox比例风险模型中的变量情况。(此时sex由于没有统计学意义已经被模型迭代删除。)

此时的生存模型如上图所示。

但由于分类协变量无意义,所以未得出根据分类变量进行分类时的双分布图。

当分类有统计学意义时,可以得到如下图像:


PSYCH统计实验室

通知公告

网络分析课程目前开放视频课啦

单次课200元/讲(学生),250元/讲(非学生)

共有四讲内容:

①横断面网络分析简介与基础

②网络分析与因子分析

③交叉滞后网络分析

④时间序列网络分析

购买后开放视频权限14天,可多次申请。

并赠送所有课程相关资料(无PPT)

如果想申请购买,请联系M18812507626


更多资讯

关注我们

文稿:莲花清瘟

排版:Little Star

责编:Wink
审核:摘星

本文由“Psych统计自习室”课题组原创,欢迎转发至朋友圈。如需转载请联系后台,征得作者同意后方可转载。

Psych统计自习室
大家好,我们是由来自北京师范大学,西南大学,天津医科大学等高校在读硕士、博士研究生组成的一个科研团队——Psych统计自习室。Psych统计自习室旨在关注心理学、精神病学领域的最前沿的系列研究,并做前沿统计知识的分享。
 最新文章