Esecuzione di un modello lineare stepwise con criterio BIC

Question

È possibile impostare un modello lineare stepwise per utilizzare i criteri BIC anziché AIC?

Ho provato questo ma calcola ancora ogni passo utilizzando valori AIC anziché 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

Aggiungere l'argomento k=log(n) al step funzione (n numero di campioni nella matrice modello)

Da ?step:

Argomenti:
...

k il multiplo del numero di gradi di libertà utilizzati per la sanzione. Solo k = 2 fornisce l'AIC originale; k = log (n) è talvolta denominato come BIC o SBC.