Как установить данные для анализа выживаемости R - PullRequest
1 голос
/ 27 сентября 2019

Извините, если вопрос выглядит глупо, но у меня есть небольшой набор данных, которым я бы хотел манипулировать с помощью функции "Survfit" пакета R "Survival", и, ну, я не знаю, чтобы установитьправильный фрейм данных, который можно использовать с помощью «Surfit»;данные выглядят следующим образом:

      time number_at_risk number_death number_censored
    1   25             10            0               2
    2   28              8            1               0
    3   33              7            1               0
    4   37              6            0               1
    5   41              5            1               0
    6   43              4            0               1
    7   48              3            0               3

И теперь, если я запускаю обычный синтаксис survfit(Surv(time, number_censored) ~ 1, data = data), он выдает предупреждение In Surv(time, number_censored) : Invalid status value, converted to NA.

Очевидно, что данные не организованы должным образом,Итак, как мне установить мой фрейм данных?Спасибо.

1 Ответ

1 голос
/ 27 сентября 2019

time должен быть вектором со временем, когда событие произошло, и status показателем, если это событие является цензурой или смертью (0/1).

В вашем примере данные должны выглядетькак это:

times = c(1,1,2,3,4,5,6,7,7,7)
status = c(0,0,1,1,0,1,0,0,0,0)

survfit(Surv(times,status)~1)
...