2013-04-10 129 views
0

我正在使用svyglm模型,我需要獲得「AIC」來比較不同模型。svyglm中的擬合模型

我的問題是,如果我用我的模型:出現

Modelo12=svyglm(formula = Asiste ~ sexo + E27 + JovenActivo + hijos + 
    +  jefe + LN_YSVL_sin_joven_prom + aniosed + climaeducativo + 
    +  icv2 + TV + Computadora + Telefono + Cable + Calefon + 
    +  DVD + Microhondas + Aire + Auto_o_moto + Secadora + Lavavajillas + 
    +  Refrigerador + Actividad_del_Jefe + Hacinamiento, family = quasibinomial(link 
    =  
    "probit"), data = Personas.con.muestra, design=diseno_personas_14_17) 

當我寫總結(Modelo12)AIC犯規。

在另一方面,我不能老是用階梯,這是錯誤:

stepwise(Modelo12) 

    Direction: backward/forward forward/backward backward forward 
    Criterion: BIC 

    Error en extractAIC.svyglm(fit, scale, k = k, ...) : 
    svyglm not fitted by maximum likelihood 

謝謝!

納塔利婭

+0

嗨natalia,歡迎來到SO。一個快速的谷歌搜索[來自'調查'包作者的這個筆記](https://stat.ethz.ch/pipermail/r-help/2012-July/319508.html) – 2013-04-10 21:55:45

回答

0

AIC是一個基於似然度量。 svyglm不使用最大似然擬合,因爲錯誤和?svyglm告訴你。在那裏你也可以得到使用regTermTest的建議,它提供Wald測試和調整的Rao-Scott似然比測試。

+0

感謝您的回答!是的,我認爲svyglm使用Horvitz Thomson估計量可能是? 我看到regTermTest會計算F,但是我怎麼能自動比較模型呢? 謝謝 – Natuk 2013-04-17 15:58:06