<font face="trebuchet ms,sans-serif">Alexandre,<br><br>O que é índice ambíguo e decisão na distribuição do índice?<br>No seu CMR, se você diz que é DIC, não precisa incluir efeito das repetições. E mesmo que fosse um caso de inclusão, elas deveriam ser codificadas como fator e não métricas. Se entendi, você questiona a apresentação dos resultados da glht(). Você tava esperando um conjunto de letras ao lado das estimativas pontuais, é isso? Isso não vai acontecer, mesmo porque, para atribuir letras sem ter problemas de interpretação, é necessário que todo par de contraste tenha o mesmo erro padrão, para que exista uma dms (diferença mínima significativa) constante. Isso só ocorre para o caso Gaussiano quando temos homocedasticidade e balanceamento. Veja na coluna Std.Err. que os erros padrões são variáveis. Isso porque na família glm a variância depende da média, logo, contraste entre médias tem variância que dependem das médias envolvidas. Na biblioteca HH existe um método gráfico para ajudar a interpretar essa saída <a href="http://www.oga-lab.net/RGM2/func.php?rd_id=HH:mmc">http://www.oga-lab.net/RGM2/func.php?rd_id=HH:mmc</a>.<br>
<br><span style="font-family: courier new,monospace;">#y1 <- c(rpois(100, lambda=5), rpois(100, lambda=10), rpois(100, lambda=5), rpois(100, lambda=7))##Geração da variável resposta</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">y1 <- c(mapply(rpois, lambda=c(5,10,5,7), MoreArgs=list(n=100)))</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#trat <- c(rep(1,100), rep(2,100), rep(3,100), rep(4,100))##Geração dos</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">trat <- gl(4,100)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#repeticao <- rep(1:100, 4) ##Imaginado um DIC colocando as repetições</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">#trat<-as.factor(trat)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># se é dic, não inclui efeito das repetições</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">m.1 <- glm(y1~trat,family="poisson")</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">anova.m.1 <- anova(m.1, test="Chi")</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">anova.m.1</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">summary(glht(m.1, linfct=mcp(trat="Tukey")))</span><br style="font-family: courier new,monospace;">
<br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Walmes Marques Zeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">fone: (+55) 41 3361 3573</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">VoIP: (3361 3600) 1053 1173</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">twitter: @walmeszeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">linux user number: 531218</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br>