web-dev-qa-db-ja.com

BIC基準を使用した段階的線形モデルの実行

AICではなくBIC基準を使用するように段階的線形モデルを設定することは可能ですか?

私はこれを試してきましたが、それでもBICではなくAIC値を使用して各ステップを計算します

null = lm(data[,1] ~ 1)
full = lm(data[,1] ~ age + bmi + gender + group)
step(null, scope = list(lower=null,upper=full),
     direction="both", criterion = "BIC")
12
user2846211

引数k=log(n)step関数に追加します(nモデル行列のサンプル数)

?stepから:

引数:
.。

kペナルティに使用される自由度の数の倍数。 k = 2のみが本物のAICを与えます。 k = log(n)は、BICまたはSBCと呼ばれることもあります。

16
rcs