本人是医学生,最近在做nomogram.
但是code总是报错,不知道哪位大神能帮我指点一下?
下面是我的代码:
n<-data.frame(y.surv.va,rad,clinical)
dd<-datadist(n)
options(datadist="dd")
y<-Surv(n$time,n$status == 1)
coxm <-cph(y~rad+N,x=T,y=T,data=n,surv=T)
dd<-datadist(n)
options(datadist=”dd”)
scoxm<-step(coxm)
surv<-Survival(scoxm)
nom <- nomogram(scoxm, fun=list(function(x) surv(36, x), function(x) surv(60, x), function(x) surv(120, x)),
lp=F, funlabel=c("3-year survival", "5-year survival", "10-year survival"), maxscale=10,
fun.at=c(0.95, 0.9, 0.85, 0.8, 0.75, 0.7, 0.6, 0.5))
运行到最后nomo图的绘制的时候总是报错。
而且在step(coxm)这一步也会报错。
相同的程序我用R自带的lung数据集就能够成功画出nomogram.
rad和N这两个参数我分别做过Cox单因素回归,都是是有统计学意义的,p值小于0.05
希望大家可以一起讨论研究一下nomogram的代码。
谢谢啦~