Hitung model Weibull
Untuk model Weibull dengan kovariat, kita dapat menghitung kurva kelangsungan hidup seperti yang kita lakukan pada model Weibull tanpa kovariat. Satu-satunya hal yang perlu kita lakukan adalah menentukan nilai kovariat untuk kurva kelangsungan hidup tertentu dalam fungsi predict(). Ini dapat dilakukan dengan argumen newdata.
Latihan ini merupakan bagian dari kursus
Analisis Survival di R
Instruksi latihan
- Hitung model Weibull untuk himpunan data
GBSG2dengan kovariathorThuntuk menganalisis pengaruh terapi hormonal terhadap kelangsungan hidup pasien. - Hitung kurva kelangsungan hidup untuk pasien yang menerima terapi hormonal.
- Lihat struktur kurva kelangsungan hidup dengan
str().
Latihan interaktif langsung praktik
Cobalah latihan ini dengan melengkapi kode contoh ini.
# 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)