Cox PH 모델에서 예측 위험률을 계산하는 방법은 무엇입니까?

다음 Cox PH 모델이 있습니다.

(시간, 이벤트) ~ X + Y + Z

내가 예측 위험 좀하고 싶습니다 요금 (내가 위험 비율에 대해 이야기하고 하지 의 특정 값을 주어 위험 비율) X, Y, Z. muhaz R 패키지 가 관측 된 위험률을 계산할 수 있다는 것을 알고 있지만 예측 모델에 관심이 있습니다.

R에서 이것을 할 수있는 방법이 있습니까?

기준 위험률을 얻는 R 생존 패키지의 기능은 basehaz입니다.

그런 다음 다양한 에 대해 계수를 곱하여 찾은 계수가 지정된 특정 위험률을 구해야합니다.$e^{\beta}$

간단한 예가 도움이 될 수 있습니다.

library(survival) #survival analysis
library(eha) #used for data 
data(oldmort) #create the data

# Create surv data set
mort <- Surv(time=oldmort$enter,time2=oldmort$exit,event=oldmort$event)

reg_fit <- coxph(formula=mort~oldmort$sex)
summary(reg_fit)

# Now get baseline curve
baseline <- basehaz(reg_fit)

# Draw baseline hazard (that's male)
plot(baseline$time, baseline$hazard, type='l',main="Hazard rates") 

# Draw female hazard
lines(baseline$time, exp(-0.1929)*baseline$hazard, col="blue") 

복사 콕스 회귀 분석의 요약 결과에서 붙여, , 암시 데이터의 여성은 위험 비율이 가지고 남성보다 낮다.$\exp(-0.1929)=0.8245$summary(reg_fit)$18\%$

함수 basehaz(이전 답변에서)를 제공 누적 위험 하지 위험 함수 (속도). 그 질문은 위험 기능에 관한 것이라고 생각합니다. 위험 함수를 추정하려면 평활화 유형 (밀도 추정과 같은)을 지정해야합니다. MuhazR 패키지는 하나의 샘플 데이터에 대해이 작업을 수행 할 수 있습니다. Cox 모델의 기준 위험에 대해이를 수행 할 수있는 기능을 알지 못합니다. 나도 이것이 필요합니다. 나는 스스로 누적 위험을 완화해야한다고 생각합니다.