Cox风险回归分析

(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))




你可能感兴趣的:(生物统计)