Kaplan-Meier法又称乘积极限法(Product-Limit method),前面我推送的《SPSS案例实践:Kaplan-Meier生存分析》给大家介绍了具体用法,今天图文主要介绍它的生命表中累积生存率(生存率)的计算过程。
上表是某个KM生存分析所得各时点生存率表格,第一组生存时间分布中有18例个案,第3列是时点(单位:天),第5列是累积生存率,具体是如何计算的呢?
第1个时点(第7天)
总18例,该时点死亡1例,此时的生存概率=17/18=0.9444,因为是第一个终点时间,此时累积生存率就是第1年的生存概率0.9444。
第2个时点(第8天)
上一时点尚有17例随访,其中1例在本期死亡,剩余16例,本期生存概率=16/17=0.9412,由于这个时点生存实际上是在上一个时点基础或条件下产生的,因此累积生存率是基于上一个生存率来计算的,因此随访到第8天的累积生存率=0.9444*(16/17)=0.8889。
在没有遇到失访个案下,以此类推进行计算。
第13个时点(第80天)
前一个即第12时间点累积生存率0.333,随访到第13时点即80天时出现一例死亡,剩余5个个案,80天时点生存概率=5/6=0.8333,截止到80天的累积生存率=0.8333*(0.333)=0.2777。
第14个时点(第83天)
第83天失访1例,从当前剩余随访对象总数中减去1例,即还剩余4例。因为不知该失访对象最终结局,所以无法计算本期本时点的生存率,确切的累积生存等仍然为上一期的0.2777,下一个失效事件发生后再行计算。
第15个时点(第84天)
上一点即83天失访1例后剩余随访对象4例,其中有1例在第84天死亡,84天时点生存概率=3/4=0.7500,随访截止到84天的累积生存率=0.2777*(3/4)=0.2082。
第16个时点(第90天)
第90天死亡1例,时点生存概率=2/3=0.6667,累积生存率=0.2082*(2/3)=0.1388。
第18时点(第186天)
随访到第186天时,最后1例死亡,当期时点生存概率0/1=0,累积生存率=0.069*0=0。本组所有个案1例失访,17例死亡。随访结束。
======
请大家根据上表,以及我给出的计算示范,自己再试算一下,加深对Kaplan-Meier生命表累积生存率的计算过程的理解。
本文完
图/文=数据小兵
大家好,我是数据小兵,一名统计学/数据分析知识博主。长期从事统计学方法应用、统计答疑咨询、统计软件教学培训与数据分析工作。这是我的博客网站,欢迎大家长期关注,感谢支持。
以下是我在网易云课堂发布的视频课程:
《SPSS统计分析视频教程》
《量表问卷数据分析》
《AMOS结构方程视频教程》
《JASP统计分析视频教程》
《正交与响应面试验数据分析》
《Gpower+PASS样本量计算》
《Stata与Meta分析入门》
《SPSS+Excel论文图表》
视频课程二维码地址:
如果没有时间看视频教程,就是想快速解决个别小问题小疑惑,欢迎直接到我的“知乎咨询”小额付费答疑解决。“知乎咨询”答疑网址如下:
https://www.zhihu.com/consult/people/723499335798489088
“知乎咨询”答疑手机二维码地址如下: