<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>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>