2013-05-19 61 views
0

我有這個模型,我需要計算「AIC」來進行模型選擇,我的問題是當我運行模型「AIC」顯示爲「NA」時。如何計算svyglm模型中的「AIC」?

MODEL= svyglm(formula = Asiste ~ afro + sexo + E27 + JovenActivo + hijos + 
jefe + LN_YSVL_sin_joven_prom + aniosed + climaeducativo + 
Mdeo + icv2 + TV + Computadora + Telefono + Internet + Cable + 
Calefon + DVD + Microhondas + Aire + Auto_o_moto + Secadora + 
Madre_ausente + Internet + TV + Lavavajillas + Refrigerador + 
climaeducativo + Actividad_del_Jefe + Hacinamiento, family = quasibinomial(link =  
"logit"), 
data = Personas.con.muestra, design = diseño_personas_14_17, 
subset = (Stratum != 0)) 

MODEL$aic 
[1] NA 

有什麼建議嗎?

謝謝!納塔利婭

+0

搜索了'survey'包文檔http://cran.r-project.org/web/packages/survey/survey.pdf,我看沒有提及「AIC」或「Akaike」。也許你想運行'AIC(MODEL)'? – Frank

+2

您認爲擬合準則能產生AIC估計值的基礎是什麼? (即使準家庭的AIC定義明確,您提出的程序仍然是統計上可疑的。)如果您懷疑我的話,請搜索R-help檔案。請務必特別注意Lumley和Harrell的帖子。 –

+5

@Frank。 'svyglm'將會有一個名爲「aic」的列表元素,但是它只針對普通家庭計算,對於準家庭來說總是NA。 –

回答

-1

你可以試試這個:

summary(MODEL)$aic