Partição de somas de quadrado usando contrastes ortogonais
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 #################################################################
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.
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.
Walmes desconsidere a ultima mensagem. Não tinha me atentado para a falta de ortogonalidade entre os contrastes (somatório de ai*bi =0). 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 ################################################################# From: tiagosouzamarcal@hotmail.com To: r-br@listas.c3sl.ufpr.br Subject: RE: [R-br] Partição de somas de quadrado usando contrastes ortogonais Date: Thu, 16 Apr 2015 22:28:44 +0300 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.
participantes (2)
-
Tiago Souza Marçal -
walmes .