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