Я пытаюсь вычислить C-статистику модели пропорциональных рисков Кокса, которая обычно идет довольно гладко с помощью dynpred :: cindex. Однако теперь я сталкиваюсь с ошибкой, с которой я никогда не сталкивался. Ниже приведен воспроизводимый пример.
Воспроизводимый пример
# 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-статистику с предоставленным набором данных, поскольку они могут быть реалистичнымиданные.
Спасибо, Стив