> library(foreign)
> es<-read.csv(file.choose(),header=T)
> library(survival)
> library(rms)
> coxm<-coxph(Surv(time,status)~.,data=es)
Warning message:
In fitter(X, Y, strats, offset, init, control, weights = weights, :
Loglik converged before variable 8 ; beta may be infinite.
> scoxm<-step(coxm)
.......
> dd<-datadist(es)
> options(datadist("dd"))
> surv<-function(x)surv(1*6,lp=x)
> nom<-nomogram(scoxm,fun=list(surv),funlabel=c("6-month survival"),maxscale=10,fun.at=c(0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1),lp=0)
Error in lims[] : 下标出界
求高手指导,出错的地方红笔指出,总是提示“下标出界”?不知道问题出在哪里?