Stata基础:生存分析

学术   2024-12-01 11:29   中国  

01

生存分析的基本介绍

生存分析 (survival analysis) 又被称为久期分析 (duration analysis) 、转换分析 (transition analysis)、风险分析 (hazrd analysis)、信度分析(reliability analysis)、失效时间分析(failure-time analysis)、历史事件分析 (event history analysis) 等,主要研究 关注事件 发生所需要的时间,具体可归纳为三类 (Paul D. Allision,2018):
·定性变量状态变化所花费的时间, 比如离婚、晋升、死亡
·定量变量出现急剧变化所经历的时间, 比如人口总和出生率锐减
·定量变量超过某个阈值所需要的时间, 比如将体重减到健康标准
生存分析的基本步骤:
·第一步 设定生存分析数据
·第二步 画生存函数、累积风险函数和风险函数
·第三步 进行参数回归及非参数回归
·第四步 比例风险假定的检验

02

数据来源

使用的数据是来自 Stata-press 官网的 hip.dta 数据(右击可以另存到本地)。该数据集主要调查一种新的可穿戴充气防护装置是否可以保护老年人,即降低老年人发生跌倒骨折的风险。邀请了 48 位年龄为 60 岁以上的女性展开一项试验,任意挑选 28 名女性发放并让她们穿戴防护装置,作为实验组, 剩下的 20 名女性则不穿防护装置作为控制组。
研究中的关键变量定义如下:
·关注事件 为骨折 (fracture)
·持续时间 为开始被调查到发生骨折(fracture = 1) 的时间, 同时也包括发生过骨折以后恢复了再次发生骨折的时间 (被试会被调查几个阶段,不同阶段之间会有时间间隔)    
·解释变量 包括:年龄 (age)、血钙(calcium) 和 穿戴防护装置 (protect)。

03

stata代码

3.1 将数据设定为生存分析格式
3.2 描述性分析 (非参数分析)
3.3 模型估计    
3.4 模型的检验
比例风险模型和 Cox PH 模型的重要假设是, 如果这个假设不成立则上述模型不成立。因此,必须要对比例风险模型设定进行检验。下面介绍三种检验方式。

资源获取方式


1、点亮文章底部“在看”图标


2、关注“社科学术汇”公众号


点击公众号右上角三个点-设为星标



3、点击“发消息” 在后台发送关键词

生存分析

按照提示获取资源哦


扫码添加客服微信|领取课程大额优惠券

   

社科学术汇
传播前沿知识信息,助力社科实证研究
 最新文章