Визуализация результатов модели Вейбулла

Анализ выживаемости в R

Heidi Seibold

Statistician at LMU Munich

Визуализация модели Вейбулла

Инструменты визуализации часто работают со ступенчатыми функциями.

Поэтому следующий код НЕ работает:

wb <- survreg(Surv(time, cens) ~ 1)
ggsurvplot(wb)

Визуализация модели Вейбулла

wb <- survreg(Surv(time, cens) ~ 1)

Кривая выживаемости:

surv <- seq(.99, .01, by = -.01)
t <- predict(wb, type = "quantile", p = 1 - surv, newdata = data.frame(1))

surv_wb <- data.frame(time = t, surv = surv, 
                      upper = NA, lower = NA, std.err = NA)

График:

ggsurvplot_df(fit = surv_wb, surv.geom = geom_line)

Давайте потренируемся!

Анализ выживаемости в R