2013-10-16 26 views
10

是否可以將逐步線性模型設置爲使用BIC標準而不是AIC?使用BIC標準運行逐步線性模型

我一直想這一點,但它仍然使用計算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") 

回答

13

參數k=log(n)添加到step功能(在模型矩陣n採樣數)

每一步

?step

參數:
...

k用於懲罰的自由度數的倍數。 只有k = 2給出真正的AIC;有時將k = log(n)稱爲 作爲BIC或SBC。

+3

'n'是必須知道的觀測值的數量,而'k'是每個參數使用的懲罰的數值。 – 2014-04-01 11:17:55