<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
</head>
<body text="#000000" bgcolor="#ffffff">
Walmes,<br>
<br>
Obrigado pela ajuda, instalei o pacote HH e vou estudar a
interpretação dos resultados;<br>
<br>
Não sei se esta correta minha interpretação mais chamo de índice
ambíguo quando o Tukey apresenta resultados do tipo ab, por
exemplo; <br>
<br>
Bom na verdade eu incluí as repetições pois futuramente pretendo
aplicar esta abordagem à dados reais advindos do campo que são em
DBC, onde algumas vezes encontro variação entre blocos e zeros
inflacionados, onde pretendo verificar se há significância ou não e
como estão influenciando na correta determinação do efeito dos
tratamentos, bom mesmo assim me esqueci de especificar com fator;<br>
<br>
Gostaria de te pedir se você conhece algum trabalho(s) de
referência(s) que utiliza desdobramentos de contrastes com GLM na
análise dos dados, pois apesar de ser adepto à abordagem da análise
de dados usando sua adequada distribuição de erros, tenho muita
dificuldade em justificar a metodologia de contrastes com GLM para
revisores de revistas na área a qual trabalho.<br>
<br>
Redobrados agradecimentos,<br>
<br>
Alexandre <br>
<br>
<br>
On 10/14/2011 11:12 AM, Walmes Zeviani wrote:
<blockquote
cite="mid:CAFU=EkY+MueCR6u1y+-TtKY2gqfqUd7dd29fyLdFZs=OvETgLA@mail.gmail.com"
type="cite"><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 moz-do-not-send="true"
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
moz-do-not-send="true" 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
moz-do-not-send="true" 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>
<pre wrap="">
<fieldset class="mimeAttachmentHeader"></fieldset>
_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>