Ожидается тип функции R predict.cox

Я использую модель Кокса для прогнозирования.

fit <- coxph(Surv(time, status) ~ age + ph.ecog + strata(inst), lung)
predict(fit,type="expected")

В документе predict.coxph говорится, что ожидаемый тип — это «ожидаемое количество событий с учетом ковариат и времени последующего наблюдения». Что означает время последующего наблюдения здесь? Как я могу это изменить?


r survival-analysis cox-regression
person Yoki    schedule 30.06.2015    source источник
comment
В случае, представленном без аргумента newdata, последующее время представляет собой просто переменную time.   -  person IRTFM    schedule 30.06.2015
comment
Спасибо, если я хочу предсказать новые данные, как я могу изменить время?   -  person Yoki    schedule 30.06.2015

Ответы (1)


arrow_upward
3
arrow_downward

Это пример «ожидаемого» для одного конкретного набора характеристик пациента с выбранным единственным временем, превышающим среднее время выживания для этой когорты, но с использованием всех уникальных кодов учреждения для набора данных «легкие».

 (pr.fit.new <- predict(fit, newdata=data.frame(age=50, time=365, status=0, 
                                ph.ecog=2,inst=unique(lung$inst[!is.na(lung$inst)]) ),
                        type="expected"))
 [1] 1.1886094 1.7679429 1.6651195 1.3117278 0.7743925 0.8501810
 [7] 1.9686954 1.2560981 2.5845069 0.6513932 1.1241870 1.1515346
[13] 1.1616482 1.9071766 4.2468797 0.4921904 0.5875339 0.5647238

Слои имеют заметно различающиеся «базовые опасности», что приводит к существенным различиям в «ожидаемых» для гипотетических лиц с одинаковыми ковариантами и периодом наблюдения.

person IRTFM    schedule 30.06.2015