医学统计学生存分析.ppt

(3) 死于与研究疾病无关的原因：终止随访时间为死亡时间。 (4) 研究终止。研究终止时观察对象仍然存活。终止随访时间为研究终止时间。 3．影响生存的有关因素：如患者年龄、病情、病程、术前健康等情况，以便分析这些因素对生存率的影响。 第三十页，共五十七页，2022年，8月28日 二、随访方式 1．全部观察对象同时接受处理措施，观察到最后一例出现结果或事先规定的随访截止时间。 图3 随访资料常见形式示意图 第三十一页，共五十七页，2022年，8月28日 2．观察对象在不同时间接受处理措施，完成 一定数量随访病例或按事先规定的时间停止随访。 图4 随访资料常见形式示意图 第三十二页，共五十七页，2022年，8月28日 生存率的估计与生存曲线： 一、乘积极限法 乘积极限法(product-limit estimate)又称Kaplan- Meier法，适用于未分组生存资料的分析。 例16-1 14例膀胱肿瘤3.0cm患者和16例膀胱肿瘤≥3.0患者的生存时间（月）如下，试估计两组各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时间，并绘制生存曲线。 第三十三页，共五十七页，2022年，8月28日 肿瘤3.0cm 14 19 26 28 29 32 36 40 42 44+ 45 53 + 54 59 + 肿瘤≥3.0cm 6 7 9 10 11 12 13 20 23 25 27 30 34 37 43 50 表3 肿瘤3.0cm组生存率及标准误的计算 生存时间 t 死亡数 dt 期初病例数 nt 截尾数 ct 死亡概率 qt 生存概率 pt 生存率 S(t) 生存率标准误 SE[S(t)] 14 1 14 0 1/14=0.0714 0.9286 0.9268 0.0688 19 1 13 0 1/13=0.0769 0.9231 0.8572 0.0935 26 1 12 0 1/12=0.0833 0.9167 0.7858 0.1097 28 1 11 0 1/11=0.0909 0.9091 0.7144 0.1207 29 1 10 0 1/10=0.1000 0.9000 0.6429 0.1281 32 1 9 0 1/9=0.1111 0.8889 0.5715 0.1323 36 1 8 0 1/8=0.1250 0.8750 0.5001 0.1336 40 1 7 0 1/7=0.1429 0.8571 0.4286 0.1323 42 1 6 0 1/6=0.1667 0.8333 0.3571 0.1281 44 0 5 1 0/5=0.0000 1.0000 0.3571 0.1281 45 1 4 0 1/4=0.2500 0.7500 0.2678 0.1233 53 0 3 1 0/3=0.0000 1.0000 0.2678 0.1233 54 1 2 0 1/2=0.5000 0.5000 0.1339 0.1130 59 0 1 1 0/1=0.0000 1.0000 0.1339 0.1130 第三十四页，共五十七页，2022年，8月28日 生存率的计算 1．生存时间t：由小到大排列，遇非截尾和截尾值相同，截尾值排后。 2．死亡数dt：与生存时间t对应。注意：截尾值对应的个体未发生“死亡”事件，故死亡数为0。 3．期初病例数nt，表示恰好在该时刻以前的病例数。如n29为10，表示恰好在29月时点前有10人存活。 4．死亡概率qt，表示t月前的观察对象恰好在t月时点死亡的概率。 第三十五页，共五十七页，2022年，8月28日 5．生存概率pt，表示t月前的观察对象恰好在t月时点存活的概率。 6．生存率S(t)。表示该人群恰好活过t 时刻的概率。它为小于和等于 t 时刻的各时点生存概率的乘积。 7．生存率的标准误SES(t) 。 第三十六页，共五十七页，2022年，8月28日 假定生存率近似服从正态分布，某时点总体生存率的（1－a）％可信区间，公式为： 总体生存率的可信区间 本例28月总体生存率的95％可信区间： 即膀胱肿瘤3.0cm患者28月生存率的95%可信区间为47.78％~95.10%。 生存曲线尾部的生存率不适合于用该法计算总体生存率的可信区间。 第三十七页，共五十七页，2022年，8月28日 由表3可见，中位生存时间估计