r - R で Cox 比例ハザード モデルのすべてのサブセット変数選択を実行できますか?

Question

データの上位の Cox 比例ハザード モデルを選択するときに、実際にはそうではないにしても、プログラム R の jumps パッケージの regsubsets に似た関数を使用しようとしています。これは可能ですか？もしそうなら、関数はすでに存在しますか？

Answer 1

あなたはすでに次のことをよく知っていると思います... 「トップモデル」の基準としてAICを使用している場合、これは妥当な出発点になります。

library(survival)
data(colon)
c1 <- coxph(Surv(time=time, event=status) ~
     as.factor(extent) + age + sex, data=colon)
step(c1)

欠損値がある場合は注意してください ( NA)。もちろん、この方法では見つからないより良いモデルが存在する可能性がありますが、潜在的な予測変数が少数であれば、それを見逃す可能性はほとんどありません。情報に基づいた意見がより信頼できる数値手法の使用に関する上記の警告 (@DWin に感謝)。