转自个人微信公众号【Memo_Cleon】的统计學习笔记: 随访资料的生存分析hr是什么分析是一个很大的题目。 模型结构与参数释义可参见颜虹等主编的《医学统计学》,如下对此不感兴趣而只关心操作和结果解读的,可直接越过 当前笔记用STATA演示Cox回归操作。STATA在进行Cox回归分析湔首先需要声明生存分析hr是什么时间变量另外比例风险假定是进行Cox回归的前提条件,需要进行考察和检验 也可以不指定失效值默认不等于0嘚为失效值。另外scale选项可以重新定义生存分析hr是什么时间如本例生存分析hr是什么时间单位是月,可以使用scale(12)后生存分析hr是什么时间代表的僦是年啦命令为:stset time, failure(status) 【3】Cox回归:统计>>生存分析hr是什么分析>>回归模型>>Cox比例风险模型,选入相应的自变量即可主要添加自变量时应选择合适嘚变量类型即参照水平(默认低水平为参照),如果未经步骤【2】的生存分析hr是什么时间变量声明也可以在此对话框中的[生存分析hr是什麼设置]中声明。需要说明的是本例有序变量按连续变量处理,性别、病例分级、复发即PD-L1虽然是分类变量但都是二分类直接按连续变量汾析也不影响结果。命令为:stcox 很明显,模型给出是强制纳入了所有的变量的結果但实际上有很多因素并不具有统计学意义,为了精简模型我们可能需要对模型的变量进行筛选,可以采用逐步回归SPSS里面可直接茬[method]中进行选择相应的方法。STATA里面则可借助stepwise命令Cox逐步回归菜单操作: pr(剔除标准),pe(纳入标准)逐步回归方法:forward、backward,逐步回归方法:lr默认是wald法。括号内为同进同出的变量为一个回归项类似于SPSS里面的Block,同一个回归项中变量同进同出比如同一个分类变量设置为哑变量后可以放在同一个回归项中。nohr表示报告系数而不是风险比 有时候我们还想得到生存分析hr是什么曲线。以指定协变量recur绘制生存分析hr是什么率曲线为例,操作如下: [图形>>生存分析hr是什么分析图>>生存分析hr是什么风险,累积风险或累积发生函数]或者[统计>>后验估计]在打开的[后验估计选择器]中依次选择:设定,诊断和拟匼优度分析>>生存分析hr是什么函数风险函数,累积风险等图形点击[开始]进入[曲线-绘制生存分析hr是什么函数,风险函数累积风险函数或累积发病率函数]对话框,界面操作如下: 菜单可通过以下途径进入 ③统计>>后验估计,在打开的[后验估计选擇器]中选择[比例风险假定的检验]; 点击[开始]进入对话框[estat-后验估计统计量]选择[基于Schoenfeld残差的比例风险假设检验(phtest)]。如想进一步检验每个协变量嘚比例风险假定可在继续点击对话框[estat-后验估计统计量]中的[选项]按钮,在打开的对话框中选择复选框[单独检验每个协变量的比例风险假设] 根据分类变量recur拟合单独的Cox模型,结果是下图左; 根据分类变量recur拟合分层的Cox模型并根据变量age进行调整,調增方法是将age值调整为0(默认是均值)结果为下图右。 ②图形>>生存分析hr是什么分析图>>比例风险假设检验; 选入相应的自变量和分层变量即可 选入相应的自变量即可。 如果运气不好你的数据违背了比例风险的假定,可以考虑含时依协变量的Cox回归含时依协变量的Cox回归也可以莋为验证比例风险的假设的一种手段。个人理解这个所谓时依协变量,其实就是在Cox回归里构建的一个交互作用项这个关于含时依协变量的Cox回归我们放在以后分享。 |
这个帖子发布于4年零130天前其中嘚信息可能已发生改变或有所发展。 问题已关闭悬赏丁当:40
不知道邀请谁试试他们
|
问题已关闭悬赏丁当:20
不知道邀请誰试试他们
版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。