<font face="trebuchet ms,sans-serif">Caro amigo sem nome,<br><br>O fato de ser glm não complica as coisas. O nível de referência segue ordem alfanumérica e não de entrada na planilha. A densidade está declarada como contínua e presente apenas com efeito de primeira ordem? Você vai ter que comparar as estimativas das médias para um valor fixo de densidade. Se usar o pacote contraste, isso é bem simples. Veja abaixo um exemplo simulado bem mediocre mas que é útil para entender o procedimento<br>
<br><span style="font-family: courier new,monospace;">tr <- gl(3,4) # níveis de tratamento</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">ds <- runif(length(tr), 1, 1.2) # densidade do solo</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">y <- rbinom(length(tr), prob=0.5, size=50) # germinação</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># tr é categórica, ds é númerica</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">m0 <- glm(cbind(y, 50-y)~tr*ds, family=binomial)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">anova(m0, test="Chisq")</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">summary(m0)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># compara os níveis de tratamento entre sí com ds=1.1</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">require(contrast)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># tr 1 contra tr 2 sendo a densidade do solo igual a 1.1</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">contrast(m0, list(tr="1", ds=1.1), list(tr="2", ds=1.1))</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># diferença em termos de preditor linear</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">predict(m0, data.frame(tr="1", ds=1.1))-</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">predict(m0, data.frame(tr="2", ds=1.1))</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># diferença em termos de resposta</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">predict(m0, data.frame(tr="1", ds=1.1), type="response")-</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">predict(m0, data.frame(tr="2", ds=1.1), type="response")</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># fazer as demais comparações. usar método Bonferroni para</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># manutenção do erro tipo 1.</span><br style="font-family: courier new,monospace;">
<br>Recomendo que você participe da nossa lista de discussão sobre R. Suas dúvidas serão bem recebidas e as respostas bem rápidas. Neste link estão informações para cadastro, envio e etiqueta na lista<br><br><span style="font-family: courier new,monospace;"><a href="http://www.leg.ufpr.br/doku.php/software:rbr">http://www.leg.ufpr.br/doku.php/software:rbr</a></span><br>
<br>À disposição.<br>Walmes.<br><br clear="all"></font><br clear="all"><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>

<br><br><div class="gmail_quote">2011/9/9  <span dir="ltr"><<a href="mailto:frdor85@hotmail.com">frdor85@hotmail.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Olá Walmes,<br>
<br>
É a primeira vez que pergunto algo aqui e gostaria muito de resolver um dúvida nas minhas análises. Faço mestrado trabalhando com emergência de plântulas herbáceas em diferentes tratamentos de serrapilheira(Adição, controle e exclusão de serrapilheira). O meu modelo inclui o teor de argila no solo e a interaçãoentre argila e os 3 tratamentos. No summary do teste não aparece os resultados referentes ao tratamento de adição de serrapilheira. Me falaram que o R compara o primeiro nível do tratamento com os outros (na minha planilha original a adição de serrapilheira é o primeiro nível dos tratamentos). Eu inverti a ordem de entrada dos tratamentos e a comparação foi a mesma!!! Como não é possível usar HDS Tukey, como posso obter todos os resultados da comparação entre os tratamentos de uma ancova, usando GLM???<br>

<br>
</blockquote></div><br>