(event = 1,指的是确实事件;0是死亡,1是缺失
1.用的是ovarian自带数据 ,首先画出生存曲线survfit(),但是曲线上不显示确实值的符号+号,怎么在函数中设置合适的缺失值???
plot()设置mark.time = TRUE
2接下来对rx进行log-rank检验
p取的是0.05,所以拒绝原假设,两组实验没有显著性区别,
若是有显著性,可以通过中位生存期,或是RR值比较两种情况对生存状况的影响;
3.log-rank只是对单因素进行了检验,建立回归模型,考虑多因素的情况
(中途换了一组数据--恶性肿瘤的预后)
从结果中可以看出,X--(4,5,6)有显著性,RR值即为exp(coef),ceof即为β i值
4,library("MASS"),fit2<-stepAIC(fit1,direction = "both")
可以看出,模型中只有2和3时是最好的
plot(survfit(fit2))