我正在尝试计算考克斯比例风险模型的C统计量，通常使用dynpred：：cindex进行得非常顺利。然而，现在我遇到了一个我以前从未遇到过的bug。请在下面找到一个可复制的示例。

可复制示例

# Produce a random dataset, columns: linear predictor (-1,+1), survival (1-60), vital status (0/1)

TestData <- setnames(as.data.frame(
    cbind(runif(n=20, min=-1, max=1),
    runif(n=20, min=1, max=60),
    round(runif(n=20, min=0, max=1)))),
    c("LP","Overall_survival","Vital_status"))

cindex(Surv(Overall_survival, Vital_status)~LP, data=TestData)

在这个可重复的代码中，我随机遇到（大约50%的时间我运行这个代码）以下错误。

遇到dynpred：：cindex错误

Error in if (time[j] > time[i]) { : missing value where TRUE/FALSE needed

虽然我不明白为什么dynpred：：cindex不应该总是能够用提供的数据集计算C统计量，因为这些数据可能是真实的数据。
谢谢，史蒂夫