Sim Walmes esta é a saída que desejo. No entanto, observe que se eu alterar a ordem dos contrastes as somas de quadrado alteram. Exemplo que você me enviou:mat.c [,1] [,2][1,] -1 2[2,] 0 -1[3,] 1 -1 Df Sum Sq Mean Sq F value Pr(>F)trat 2 91.5 45.75 0.915 0.435 trat: C.orto.1 1 78.5 78.47 1.569 0.242 # soma de quadrados correta trat: C.orto.2 1 13.0 13.03 0.261 0.622 # esperava para este contrate (2 - 1 -1) o valor de 34.4Residuals 9 450.0 50.00 Agora invertendo a ordem de declaração dos contrastes: mat.c [,1] [,2][1,] 2 -1[2,] -1 0[3,] -1 1 Df Sum Sq Mean Sq F value Pr(>F)trat 2 91.5 45.75 0.915 0.435 trat: C.orto.1 1 34.4 34.41 0.688 0.428 # soma de quadrados correta trat: C.orto.2 1 57.1 57.09 1.142 0.313 # esperava para este contrate (-1 + 0 + 1) o valor de 78.5Residuals 9 450.0 50.00 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 ################################################################# Date: Thu, 16 Apr 2015 16:06:33 -0300 From: walmeszeviani@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Partição de somas de quadrado usando contrastes ortogonais Seria isso? str(exp) contrasts(exp$trat) mat.c <- matrix(c(-1,0,1,2,-1,-1), nc=2) contrasts(exp$trat) <- mat.c m0 <- aov(y~trat, data=exp) summary(m0, split=list("trat"=list("C.orto.1"=1, "C.orto.2"=2))) À disposição. Walmes. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.