Функция IDI.INF для сравнения двух моделей регрессии Кокса, ошибка?

Question

Я хочу использовать функцию IDI.INF из пакета aliveIDINRI для сравнения производительности двух моделей регрессии Кокса с использованием базы данных из 15000 участников. Я использовал этот скрипт:

D=subset(Dataset, select=c("time","status","Var1","Var2","Var3", "Var4", "Var5", "Var6", "Var7", "Var8", "Var9","Var10","Var11","Var12","Var13"))
D$status=as.numeric(D$status==2)
D=D[!is.na(apply(D,1,mean)),] ; dim(D)
Dataset=D[1:15930,]

outcome=D[,c(1,2)]
covs1<-as.matrix(D[,c(-1,-2)])
covs0<-as.matrix(D[,c(-1,-2,-13)])

t0=80
x <- IDI.INF(outcome, covs0, covs1, t0, npert=200)

В конце я получаю следующее сообщение об ошибке:

Error in unoecdf(cc, pdiff[case], Wi[case] * PTB.Vi[case]) : 
  NA/NaN/Inf in foreign function call (arg 5)

Кто-нибудь знает, что влечет за собой это сообщение об ошибке? Я попытался использовать разные переменные и сократить количество переменных, которые я использую, а также изменить время t (0) среди прочего без успеха.

Answer 1

Если ваша переменная status является дихотомической (то есть 0, 1), то измените D$status=as.numeric(D$status==2) на D$status=as.numeric(D$status==1).