最近我用R的survIDINRI包做模型的NRI和IDI计算,但是因为变量是factor,于是进行了哑变量处理,但是处理后运行命令报错,代码和报错内容如下,求助大神帮忙解释是哪里出了问题。
代码
# 对variable2进行哑变量处理
MM<-model_matrix(biomarkers,deadfollowtime~variable2-1)
new<-cbind(std[,1:3],MM)
std<-as.matrix(subset(data2,select = c(“time”,”event”,”variable1”)))
n<-nrow(data2)
covs1<-as.matrix(new[,c(-4,-5,-7)])
covs2<-as.matrix(std[,c(-1,-2)])
t0<-max(data2[,3])
x<-IDI.INF(new[,1:2],covs2,covs1,t0,npert = 300)
IDI.INF.OUT(x)
报错
> x<-IDI.INF(new[,1:2],covs2,covs1,t0,npert = 300)
Error in unoecdf(cc, pdiff[case], Wi[case] * PTB.Vi[case]) :
外接函数调用时不能有NA/NaN/Inf(arg3)