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")
引数k=log(n)
をstep
関数に追加します(n
モデル行列のサンプル数)
?step
から:
引数:
.。kペナルティに使用される自由度の数の倍数。 k = 2のみが本物のAICを与えます。 k = log(n)は、BICまたはSBCと呼ばれることもあります。