Hola Eric,
Según entiendo, el Anova SSIII es una forma de descomponer la suma
de cuadrados para contrastar efectos de los factores (principalmente
para tratar de evitar complicaciones a los usuarios...). El ajuste del
modelo (las estimaciones de los parámetros) no dependen de este
contraste (se hacen a priori). Por tato para mí solo tiene sentido
calcular residuos del modelo ajustado, i.e. de tu objeto inf.ninf.glm
(que por cierto podrías haberlo creado con lm en lugar de glm...).
Un saludo, Rubén.
El 20/02/2016 a las 22:45, eric escribió:
Estimada Comunidad, hice un Anova (SS III) usando la biblioteca CAR y
necesito obtener los residuos del ajuste ... el argumento de Anova fue
un modelo ajustado con glm(), como indica el siguiente codigo:
inf.ninf.glm <- glm(CH ~ pa.+lam+ps.+age+area+expo+inf+P+ppacum1mes,
data=hum[hum$CH < 45,])
inf.ninf.Aov <- Anova(inf.ninf.aov, type=3, error.estimate="pearson")
mi pregunta ahora, es como extraigo los residuos desde inf.ninf.Aov,
para no hacer calculos manuales estimando y calculando diferencias y
todo eso.
He intentado con residuals(), residuals.glm() y otros, pero no me
resulta y obtengo el siguiente error:
Error in residuals.glm(inf.ninf.Aov, type = "pearson") :
could not find function "mu.eta"
Alguna pista ?
Muchas gracias,
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
