<html><body><div style="color:#000; background-color:#fff; font-family:Courier New, courier, monaco, monospace, sans-serif;font-size:10pt"><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 10pt;">Boa noite a todos,</div><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 10pt;"><br></div><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;">Tenho um conjunto de dados em arranjo fatorial que não atendem as pressuposições de normalidade e homocesdacidade de variâncias. Para contornar esse problema, optei por fazer a análise do mesmo por GLM; no entanto,como houve efeito significativo da interação, fiquei na dúvida de como proceder para fazer o desdobramento e fazer comparações múltiplas entre os tratamentos.</div><div style="font-family: 'Courier New',
courier, monaco, monospace, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><br></div><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;">Pesquisando na internet, encontrei um post que recomendava fazê-lo por meio de constrastes utilizando o pacote multcomp. Assim o fiz (como vcs podem conferir pelo código abaixo). No entanto, como migrei recentemente para o R, surgiu a dúvida se isso que fiz está correto ou não. Assim, gostaria da opinião de vcs sobre o código ou se tem uma outra forma de fazê-lo.</div><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><br></div><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 13px;
color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><br></div><div style="background-color: transparent;"><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">da <- expand.grid(rep=1:4, isol=factor(1:6), cult=factor(1:2))</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">da$y <- c(4.5, 25.5, 8.5, 12, 123, 19.5, 111.5, 83.5, 10, 4, 4.5,</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"> 90.5, 46.5, 61.5, 20.5, 33.5, 47.5, 39, 131.5, 68.5, 4,</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"> 10, 13.5, 10, 9, 0, 0, 0, 2, 0, 0, 0, 0, 12.5, 4.5,
6.5,</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"> 1, 0, 0, 1, 0, 0, 0.5, 0, 2, 2, 0.5, 0.5)</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"><br></font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">g0 <- glm(y~isol + cult + isol:cult, family=poisson, data=da)</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">summary(g0)</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">anova(g0, test="Chisq")</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif"
size="2"><br></font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">isolcult <- interaction(da$isol, da$cult, drop=TRUE)</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">m0 <- glm(da$y~isolcult, family=poisson)</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">require(multcomp)</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">glht.m0 <- glht(m0, mcp(isolcult = "Tukey"))</font></div><div style="background-color: transparent;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2">summary(glht.m0)</font></div><div style="background-color: transparent; color: rgb(0, 0, 0); font-size: 13px; font-family: 'Courier New',
courier, monaco, monospace, sans-serif; font-style: normal;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"><br></font></div><div style="background-color: transparent; color: rgb(0, 0, 0); font-size: 13px; font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-style: normal;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"><br></font></div><div style="background-color: transparent; color: rgb(0, 0, 0); font-size: 13px; font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-style: normal;"><font face="Courier New, courier, monaco, monospace, sans-serif" size="2"><div style="font-size: 13px;">Certo de contar com a boa vontade de todos, agradeço antecipadamente.</div><div style="font-size: 13px;"><br></div><div style="font-size: 13px;">Att.,</div><div style="font-size: 13px;"><br></div><div style="font-size: 13px;">--</div><div style="font-size: 13px;">Wirton Macedo
Coutinho</div><div style="font-size: 13px;">Pesquisador Fitopatologia</div><div style="font-size: 13px;">Embrapa Algodão</div><div style="font-size: 13px;">Rua Oswaldo Cruz, 1143, Centenário</div><div style="font-size: 13px;">Campina Grande PB</div><div style="font-size: 13px;">CEP 28428-095</div></font></div></div><div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"> </div></div></body></html>