问题如题,报错的为”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) ;