求：用R软件做列线图的脚本

doctorZhang

用SPSS做的COX分析，想做个评分系统，看文献里用R软件做列线图可以实现，但确实没学过，照猫画虎的写程序，怎么也搞不定；有没有哪位大侠给写个现成的模板？我下载的R是Win7，3.2.2版本的。
下面是我在R软件里写的几条
> 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")
> sur<-survival(scoxm)
错误: 没有"survival"这个函数。

出现上面这个错误，就进行不下去了。

哪位大侠帮忙给看看。多谢。

englianhu

1. > ls('package:survival')
   
2. > ls.str('package:survival')

复制代码

查看程序包内的函数与数据。

englianhu

1. > ls('package:survival')
   
2. > ls.str('package:survival')

复制代码

查看程序包内的函数与数据。

doctorZhang

englianhu 发表于 2015-12-31 14:43
查看程序包内的函数与数据。

> 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(3*12,lp=x)
>nomogram(scoxm,fun=list(surv),lp=F,funlabel=c('3-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[] : 下标出界***************
请问这个个“下标出界”怎么处理？谢谢。