, надеюсь, кто-нибудь может мне помочь.
Я создал зависимые от времени кривые ROC для конкретного маркера, используя пакет timeROC. Я хочу найти оптимальное ограничение для этого маркера в определенное время, предпочтительно используя «индекс Юдена». К сожалению, я не могу найти какую-либо функцию в R, которая может помочь мне определить эту точку вырезания для кривых ROC, зависящих от времени.
Пакеты 'cutpointr' и 'оптимальной точки резки' работают только для обычных ROC и не зависят от времениROC (если я прав).
в основном мои данные выглядят так:
- маркер <- c (1.774, 1.731, 1.667, 1.751, 1.551, 1.693, 1.765,1.717, 1.769, 1.696, 1.680, 1.607, 1.575, 1.851, 1.665, 1.710, 1.790, 1.503, 1.571, 1.812) </li>
- ID <- c (1:20) </li>
- time_change <-с (3,03, 3,0, 3,04, 2,0, 10,3, 0,98, 3,09, 1,5,6,96, 6,96, 7,15, 3,17, 4,1, 9,05, 1,98, 11,1, 3,1, 2,99, 4,98, 6,5) </li>
prog_var <- c (0,0,0,1,0,0,0,1,1,1,1,0,0,0,0,0,0,0,1,0) </p>
df <- data.frame (ID, prog_var, marker, time_change) </p>
и моя формула для создания зависимых от времени ROC:
- ROCMarker <- timeROC (T = df $ time_change, дельта = df $ prog_var, маркер = df $ маркер, причина = 1, времена = c (1,2,3)) </li>
ROCMarker содержит: TP, FP, AUC, время, совокупная заболеваемость, SurvProb, N и некоторую другую информацию, которая мне действительно не нужна.
Может кто-нибудь дать мне какой-нибудь совет? Спасибо!