走出图形操作的圈子,在脑袋里思考成型
我要画个热图,查个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