生存分析

走出图形操作的圈子,在脑袋里思考成型

我要画个热图,查个5年生存率

1.洞悉你的数据,精雕细琢

2.选个R包,写个函数,查看结果

做生存分析,用survival包,surfit函数

rt=read.table("riskTrain.txt",header=T,sep="\t") 
fit <- survfit(Surv(futime, fustat) ~ risk, data = rt) 
summary(fit) #查看总体生存率
或res.sum <- surv_summary(fit) # 查看每个时间节点的生存率
view(res.sum) # 查看每个时间节点的生存率,比如我要看3、5年生存率,那么我就需要这么做

surfit函数的本质是K-M分析,想想一下之前做的SPSS操作,时间、状态、因子,其实是一个道理。


K-M.png

如果想做整体的、不分因子的生存分析,需要做的就是把因子risk,换成1

fit <- survfit(Surv(time, status)~1, data=rt)

更多的生存分析知识,明天再跟这个帖子学习https://www.jianshu.com/p/4ad9ba730719

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容