Re: [R-es] obtener residuos de una Anova con biblioteca CAR

[rubenfcasal]

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
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] obtener residuos de una Anova con biblioteca CAR

[eric]

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 
*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@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.


___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es