R语言cox实现

最新推荐文章于 2024-03-11 12:30:27 发布

shlay

最新推荐文章于 2024-03-11 12:30:27 发布

阅读量8.3k

点赞数 5

CC 4.0 BY-SA版权

分类专栏： R语言文章标签： R语言生存分析 COX实现

本文链接：https://blog.youkuaiyun.com/s1164548515/article/details/89574460

R语言专栏收录该内容

8 篇文章

订阅专栏

R语言cox实现

加载包


library("survival")

没有的先安装

install.packages('survival')

示例

testdata <- data.frame(time=c(3,5,6,7,9),status=c(1,0,1,0,1),z=c(0,1,1,0,1))
resCox<-coxph(Surv(time,status)~z,data=testdata)
summary(resCox)

结果

Call:
coxph(formula = Surv(time, status) ~ z, data = testdata)

  n= 5, number of events= 3 

     coef exp(coef) se(coef)      z Pr(>|z|)
z -0.5493    0.5774   1.4179 -0.387    0.698

  exp(coef) exp(-coef) lower .95 upper .95
z    0.5774      1.732   0.03585     9.297

Concordance= 0.667  (se = 0.222 )
Likelihood ratio test= 0.15  on 1 df,   p=0.7
Wald test            = 0.15  on 1 df,   p=0.7
Score (logrank) test = 0.15  on 1 df,   p=0.7