survreg()로 생존곡선 추정하기
이제 GBSG2 데이터의 유방암 환자 생존을 Weibull 모형으로 추정해 보겠습니다.
Weibull 분포에는 생존곡선의 형태를 결정하는 두 개의 모수가 있습니다.
이 연습 문제에서는 survival 패키지와 GBSG2 데이터가 이미 로드되어 있습니다.
이 연습은 강의의 일부입니다
R로 배우는 생존 분석
연습 안내
- 유방암 환자에 대해 Weibull 모형을 추정하세요.
type = "quantile"옵션과 함께predict()함수를 사용해 모형에서 추정된 생존곡선을 계산하세요.- 시간 지점과 해당 생존확률을 포함하는
data.frame을 만드세요. head()함수를 사용해 결과의 앞부분을 확인하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# 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(___)