IDI.INF 函数比较两个 Cox 回归模型，错误？答案

【问题标题】：IDI.INF function to compare two Cox regression models, error?IDI.INF 函数比较两个 Cox 回归模型，错误？
【发布时间】：2019-01-22 20:32:38
【问题描述】：

我想使用 survIDINRI 包中的 IDI.INF 函数来比较两个 Cox 回归模型（使用 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) 等但没有成功。

【问题讨论】：

标签： r cox-regression

【解决方案1】：

如果您的 status 变量是二分变量（即 0、1），则将 D$status=as.numeric(D$status==2) 更改为 D$status=as.numeric(D$status==1)。

【讨论】：

【解决方案2】：

就我而言，它已将数据设置为小标题：

# include in IDI.INF help example code to invoke "bug"

library(tibble)  
mydata=as_tibble(mydata)

【讨论】：

【解决方案3】：

在我的例子中，问题是通过从数据中省略 NA 来解决的。

【讨论】：