[R-br] teste de médias para interação

[Daniele Portela de oliveira]

Benilton eu usei esse modelo e os dados estão em anexo. 
 
# Poisson dados balanceados
mod.p1 <- glm(Nasc ~ avali*trat,family=poisson(link = "log"),contrasts=list(avali="contr.sum",trat="contr.sum"))
an.p1 <- anova(mod.p1,test="Chisq")
an.p1
 
O tratamento (trat), o período de avaliaçao (avali) e a interação foram significativos, eu gostaria de fazer
o teste para a interação.
O teste de médias para trat e avali eu usei o multcomp(abaixo) mas para a interação o mesmo teste não deu certo.
 
#Teste de médias para Avali
library(multcomp)
glht.mod1 <- glht(mod.p1, linfct=mcp(avali="Tukey"), alternative = "two.sided" )
summary(glht.mod1,test=Chisqtest())
confint(glht.mod1)
cld(glht.mod1)

Eu vou tentar montar os contrastes para fazer o teste.
 


Daniele Portela de Oliveira
Mestranda Produção Animal - UEM
Cel: 44 9949 3525
       67 9617 0965
 

De: Benilton Carvalho <beniltoncarvalho em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Daniele Portela de oliveira <dane_portela em yahoo.com.br>
Enviadas: Terça-feira, 1 de Novembro de 2011 12:58
Assunto: Re: [R-br] teste de médias para interação

Se vc se refere a "teste de hipoteses envolvendo termos de interacao
de um GLM", voce deve construir a matriz de contrastes de interesse e,
entao, os exemplos da ajuda do glht() podem ser aplicados diretamente.

Vc e' bem-vinda a dar um exemplo reproduzivel do seu caso e,
certamente, tera' respostas mais precisas que a minha (ne, Walmes?).

b
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111101/f4b7cbbc/attachment.html>
-------------- Próxima Parte ----------
Um texto embutido e sem conjunto de caracteres especificado foi limpo...
Nome: dados.txt
Url: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111101/f4b7cbbc/attachment.txt>