Cox回归分析及其SPSS操作方法概述
我们先回顾一下生存分析的KM法和寿命表法(回复数字26和27可以查看KM法的详细内容),其共同点是只能分析一种因素与生存率的关系,Log-Rank法也是比较一个因素两种水平间的生存差别,如果生存率的影响因素有很多,我们怎么避免其它混杂因素的影响呢?我们可以使用回归分析方法,但如果使用logistic回归,也是只能观察影响因素与结局的关联,没有考虑结局发生的时间因素。Cox回归可以解决这个问题。Cox回归一般模型假设为
其中h(t,X)是在时刻t的风险函数又可称瞬时死亡率,h(0,t)是基线风险率,其它与logistic回归模型相同。βj大于0则xj越大,病人死亡风险越大,βj小于0则xj越大,病人死亡风险越小,βj等于0则xj越与死亡率没有影响。Exp(β)为危险比(HR)或相对危险度(RR)。
下面以一个例子说明在SPSS中作Cox回归如何操作。
我们想观察乳腺癌的生存率及其影响因素,收集了1207例病例并进行了随访。观察的因素包括年龄(age)、病理肿瘤大小(pathsize)、腋窝淋巴结个数(lnpos)、组织学分级(histgrad)、雌激素状态(er)、孕激素状态(pr)和淋巴结转移(ln_yesno)等。time为随访时间,status为生存状态。
在SPSS菜单里点击“分析”-“生存函数”-“Cox回归”,在弹出的对话框里,将”time”和” status”分别选入时间和状态对话框,点击“定义事件”,填写“1”,将不同的影响因素选入协变量框中,方法可以选“向后:LR”(各种方法差别不大,可以自由选择)。
如果有多分类变量需要设置哑变量,可以点击右上角“分类”,将要设置哑变量的变量选入右边框中。并可以选择以第一个或者最后一个作为参照。
在右上角点击“选项”,可以选择“CI用于exp(B)”,用于计算HR的95%置信区间。
最后点击确定可看到Cox回归分析结果。
结果中第一个表给出病例纳入情况,如下图,数据共1207个病例,但最后一共纳入590例,其中40例出现事件(即死亡),另外617例因为有缺失值被排除。从中可以看出,数据质量不太好,有缺失值的病例占一半以上且有观察终点的病例只有40例。
下面的表中是哑变量编码情况,histgrad中“1”被编码为“0”“0”,即histgrad中“2”“3”均以“1”为参照。
下面再看主要的结果,即“方程中的变量”表。本表列出了多个步骤,在步骤1中,全部我们纳入的变量都进入分析,从前往后分别是模型系数(B)、系数标准误(SE)、Wald检验值,自由度(df)、p值,HR值(Exp(B))及其置信区间。接下来看步骤2,其相对于步骤1少了一个变量er。即步骤2中删除了步骤1中的P值最大的变量。同理依次删除p值最最大的变量。
下面我们看最后一步,即步骤5.经过筛选,只剩下三个变量,即认为这三个变量对生存率的影响,其中病理肿瘤大小对应的HR为1.566,大于1,即认为病理肿瘤越大,生存时间越短;同理腋窝淋巴结个数越多,生存时间越短;孕激素状态对应的HR为0.511,小于1,即有孕激素时生存时间越长。
需要说明的是Cox回归分析是比例风险模型,即模型假设在任一时间点两组的危险比是相同的。如下图所示:
而下图所示则不符合比例风险模型,不能作简单Cox回归。如果想作回归分析,可以咨询相关统计专家或查看专业书籍。
数据分析咨询请扫描二维码
若不方便扫码,搜微信号:CDAshujufenxi
在当今数字化时代,数据分析师的重要性与日俱增。但许多人在踏上这条职业道路时,往往充满疑惑: 如何成为一名数据分析师?成为 ...
2025-04-02最近我发现一个绝招,用DeepSeek AI处理Excel数据简直太爽了!处理速度嘎嘎快! 平常一整天的表格处理工作,现在只要三步就能搞 ...
2025-04-01你是否被统计学复杂的理论和晦涩的公式劝退过?别担心,“山有木兮:统计学极简入门(Python)” 将为你一一化解这些难题。课程 ...
2025-03-31在电商、零售、甚至内容付费业务中,你真的了解你的客户吗? 有些客户下了一两次单就消失了,有些人每个月都回购,有些人曾经是 ...
2025-03-31在数字化浪潮中,数据驱动决策已成为企业发展的核心竞争力,数据分析人才的需求持续飙升。世界经济论坛发布的《未来就业报告》, ...
2025-03-28你有没有遇到过这样的情况?流量进来了,转化率却不高,辛辛苦苦拉来的用户,最后大部分都悄无声息地离开了,这时候漏斗分析就非 ...
2025-03-27TensorFlow Datasets(TFDS)是一个用于下载、管理和预处理机器学习数据集的库。它提供了易于使用的API,允许用户从现有集合中 ...
2025-03-26"不谋全局者,不足谋一域。"在数据驱动的商业时代,战略级数据分析能力已成为职场核心竞争力。《CDA二级教材:商业策略数据分析 ...
2025-03-26当你在某宝刷到【猜你喜欢】时,当抖音精准推来你的梦中情猫时,当美团外卖弹窗刚好是你想吃的火锅店…… 恭喜你,你正在被用户 ...
2025-03-26当面试官问起随机森林时,他到底在考察什么? ""请解释随机森林的原理""——这是数据分析岗位面试中的经典问题。但你可能不知道 ...
2025-03-25在数字化浪潮席卷的当下,数据俨然成为企业的命脉,贯穿于业务运作的各个环节。从线上到线下,从平台的交易数据,到门店的运营 ...
2025-03-25在互联网和移动应用领域,DAU(日活跃用户数)是一个耳熟能详的指标。无论是产品经理、运营,还是数据分析师,DAU都是衡量产品 ...
2025-03-24ABtest做的好,产品优化效果差不了!可见ABtest在评估优化策略的效果方面地位还是很高的,那么如何在业务中应用ABtest? 结合企业 ...
2025-03-21在企业数据分析中,指标体系是至关重要的工具。不仅帮助企业统一数据标准、提升数据质量,还能为业务决策提供有力支持。本文将围 ...
2025-03-20解锁数据分析师高薪密码,CDA 脱产就业班助你逆袭! 在数字化浪潮中,数据驱动决策已成为企业发展的核心竞争力,数据分析人才的 ...
2025-03-19在 MySQL 数据库中,查询一张表但是不包含某个字段可以通过以下两种方法实现:使用 SELECT 子句以明确指定想要的字段,或者使 ...
2025-03-17在当今数字化时代,数据成为企业发展的关键驱动力,而用户画像作为数据分析的重要成果,改变了企业理解用户、开展业务的方式。无 ...
2025-03-172025年是智能体(AI Agent)的元年,大模型和智能体的发展比较迅猛。感觉年初的deepseek刚火没多久,这几天Manus又成为媒体头条 ...
2025-03-14以下的文章内容来源于柯家媛老师的专栏,如果您想阅读专栏《小白必备的数据思维课》,点击下方链接 https://edu.cda.cn/goods/sh ...
2025-03-13以下的文章内容来源于刘静老师的专栏,如果您想阅读专栏《10大业务分析模型突破业务瓶颈》,点击下方链接 https://edu.cda.cn/go ...
2025-03-12