R 重回帰分析におけるstep関数の項数を固定したい

Question

前提・実現したいこと

説明変数が50個ある重回帰分析モデルをステップワイズ法にてモデルの最適化を行っているのですが，実装の都合上で説明変数が8個のモデルしか実装できないので，指定した説明変数の個数での変数選択を行いたいのです。

試したこと

関数step()を利用していますが，それでは説明変数の個数が指定できないので困っております。

補足情報（FW/ツールのバージョンなど）

R 4.1.1

Answer 1

あまり詳しくないのですが、私が理解した範囲ではステップワイズ法とは「変数の数を増やしたり減らしたりしながら」最適なモデルを探す手法。最初から変数の数が決まっているという状況にはそもそも合致していないように思います。

別の方法として、パッケージ bestglmにbestglm 関数というのがあり、こちらは総当たり法をサポートしています。マニュアルを見る限り、nvmaxというオプションで説明変数の個数の最大数を指定できるようですから、これを試してみるのはどうでしょうか。

bestglm function - RDocumentation

Comments

[MZtech]

KojiDoi様，回答ありがとうございました。
無事，目的のモデル実装ができました。非常に感謝しております。