R语言 如何在svyglm模型中计算“AIC”?

nnsrf1az  于 2023-09-27  发布在  其他
关注(0)|答案(2)|浏览(180)

我有这个模型,我需要计算“AIC”来进行模型选择,我的问题是,当我运行模型时,“AIC”显示为“NA”。

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

有什么建议吗?
谢谢!!娜塔莉亚

sg24os4d

sg24os4d1#

现在在survey包中有一个AIC.svyglm函数,基于Lumley & Scott开发的新方法(2015)。基本上,$2p$惩罚需要取决于权重(更高的权重意味着更高的惩罚)和由于聚类而增加的方差(更多的不确定性意味着更高的惩罚)

e3bfsja2

e3bfsja22#

你可以试试这个:

summary(MODEL)$aic

相关问题