[R-br] Barras de erro em barplot

Possuo uma planilha distribuída da seguinte forma: cidade       pH     desvpad..... Brasília       4,5     0,5 Itirapina       4,3     0,2 Como faço p plotar acrescentar o desvio padrão superior e inferios a minha figura partindo de barplot(pH~cidade,data=dados, .....) Muito obrigado.