Boa tarde pessoal, 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(). Segue a baixo um cmr com minhas tentativas iniciais. CMR 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 #################################################################