R语言做nomo图报错 Error in lims[[i]] : subscript out of bounds

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，我自己才接触R语言不久，现在在尝试做nomo图，但是一直有报错，代码如下：
dd<-datadist(train.plus.rad_score)
options(datadist="dd")
nom<-nomogram(crmodel, fun=list(function(x) surv(24, x)), lp=F, funlabel=c("6 month survival"), maxscale=10, fun.at=c(0.95, 0.85,0.75, 0.65, 0.5))

之后就有报错
Error in lims[[i]] : subscript out of bounds
In addition: Warning message:
In structure(limits, class = "data.frame", row.names = c("Low:effect",  :
  Calling 'structure(NULL, *)' is deprecated, as NULL cannot have attributes.
  Consider 'structure(list(), *)' instead.

请问各位有遇到这个问题吗？非常感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：script Bounds scrip Error Bound

这个报错的意思是你lims调用的地方超出了定义列表的长度。比如你一个list有5个元素，你调用了第20个。

owenqi 发表于 2019-11-27 10:47
这个报错的意思是你lims调用的地方超出了定义列表的长度。比如你一个list有5个元素，你调用了第20个。

非常感谢，解决啦！

tomato7 发表于 2019-11-28 10:13
非常感谢，解决啦！

您好，我也遇到了相同问题，能向您请教一下怎么解决的么？

tomato7 发表于 2019-11-28 10:13
非常感谢，解决啦！

您好，请问这种情况怎么解决呢？

kaimuo 发表于 2020-3-17 11:04
您好，请问这种情况怎么解决呢？

需要重新用其他函数重建一下，才行

tomato7 发表于 2020-3-23 09:41
需要重新用其他函数重建一下，才行

谢谢您的回复，那请问是建模的函数需要重建么？还是怎么弄呢？十分感谢呀！方便留一个邮件或微信的联系方式么~冒昧了~

tomato7 发表于 2020-3-23 09:41
需要重新用其他函数重建一下，才行

谢谢您的回复，那请问是建模的函数需要重建么？还是怎么弄呢？十分感谢呀！方便留一个邮件或微信的联系方式么~冒昧了~

kaimuo 发表于 2020-3-23 16:46
谢谢您的回复，那请问是建模的函数需要重建么？还是怎么弄呢？十分感谢呀！方便留一个邮件或微信的联系方 ...

很可能打包的数据框中各别变量（列）中有缺失值“NA”，
检查一下。

绘制列线图的打包数据框中各变量（列）中不能有缺失值，
即便不用于拟合模型的变量（列）也不要有缺失值。

kaimuo 发表于 2020-3-23 16:46
谢谢您的回复，那请问是建模的函数需要重建么？还是怎么弄呢？十分感谢呀！方便留一个邮件或微信的联系方 ...

用coxph()函数会出现这个错误，
但是用cph()就不会报错。

具体例子的代码见下：

1. install.packages("rms")
   
2. library(rms)
   
3. 
   
4. lung1 = datadist(lung)
   
5. options(datadist = "lung1")
   
6. 
   
7. fitlung = cph(Surv(time, status) ~ age + sex + ph.karno,
   
8.                 data = lung, x = TRUE, y = TRUE, surv = TRUE)
   
9. #cph()中要加 surv=TRUE
   
10. surv = Survival(fitlung)
    
11. surv1 = function(x)surv(365, x)
    
12. surv2 = function(x)surv(730, x)
    
13. 
    
14. nomolung = nomogram(fitlung,
    
15.                     fun = list(
    
16.                                surv1, surv2), lp=F,
    
17.                     fun.at = c(0.05, seq(0.1, 0.9, by = 0.05), 0.95),
    
18.                     funlabel = c('1 year survial', '2-year survival'))
    
19. plot(nomolung)

复制代码