[R-br] Partição de somas de quadrado usando contrastes ortogonais

[Tiago Souza Marçal]

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 ################################################################# 
 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150416/8fe1c200/attachment.html>