看到丫姐画完这篇1区28分JCO里所有的图这篇推文后，对里面出现的nomogram列线图以及用于评价模型的图十分眼馋，于是决定试着画画看
1. 使用rms包的cph()建立cox模型
library(survival)
library(rms)
data(package="survival")
dd<-datadist(lung)
options(datadist="dd")
f <- cph(Surv(time,status)~ age+sex+ph.karno, data = lung, x=T, y = T, surv = T)
survival <- Survival(f)
survival1 <- function(x)survival(730,x)
nom <- nomogram(f, fun = survival1, fun.at = c(0.1,seq(0.1,0.9,by = 0.2), 0.9), funlabel = "2 years survival")
plot(nom)
结果
2. 使用survival包coxph()建立模型
library(regplot)
library(survival)
data(pbc)
pbc$catbili <- cut(pbc$bili,breaks=c(-Inf, 2, 4, Inf),
labels=c("low","medium","high"))
pbc$died <- pbc$status==2
pbccox <- coxph(formula = Surv(time,died) ~ age + catbili + sex +
copper +stage + trt,data=pbc)
regplot(pbccox,observation=pbc[1,], failtime = c(365,1095,1825), prfail = TRUE )
结果
大家看看医咖会的这篇吧手把手教你用R画列线图(Nomogram)及解读结果
好像不是很难吼，大家画起来
