Uma tentativa...
### <code r>
dat <- structure(list(time = c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 4L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
8L), pH = c(5.8, 5.8, 5.8, 4.5, 4.5, 4.5, 4, 4, 4, 5, 5, 5, 5.8,
5.8, 5.8, 4.5, 4.5, 5.8, 4.5, 4.5, 4.5, 4, 4, 4, 5, 5, 5), sinal = c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), SHAM = c(0.002, 0.003,
0.008, 0.004, 0.002, 0.006, 0.043, 0.024, 0.039, 0.013, 0.008,
0.004, 0.007, 0.001, 0.007, 0.017, 0.005, 0.019, 0.184, 0.188,
0.198, 0.119, 0.134, 0.231, 0.072, 0.135, 0.123)), .Names = c("time",
"pH", "sinal", "SHAM"), class = "data.frame", row.names = c(NA,
-27L))
dat$new <- with(dat, paste0(format(dat$pH, dig=2), " (", sinal, ")")); dat
require(sciplot)
bargraph.CI(x.factor=new, response=SHAM, group = time, data = dat,
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)")