Como calcular as taxas de risco previstas de um modelo Cox PH?

Eu tenho o seguinte modelo Cox PH:

(Hora, Evento) ~ X + Y + Z

Gostaria de obter os perigos previstos taxas (eu estou falando sobre taxas de risco NÃO taxas de risco) dado valores específicos de X, Y, Z. Sei que o pacote muhaz R pode calcular as taxas de risco observadas, mas estou interessado no modelo previsto.

Existe uma maneira de fazer isso no R?

A função no pacote de sobrevivência R para obter a taxa básica de risco é basehaz.

$e^{\beta}$

Um exemplo simples pode ajudar:

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\%$

A função basehaz(da resposta anterior) fornece o risco cumulativo , não a função de risco (a taxa). Eu acredito que essa pergunta era sobre a função de risco. A estimativa da função de risco exigiria a especificação do tipo de suavização (como na estimativa de densidade). O Muhazpacote R pode fazer isso para uma amostra de dados. Não conheço uma função que possa fazer isso para o risco de linha de base no modelo Cox. Eu também preciso disso. Acho que precisarei diminuir o risco cumulativo por conta própria.