Add the argument k=log(n)
to the step
function (n
number of samples in the model matrix)
From ?step
:
Arguments:
...k the multiple of the number of degrees of freedom used for the penalty. Only k = 2 gives the genuine AIC; k = log(n) is sometimes referred to as BIC or SBC.