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

[Tiago Souza Marçal]

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