[R-br] multcomp/doBy comparação multipla com interação

Fernando A. Souza nandodesouza em gmail.com
Segunda Dezembro 3 16:10:32 BRST 2012


Caros colegas
Estou fazendo uma análise de um experimento DIC com arranjo fatorial 
(Gest e Manej). A resposta medida é o peso vivo. A análise de variância 
mostrou que há efeito de interação entre os dois fatores avaliados. Eu 
gostaria de saber onde está esta diferenças. A minha dúvida é o seguinte 
como fazer a matrix de contrastes para a interação? estou utilizando o 
pacote doBy.

A outra dúvida é a seguinte: o R utiliza como padrão o contraste 
tratamento o qual compara todos os tratamento como o grupo controle( em 
meu caso Gest=0 e Manej=1) como faço para alterar esta forma de 
contraste, para um que compare todos os tratamentos entre si ?

>dput(abate[,1:5])
structure(list(Animal = c(1L, 2L, 3L, 4L, 5L, 7L, 8L, 9L, 12L,
14L, 15L, 17L, 18L, 19L, 21L, 22L, 23L, 25L, 26L, 28L, 30L, 32L,
34L, 35L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 46L, 47L, 48L,
49L, 50L, 53L, NA, NA), Gest = structure(c(4L, 2L, 3L, 2L, 2L,
4L, 3L, 4L, 2L, 2L, 3L, 4L, 1L, 4L, 2L, 1L, 4L, 1L, 2L, 4L, 3L,
3L, 1L, 4L, 3L, 1L, 2L, 3L, 1L, 1L, 3L, 1L, 2L, 3L, 1L, 2L, 1L,
1L, NA, NA), .Label = c("0", "100", "130", "140"), class = "factor"),
     Manej = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L,
     2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L,
     1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, NA, NA
     ), .Label = c("1", "2"), class = "factor"), rep = c(2L, 3L,
     1L, 2L, 4L, 3L, 1L, 2L, 2L, 5L, 4L, 3L, 1L, 4L, 1L, 2L, 4L,
     4L, 3L, 1L, 3L, 5L, 6L, 1L, 2L, 5L, 1L, 2L, 4L, 3L, 3L, 3L,
     4L, 4L, 5L, 5L, 1L, 2L, NA, NA), PV = c(49000L, 40800L, 51800L,
     41900L, 35000L, 40500L, 35400L, 42300L, 36600L, 36600L, 38400L,
     46400L, 26900L, 50200L, 35500L, 35500L, 48400L, 40500L, 32500L,
     42800L, 36500L, 50000L, 31200L, 33500L, 46500L, 40900L, 36500L,
     31200L, 31500L, 35500L, 38700L, 30000L, 37900L, 42700L, 35500L,
     42000L, 34500L, 36700L, NA, NA)), .Names = c("Animal", "Gest",
"Manej", "rep", "PV"), class = "data.frame", row.names = c(NA,
-40L))

abate$Gest<-factor(abate$Gest)
abate$Manej<-factor(abate$Manej)
C<-lm(PV~Gest*Manej,abate)
Anova(C,type="III")


-- 
Fernando Antônio de Souza
Zootecnista, Msc. Nutricão Animal, UFMG
lattes: http://lattes.cnpq.br/6519538815038307

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121203/6444c0b4/attachment.html>


Mais detalhes sobre a lista de discussão R-br