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 adalah bagian dari kursus
Analisis Survival di R
Petunjuk 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 praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# 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)