博士号取得のために行っている生存分析のハザード関数をプロットして、2 つの異なる条件のハザード率を比較しようとしています。
予測変数の両方のレベルの平滑化されたハザード線の信頼区間を取得するために、コードを意図したとおりに機能させる方法を見つけることができません (参照用に、図 4、ページ 7)。
参照用にコードを追加しています：
fitt<-bshazard(Surv(time,event) ~ session.type,data=data,lambda=10,nbin=60)
plot(fitt,overall=FALSE, col=1, conf.int = TRUE)
関数「overall=FALSE」は、プロットから結果を推定する必要がある信頼区間を含まない 2 つの平滑化されたハザード曲線を提供します。コードから取得したプロットの画像を次に示します。
各時間間隔の値を知るために、時間表のハザード率 (上限と下限の信頼区間) を取得する方法を誰かが知っていれば、大いに役立ちます。
助けてくれる人に感謝します！