使用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.