原代碼如下：
tempfit <- cph(fit$formula,data=https://bbs.csdn.net/topics/mydata,x=T,y=T,surv = T)
surv<-Survival(tempfit)
surv1<-function(x)surv(365,x)
surv2<-function(x)surv(365*2,x)
surv3<-function(x)surv(365*3,x)
surv4<-function(x)surv(365*4,x)
surv5<-function(x)surv(365*5,x)
surv6<-function(x)surv(365*6,x)
surv7<-function(x)surv(365*7,x)
surv8<-function(x)surv(365*8,x)
surv9<-function(x)surv(365*9,x)
surv10<-function(x)surv(365*10,x)
nom <- nomogram(tempfit, fun=flist[1:yr],
fun.at=c(seq(.1,.9, by=.1)),
funlabel=fl)
產生了10年的surv，根據輸入的yr在列線圖中產生yr個生存概率線，但如果超過10年呢？如果換成月難道要一個個敲出來嗎，不可能的
于是我使用如下的for回圈想要生成每個surv[i]
su <- list()
for(i in 1:10){
su[i]<-function(x)surv(365*i,x)
}
flist <- list(su1,su2,su3,su4,su5,su6,su7,su8,su9,su10)
但是出錯如下
Error in su[i] <- function(x) surv(365 * i, x) :
cannot coerce type 'closure' to vector of type 'list'
一開始搜索以為surv沒被賦值，被系統當做函式的surv了，我就改成了su，但還是錯誤，我覺得是su[i]<-function(x)surv(365*i,x)這句的右邊surv(365*i,x）這個surv出錯，但是又不知道怎么改，球球老哥們幫幫忙
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/280393.html
標籤：其他開發語言
上一篇：ABAP常見陳述句解釋（）