Estoy corriendo el siguiente modelo lineal general en R:
model1<-lm(leche ~ MT1 + ednp + cnv3 + poly(DIM, 3) + HTD, data=misdatos)
MT1
, ednp
, cnv3
y HTD
son factores, los he predefinido como tal; mientras que DIM
es una variable continua que utilizo como covariable en el modelo. Ejecutando el modelo de esta manera, mi anova considera todos los efectos
Analysis of Variance Table
Response: leche
Df Sum Sq Mean Sq F value Pr(>F)
MT1 1 2788946 2788946 67.049 2.703e-16 ***
ednp 9 87481407 9720156 233.680 < 2.2e-16 ***
cnv3 2 4457347 2228674 53.579 < 2.2e-16 ***
poly(DIM, 3) 3 947794411 315931470 7595.244 < 2.2e-16 ***
HTD 1012 762410521 753370 18.112 < 2.2e-16 ***
Residuals 57613 2396468010 41596
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Sin embargo, cuando cambio el orden de mis efectos a
model2<-lm(leche ~ HTD + MT1 + cnv3 + poly(DIM, 3) + ednp, data=misdatos)
el resultado omite la salida para MT1
, lo que creo que carece de sentido porque en teoría el orden no debería alterar el resultado, ¿no?
Analysis of Variance Table
Response: leche
Df Sum Sq Mean Sq F value Pr(>F)
HTD 1013 1350638499 1333306 32.054 < 2.2e-16 ***
ednp 9 45711618 5079069 122.105 < 2.2e-16 ***
cnv3 2 1413259 706630 16.988 4.211e-08 ***
poly(DIM, 3) 3 407169256 135723085 3262.891 < 2.2e-16 ***
Residuals 57613 2396468010 41596
¿Qué puede estar pasando? Esto ya me crea inseguridad en los resultados que obtengo y no se me ocurre que puede estar pasando. ¿Alguna idea?