Ước lượng mô hình Weibull
Với mô hình Weibull có biến đồng biến, bạn có thể tính đường cong sống sót giống như khi làm với mô hình Weibull không có biến đồng biến. Điều duy nhất cần làm là chỉ định giá trị của các biến đồng biến cho đường cong sống sót đó trong hàm predict(). Bạn thực hiện điều này bằng đối số newdata.
Bài tập này là một phần của khóa học
Phân tích Sống sót trong R
Hướng dẫn bài tập
- Ước lượng mô hình Weibull cho tập dữ liệu
GBSG2với biến đồng biếnhorThđể phân tích tác động của liệu pháp hormone lên thời gian sống của bệnh nhân. - Tính đường cong sống sót cho những bệnh nhân nhận liệu pháp hormone.
- Xem cấu trúc của đường cong sống sót bằng
str().
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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)