全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4234 4
2017-08-24
  问题如题,报错的为”Error in if (pest$IDI > 0) { : missing value where TRUE/FALSE needed“,具体代码如下,不知道到底是哪里的问题,希望有做过的大神帮忙看一下。谢谢~~~~

dat<-read.csv("Data2.csv",head=T,sep=",")
DA=subset(dat, select=c("time","status","age","sex"))
DA=DA[!is.na(apply(DA,1,mean)),] ; dim(DA)
t0=5*365
mydata=DA
indata1=mydata;
indata0=mydata[,-4]; n=nrow(DA);
covs1<-as.matrix(indata1[,c(-1,-2)])
covs0<-as.matrix(indata0[,c(-1,-2)])
indata3<-as.matrix(mydata[,1:2])
#--- inference ---
x<-IDI.INF(indata3, covs0, covs1, t0, npert=200) ;
#--- results ---
IDI.INF.OUT(x) ;

二维码

扫码加我 拉你入群

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

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

全部回复
2017-8-25 09:43:15
自己顶一下,求助求助
二维码

扫码加我 拉你入群

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

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

2018-12-7 17:36:29
楼主问题解决了么,是不是你的数据里有缺失值?我想咨询一下你带入模型的数据是factor还是numeric?
二维码

扫码加我 拉你入群

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

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

2019-5-7 19:07:26
我也遇到同样的问题,说明书里面案例可以正常运行,但是自己的数据试了半天都没有结果(包括R升级、数据类型统一为num),感觉是不是数据的问题,不知道最后楼主怎么解决的?
思考了8个小时,终于知道哪里出问题了。很关键的一个地方是t0的设置。
如果大家的生存数据是按月来随访的,t0应该和你的最长随访时间差不多,应该是30、40、50之类的。
如果生存数据是按年来随访,t0应该为1、2、3、4、5之类的。
如果是按天来随访,那么t0为90、365、365*n等。
希望对大家有用!

二维码

扫码加我 拉你入群

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

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

2023-6-11 19:33:34
同样在用 survIDINRI包求NRI和IDI,当模型中包含分类变量时,说明书推荐使用model.matrix修改哑变量,我也删除了额外产生的intercept变量,可修改之后还是不停报错,各位有什么办法么?
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群