<div>Obrigado mais uma vez pela informação.<br><br></div>
<div class="gmail_quote">Em 1 de junho de 2011 00:52, Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>></span> escreveu:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote"><font face="trebuchet ms,sans-serif">Eduardo,<br><br>O código abaixo mostra como não é difícil analisar os seus dados pela abordagem que discuti<br>
<br><span style="FONT-FAMILY: courier new, monospace">#------------------------------------------------------------------------------------------</span><br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace"># experimento com 11 tratamentos e 1 controle, resposta binomial n=25, 4 repetições</span><br style="FONT-FAMILY: courier new, monospace">
<br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace">trat <- c("controle", paste("trat", 1:11, sep=""))</span><br style="FONT-FAMILY: courier new, monospace">
<span style="FONT-FAMILY: courier new, monospace">da <- expand.grid(trat=trat, rept=1:4)</span><br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace">da$germinadas <- rbinom(nrow(da), size=25, prob=0.8)</span><br style="FONT-FAMILY: courier new, monospace">
<span style="FONT-FAMILY: courier new, monospace">levels(da$trat) # controle é o nível de referência</span><br style="FONT-FAMILY: courier new, monospace"><br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace"># ajuste um modelo linear generalizado</span><br style="FONT-FAMILY: courier new, monospace">
<span style="FONT-FAMILY: courier new, monospace">g0 <- glm(cbind(germinadas, 25-germinadas)~trat, data=da, family=binomial)</span><br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace">anova(g0, test="Chisq") # testa hipótese sobre os efeitos principais</span><br style="FONT-FAMILY: courier new, monospace">
<span style="FONT-FAMILY: courier new, monospace">summary(g0) # cada efeito é o contraste controle vs trat_i</span><br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace"># usar correção de bonferroni nas comparações com controle</span><br style="FONT-FAMILY: courier new, monospace">
<br style="FONT-FAMILY: courier new, monospace"><span style="FONT-FAMILY: courier new, monospace">#------------------------------------------------------------------------------------------</span><br style="FONT-FAMILY: courier new, monospace">
<br style="FONT-FAMILY: courier new, monospace">Se todo mundo continuar fazendo as coisas como todos já fizeram não vamos evoluir, certo?<br><br>À disposição.<br><font color="#888888">Walmes.<br><br clear="all"></font></font>
<div>
<div></div>
<div class="h5"><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/~walmes" 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></div></div><br>_______________________________________________<br>R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>José Eduardo Vargas Lopes de Aráujo<br>Engenheiro Agrônomo<br>Contato:(38)8811-6985                         <br><a href="mailto:duvargas@gmail.com">duvargas@gmail.com</a><br>
<br>