O código que usei foi este:

require(sciplot)
bargraph.CI(pH, SHAM, group = time, data = sham, 
xlab = "Tratamento", ylab = "Sham",
cex.lab = 1.2, x.leg = 1, 
col = c("gray18", "gray37", "gray60"),
density = c(0,0.8),  legend.text = NULL)

legend("topright", c("2","4","8"), 
       pch = c(15,15,15),
col=c("gray18","gray37","gray60"), 
title="Time (h)")

sendo igual ao que você enviou Éder.
Porém eu ainda gostaria de acrescentar os niveis de mais um fator no eixo x.


2014-01-27 Éder Comunello <comunello.eder@gmail.com>
Um outro exemplo, dessa vez usando sciplot::bargraph.CI...

### <code r>
dURL <- 'http://www.ualberta.ca/~ahamann/teaching/renr480/Lab5d.zip'
download.file(dURL, basename(dURL), mode='wb')
unzip(basename(dURL),'lentils.csv')
dat <- read.table('lentils.csv', head=T, sep=',')
attach(dat)

require(sciplot)
bargraph.CI(x.factor=FARM, group=VARIETY, response=YIELD, legend=T, 
 xlab="Location", ylab="Yield (kg/ha)", ylim=c(0,800), 
 cex.names=1, cex.lab=1, 
 x.leg=6, y.leg=700, cex.leg=1, leg.lab=c("Var A","Var B","Var C"), 
 col=grey.colors(3), 
 uc=T, lc=T, err.width=0.1, err.col="black", err.lty=1) 
### </code>

Fonte/Source: http://www.ualberta.ca/~ahamann/teaching/renr480/Lab5D.pdf

Éder Comunello <comunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]





_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.



--

Natália Martins
Contato: (35) 91812482