予測にはcoxモデルを使用しています。
fit <- coxph(Surv(time, status) ~ age + ph.ecog + strata(inst), lung)
predict(fit,type="expected")
predict.coxph doc では、予想されるタイプは「共変量とフォローアップ時間が与えられたイベントの予想数」であると書かれています。ここでフォローアップ時間とはどういう意味ですか? どうすれば変更できますか?
予測にはcoxモデルを使用しています。
fit <- coxph(Surv(time, status) ~ age + ph.ecog + strata(inst), lung)
predict(fit,type="expected")
predict.coxph doc では、予想されるタイプは「共変量とフォローアップ時間が与えられたイベントの予想数」であると書かれています。ここでフォローアップ時間とはどういう意味ですか? どうすれば変更できますか?
これは、このコホートの平均生存時間よりも長くなるように選択された単一の時間が選択された、患者機能の特定のセットの「期待」の例ですが、「肺」データセットのすべての一意の施設コードを使用しています。
(pr.fit.new <- predict(fit, newdata=data.frame(age=50, time=365, status=0,
ph.ecog=2,inst=unique(lung$inst[!is.na(lung$inst)]) ),
type="expected"))
[1] 1.1886094 1.7679429 1.6651195 1.3117278 0.7743925 0.8501810
[7] 1.9686954 1.2560981 2.5845069 0.6513932 1.1241870 1.1515346
[13] 1.1616482 1.9071766 4.2468797 0.4921904 0.5875339 0.5647238
層には著しく異なる「ベースラインハザード」があり、共変量と観察期間が同じである仮説上の人物の「期待値」に重大な違いが生じます。