Ước lượng đường cong sống sót với survreg()
Giờ bạn sẽ ước lượng khả năng sống sót của bệnh nhân ung thư vú trong dữ liệu GBSG2 bằng mô hình Weibull.
Phân phối Weibull có hai tham số, quyết định dạng của đường cong sống sót.
Gói survival và dữ liệu GBSG2 đã được nạp sẵn cho bạn trong bài tập này.
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ột mô hình Weibull cho bệnh nhân ung thư vú.
- Tính đường cong sống sót ước lượng từ mô hình bằng hàm
predict()vớitype = "quantile". - Tạo một
data.framechứa các mốc thời gian và xác suất sống sót tương ứng. - Xem vài dòng đầu của kết quả bằng hàm
head().
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
wb <- ___(___(___, ___) ~ ___, data = ___)
# Retrieve survival curve from model probabilities
surv <- seq(.99, .01, by = -.01)
# Get time for each probability
t <- predict(___, type = ___, p = ___, newdata = data.frame(1))
# Create data frame with the information
surv_wb <- data.frame(time = ___, surv = ___)
# Look at first few lines of the result
head(___)