r - 使用 BIC 准则运行逐步线性模型

Question

是否可以设置逐步线性模型以使用 BIC 标准而不是 AIC？

我一直在尝试这个，但它仍然使用 AIC 值而不是 BIC 计算每个步骤

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

Answer 1

将参数添加k=log(n)到step函数（n模型矩阵中的样本数）

来自?step：

论据：
...

k用于惩罚的自由度数的倍数。只有 k = 2 给出真正的 AIC；k = log(n) 有时称为 BIC 或 SBC。