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, -- Forest Engineer Master in Environmental and Natural Resource Economics Ph.D. student in Sciences of Natural Resources at La Frontera University Member in AguaDeTemu2030, citizen movement for Temuco with green city standards for living Nota: Las tildes se han omitido para asegurar compatibilidad con algunos lectores de correo.
Javier Marcuzzi
2016-Feb-22 11:41 UTC
[R-es] obtener residuos de una Anova con biblioteca CAR
Estimado Eric Sirve esto http://www.uv.es/ayala/docencia/cda/notaR/notaR025.pdf? Javier Rubén Marcuzzi De: eric Enviado: sábado, 20 de febrero de 2016 18:45 Para: Lista R Asunto: [R-es] obtener residuos de una Anova con biblioteca CAR 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, -- Forest Engineer Master in Environmental and Natural Resource Economics Ph.D. student in Sciences of Natural Resources at La Frontera University Member in AguaDeTemu2030, citizen movement for Temuco with green city standards for living Nota: Las tildes se han omitido para asegurar compatibilidad con algunos lectores de correo. _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es [[alternative HTML version deleted]]
Voy a probar Javier, y contar si funciona, muchas gracias, eric. On 02/22/2016 08:41 AM, Javier Marcuzzi wrote:> Estimado Eric > > Sirve esto http://www.uv.es/ayala/docencia/cda/notaR/notaR025.pdf? > > Javier Rubén Marcuzzi > > *De: *eric <mailto:ericconchamunoz en gmail.com> > *Enviado: *sábado, 20 de febrero de 2016 18:45 > *Para: *Lista R <mailto:r-help-es en r-project.org> > *Asunto: *[R-es] obtener residuos de una Anova con biblioteca CAR > > 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, > > -- > > Forest Engineer > > Master in Environmental and Natural Resource Economics > > Ph.D. student in Sciences of Natural Resources at La Frontera University > > Member in AguaDeTemu2030, citizen movement for Temuco with green city > > standards for living > > Nota: Las tildes se han omitido para asegurar compatibilidad con algunos > > lectores de correo. > > _______________________________________________ > > R-help-es mailing list > > R-help-es en r-project.org > > https://stat.ethz.ch/mailman/listinfo/r-help-es >-- Forest Engineer Master in Environmental and Natural Resource Economics Ph.D. student in Sciences of Natural Resources at La Frontera University Member in AguaDeTemu2030, citizen movement for Temuco with green city standards for living Nota: Las tildes se han omitido para asegurar compatibilidad con algunos lectores de correo.
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,
>
>
>