相互作用を含むサバイバル コックス モデルを作成しましたcovariate * time
(不均衡が検出されました)。モデルから生存予測を最も簡単に得るにはどうすればよいか、今疑問に思っています。
私のモデルが指定されました:
coxph(formula = Surv(event_time_mod, event_indicator_mod) ~ Sex +
ageC + HHcat_alt + Main_Branch + Acute_seizure + TreatmentType_binary +
ICH + IVH_dummy + IVH_dummy:log(event_time_mod)
そして今、私は予測を行っている変数の組み合わせを使用survfit
して提供する予測を得ることを望んでいました:new.data
survfit(cox, new.data=new)
event_time_mod
モデルの右側にあるように、 に渡される新しいデータ フレームでそれを指定する必要がありますsurvfit
。これevent_time
は、予測の個々の時間に設定する必要があります。event_time_mod
の正確な時刻を 指定する簡単な方法はありますsurvfit
か? または、モデルから予測を達成するための他のオプションはありますか?
もちろん、予測とevent_time_mod
正しい値への設定に明確な時間があるので、新しいデータ フレームに同じ数の行を作成することもできますが、それは非常に面倒で、もっと良い方法があるに違いないと思いました。