Средняя выживаемость

Я пытаюсь использовать пакет выживания в R, чтобы запустить тест логарифмического ранга для средних показателей выживаемости - возможно ли это? Я использую код:

survdiff(surv(Age, Survival) ~ Treatment, data = GBM)

где Возраст = возраст на момент постановки диагноза; Выживание = сколько лет они прожили; Лечение = получали ли они лечение (кодируется как 1 или 2)

Если это правильно, почему я получаю сообщение об ошибке:

Error in survdiff.fit(y, groups, strata.keep, rho) : 
  There is only 1 group
In addition: Warning message:
In Surv(GBM$Age, GBM$Survival) : Invalid status value, converted to NA

Нужно ли кодировать лечение как 0 или 1 вместо 1 или 2?


r survival-analysis
person user2726449    schedule 23.11.2013    source источник
comment
ваш индикатор выживания/события неверен. иди посмотри на ?surv   -  person rawr    schedule 24.11.2013
comment
Вы правы... Я пропустил событие, и оно, похоже, сработало... не уверен, правильно ли оно вычисляется   -  person user2726449    schedule 24.11.2013

Ответы (1)


arrow_upward
0
arrow_downward

В pkg:survival нет функции surv. Это Surv.

?Surv

Не имеет значения, кодируется ли ковариата 0/1 или 1/2. Оценка будет описывать отношение риска, зависящее от разницы в одну единицу. Если вы хотите настроить вызов Surv для обработки возраста и продолжительности, тогда время события должно быть «Возраст + выживание». И, конечно же, вам понадобится индикатор событий.

person IRTFM    schedule 24.11.2013