|
> library(foreign)
> lc<-read.spss('E:/skip.sav',use.value.labels=T,to.data.frame=T)
> library(survival)
> coxm<-coxph(Surv(DFS,DFSstatus)~Gender+Age+Excision+Tstatus+Pathology+ECE+N2stations+Skip+PRtotal+PRN2,data=lc)
> scoxm<-step(coxm)
> library(rms)
> dd<-datadist(lc);options(datadist="dd")
> surv<-survfit(scoxm)
> surv<-function(x)surv(5*12,lp=x)
> nomogram(scoxm,fun=list(surv),lp=F,funlabel=c('5-year survival'),maxscale=10,fun.at=c(0.95,0.9,0.85,0.8,0.75,0.7,0.6,0.5))
Error in lims[[i]] : 下标出界
求助各位大侠。这段语言鼓捣了一个月了。。。。。。。 |
|