生存分析及cox比例风险模型.ppt

生存分析 生存分析的主要内容 生存过程的描述 生存时间的分布特点 估计生存率 估计平均生存时间 绘制生存率曲线 生存过程的统计推断 生存时间的影响因素研究 生存分析需具备的资料 生存时间或称时间变量 观察结果或称截尾变量 暴露变量及混杂因素 肾移植例 生存资料的定义(1) 将数据定义为生存资料格式 stset 时间变量 [，failure(截尾变量)] failure (截尾变量) stata以所有不等于0的非缺失值为预期的出现结果； 若failure选择项缺省，则stata默认所有纪录均出现结果 生存资料的定义(2) 四个变量---定义数据库后，stata系统自动产生的 _st /* 数据库中该条记录是否被定义为生存资料 _d /* 该条记录是否出现预期结果 _t /* 观察对象被随访的时间 _t0 /* 观察对象第一次被观察到的时间(开始过程的时间为0) stset time,failure(outcome) failure event: outcome != 0 outcome . obs. time interval: (0, time] exit on or before: failure ------------------------------------------------------------------ 33 total obs. 0 exclusions ------------------------------------------------------------------ 33 obs. remaining, representing 18 failures in single record/single failure data 11370 total analysis time at risk, at risk from t = 0 earliest observed entry t = 0 last observed exit t = 1045 生存资料的描述 stci ,[by(分组变量)选择项] p(#) /* 生存时间的百分位数 level(#) /* 可信区间的可信度 media /* 计算中位生存时间，为默认选项 rmean /* 计算平均生存时间 生存率的估计 输出生存率、生存概率的标准误等统计量 sts list [,by(分组变量) 选择项] failure /* 输出死亡函数 (1-S(t)) level /* 规定所输出可信区间的可信度 绘制生存（死亡）曲线(1) 绘制Kaplan-Meier生存(死亡)曲线 sts graph [,by(分组变量) separate 选择项] failure /* 指定绘制“死亡”曲线，缺省时绘制生存曲线 Gwood /* 绘制生存(或死亡)曲线及其Greenwood可信区间 lost /* 在曲线上标出该时间点截尾值例数 绘制生存（死亡）曲线(2) 绘制Kaplan-Meier生存曲线与Cox预测曲线 stcoxkm [,by(分组变量) separate 绘图命令选择项] 生存率的比较 sts test 分组变量 [,选择项] logrank /* 进行log-rank检验，为缺省值 例8.8 某临床试验比较A，B两治疗方案对某病的治疗效果，A组(group=0)12人，B组(group=1)13人。病人分组后检验其肾功能(kidney)，功能正常者记为0，不正常者记为1；治疗后生存时间为stime(天) ；问不同治疗方案及肾功能对病人的生存时间是否有影响? Cox比例风险模型 (1) 估计Cox比例风险模型的命令格式为 stcox [协变量] [,nohr 选择项] nohr /* 指定输出回归系数