Calcola il modello di Weibull
Per un modello di Weibull con covariate, possiamo calcolare la curva di sopravvivenza proprio come abbiamo fatto per il modello di Weibull senza covariate. L’unica cosa da fare è specificare i valori delle covariate per una data curva di sopravvivenza nella funzione predict(). Questo si può fare con l’argomento newdata.
Questo esercizio fa parte del corso
Analisi di sopravvivenza in R
Istruzioni dell'esercizio
- Calcola un modello di Weibull per il dataset
GBSG2con la covariatahorThper analizzare l’effetto della terapia ormonale sulla sopravvivenza dei pazienti. - Calcola la curva di sopravvivenza per i pazienti che ricevono la terapia ormonale.
- Dai un’occhiata alla curva di sopravvivenza con
str().
esercizio interattivo pratico
Prova questo esercizio completando questo codice di esempio.
# 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)