Rodrigo, Você pode usar coisas prontas no R, como http://addictedtor.free.fr/graphiques/graphcode.php?graph=54 http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=72 Ou pode construir com os recursos gráficos do R, de onde é possível ter mais controle, veja da <- data.frame(trat=gl(5,8)) da$y <- as.numeric(da$trat)+rnorm(nrow(da)) m0 <- lm(y~trat, da) new <- data.frame(trat=levels(da$trat)) new$pred <- predict(m0, newdata=new, interval="confidence") str(new) ylim <- c(0, max(new$pred)*1.05) bp <- barplot(new$pred[,1], ylim=ylim) arrows(bp, new$pred[,2], bp, new$pred[,3], code=3, angle=90) box() À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================