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

Question

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

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に感謝します）は、情報に基づいた意見がより信頼できる数値方法を使用することについて話します。