estou querendo pré-determinar contrastes ortogonais na anova para estimar suas somas de quadrado.
Entretanto, a ordem com que eu insiro os contrastes está diferenciando as somas de quadrado (acredito que seja função do ajuste sequencial)
Alguém sabe como determinar as somas de quadrado de 2 ou mais contrastes simultâneos usando a função aov() e lm().
exp <- structure(list(trat = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L,
2L, 3L, 3L, 3L, 3L), .Label = c("1", "2", "3"), class = "factor"),
y = c(35.1246524994234, 45.2412218795633, 38.6536395943239,
44.9890868406512, 46.8437044071867, 31.5689442940656, 43.7282923393773,
38.1832126839309, 33.8260517315041, 34.600238427063, 46.2089275518832,
24.3185159526907)), .Names = c("trat", "y"), row.names = c(NA,
-12L), class = "data.frame")
attach(exp)
m0 <- aov(y ~ trat)
anova(m0)
mat.c <- matrix(c(-1,0,1,2,-1,-1),nc=2)
c1 <- C(trat,c(-1,0,1),1)
c2 <- C(trat,c(2,-1,-1),1)
contrasts(trat) <- mat.c
m1 <- aov(y ~ c1 + trat) # Somente a soma de quadrados de c1 (78.47) correta !!!!!!!!!!!
anova(m1)
m2 <- aov(y ~ c2 + trat) # Somente a soma de quadrados de c2 (34.41) correta !!!!!!!!!!!
anova(m2)
Desde já agradeço por qualquer contribuição.
Att.
Tiago.
#################################################################
Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES
Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas
#################################################################