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")