寿命表法计算中位生存期
寿命表法计算中位生存期
生存分析中，一般以中位生存期的显著差异作为一个影响因素重要性的判断标准。
一般拿到生存数据之后，可以基于生存状态，和生存时间进行KM分析，获取中位生存时间。
比如下面的生存分析表（部分截图）：
可以采用如下代码进行生存分析：
# 估计KM生存曲线 survial <- Surv(clin$time,clin$status==1) kmfit <- survfit(survial~1) kml <- summary(kmfit) # 查看中位生存期 kmfit #Call: survfit(formula = survial ~ 1) # # n events median 0.95LCL 0.95UCL ## 307 60 2888 2052 NA
从结果可知，中位生存时间为2888天。生存曲线如下：
# 绘制生存曲线 plot(kmfit,conf.int = F,mark.time = F) #绘制一条中位生存期对应红线 abline(h=0.5,lty=2,col="red")
如果您对生存分析感兴趣，可以学习我的TCGA系列课程：
- 发表于 2018-11-16 15:08
- 阅读 ( 9615 )
- 分类：GWAS
