<font face="trebuchet ms,sans-serif">Rodrigo,<br><br>Você pode usar coisas prontas no R, como<br><br><a href="http://addictedtor.free.fr/graphiques/graphcode.php?graph=54">http://addictedtor.free.fr/graphiques/graphcode.php?graph=54</a><br>
<a href="http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=72">http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=72</a><br><br>Ou pode construir com os recursos gráficos do R, de onde é possível ter mais controle, veja<br>
<br><span style="font-family: courier new,monospace;">da <- data.frame(trat=gl(5,8))</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">da$y <- as.numeric(da$trat)+rnorm(nrow(da))</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">m0 <- lm(y~trat, da)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">new <- data.frame(trat=levels(da$trat))</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">new$pred <- predict(m0, newdata=new, interval="confidence")</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">str(new)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">ylim <- c(0, max(new$pred)*1.05)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">bp <- barplot(new$pred[,1], ylim=ylim)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">arrows(bp, new$pred[,2], bp, new$pred[,3], code=3, angle=90)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">box()</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>