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@gmail.com>
Para: r-br@listas.c3sl.ufpr.br; Daniele Portela de oliveira <dane_portela@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