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") 

、コックス回帰の要約結果からコピーして貼り付けられたデータでは女性は危険率があることを暗示し、低い男性よりもを。18 $\exp(-0.1929)=0.8245$summary(reg_fit)$18\%$

関数basehaz（前の回答から）は、ハザード関数（レート）ではなく、累積ハザードを提供します。その問題はハザード機能に関するものだと思います。ハザード関数を推定するには、平滑化のタイプを指定する必要があります（密度推定など）。MuhazRパッケージは、一つのサンプルデータのためにこれを行うことができます。Coxモデルのベースラインハザードに対してこれを実行できる関数を知りません。これも必要です。自分で累積ハザードをなめらかにする必要があると思います。