Calcular un modelo Weibull
Para un modelo Weibull con covariables, podemos calcular la curva de supervivencia igual que para el modelo Weibull sin covariables. Lo único que tenemos que hacer es especificar los valores de las covariables de la curva de supervivencia en la función predict(). Esto se puede hacer con el argumento newdata.
Este ejercicio forma parte del curso
Análisis de supervivencia en R
Instrucciones del ejercicio
- Calcula un modelo Weibull para el conjunto de datos
GBSG2con la covariablehorThpara analizar el efecto de la terapia hormonal en la supervivencia de las pacientes. - Calcula la curva de supervivencia para las pacientes que reciben terapia hormonal.
- Echa un vistazo a la curva de supervivencia con
str().
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Weibull model
wbmod <- survreg(Surv(___, ___) ~ ___, data = GBSG2)
coef(wbmod)
# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
t_yes <- predict(___, type = "quantile", p = ___,
newdata = data.frame(horTh = ___))
# Take a look at survival curve
str(t_yes)