PH    <- c(4.5, 4.3)

error <- c(0.5,0.2)

cidade <- c('Brasília','Itirapina')

da <- data.frame(cidade,PH,error)

library(ggplot2)

limits <- aes(ymax=PH + error, ymin=PH - error)

p <- ggplot(da, aes(cidade,PH))

p + geom_bar(fill='blue') + geom_errorbar(limits)      

Ultimamente estou ficando muito satisfeito com esta library. Para maiores detalhes,

Allaman

(S,f,P)