Tengo un dataframe
con varios valores faltantes y quiero realizar un análisis de la varianza (ANOVA) para comparar dos modelos logit binomial:
- Modelo A: contiene un set de variables.
- Modelo B: Contiene las mismas variables que el Modelo A más 3 variables de estudio.
Importamos los datos:
Modelo A:
modelo_logit <- glm(SAP ~ sexo + edad + peso + niv_est + enf_cron + sit_lab +
frec_act_fis + ingreso_eq + GHQ_12,
data = datos_modelo, family = binomial(link = "logit"),na.action = "na.omit")
Modelo B:
modelo_logit_viv <- glm( SAP ~ sexo + edad + niv_est + enf_cron + sit_lab +
frec_act_fis + ingreso_eq + GHQ_12 +
n_dormitorios + cont_indus + delincuencia, # variables de estudio
data = datos_modelo, family = binomial(link = "logit"),na.action = "na.omit")
Para poder realizar un ANOVA ejecuto: anova(modelo_logit,modelo_logit_viv)
Y obtengo el siguiente error:
Error in anova.glmlist(c(list(object), dotargs), dispersion = dispersion, : models were not all fitted to the same size of dataset
Ambos modelos deben ser ajustados por el mismo por el mismo set de datos, pero como hay varios valores faltantes, el modelo B tiene más NA's
que el modelo A (ya que el modelo B contiene más variables que incrementan el número de observaciones que han de eliminarse con respecto al modelo A).
Mi pregunta entonces es: ¿Cómo puedo obtener el dataframe que se emplea en el modelo B (modelo_logit_viv
) para poder estimar el modelo A (modelo_logit
) y tener así el mismo dataframe para estimar ambos modelos y efectuar a continuación el ANOVA? Tiene que haber algún elemento dentro del objeto generado por glm(*)
que contenga el dataframe
que se ha empleado para estimar una vez eliminado los NA's
, pero no lo encuentro.