[R-br] Inserir desvio padrão em gráfico de barras

Rodrigo Muniz uenf.rodrigo em gmail.com
Segunda Novembro 4 12:20:58 BRST 2013


Prezados membros da lista gostaria de pedir o auxílio de vocês.

A partir dos seguintes comando:

dados=cbind(c("Mata Nativa", "Citrus 37", "Citrus 37", "Citrus 24", "Citrus
18", "Citrus 31","Pasto 37"),
      tapply(arenosa[1:21, "ECA"], arenosa[1:21, "Trat"], mean),
      tapply(arenosa[1:21, "ECB"], arenosa[1:21, "Trat"], mean),
id=1:nrow(tapply(arenosa[1:21, "ECA"], arenosa[1:21, "Trat"], mean)))

colnames(dados)<-c("Classes", "0-10", "10-20", "id")
dados= as.data.frame(dados)
barchart(value~Classes, groups=variable, data=melt(dados, measure.var=2:3),
         auto.key=list(columns=2), horizon=FALSE, ylab="Estoque Carbono",
main="Área 3")

Eu gero a seguinte tabela de dados com o respectivo gráfico de barra.
       Classes             0-10    10-20    id
301  Mata Nativa     10.04     7.6        1
312   Citrus 37        12.36     23.12     2
322   Citrus 37        15.32     27.58     3
332   Citrus 24        4.38       8.85       4
342   Citrus 18        10.43     21.57     5
352   Citrus 31         9.19      18.03     6
364    Pasto 37        6.69      12.53     7

Eu gostaria de inserir neste gráfico de barras o desvio padrão de cada
barra.
Como eu poderia fazer isso?

Segue a tabela com os respectivos desvios padrão.
                      0-10         10-20
Mata Nativa 3.4250012     0.7715569
Citrus 37     2.0617064     4.1486182
Citrus 37     2.8501287     4.7112737
Citrus 24     0.1931321     0.3666515
Citrus 18     1.7760349     4.9343726
Citrus 31     0.6730775     1.6246333
Pasto 37      2.4274060    2.8266116

Desde já agradeço.

-- 
*Rodrigo A. Muniz*
Eng. Agrônomo. Ms Produção Vegetal (UENF)
Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP)
Cel (19) 8300-4333 (Pessoal)
Tel (19) 3375-1196 (Residêncial)
Tel (19) 3447-8537 (Trabalho)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131104/0a967500/attachment.html>


Mais detalhes sobre a lista de discussão R-br