Выполнение ступенчатой ​​линейной модели с критерием BIC

Можно ли установить ступенчатую линейную модель для использования критериев 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")

Ответ 1

Добавьте аргумент k=log(n) в функцию step (n количество выборок в матрице модели)

От ?step:

Аргументы:
...

k кратное количеству степеней свободы, используемых для наказания. Только k = 2 дает подлинную AIC; k = log (n) иногда упоминается как BIC или SBC.