[R-br] Adicionar os desvios padrão num gráfico

[Wenceslau]

Caros membros da lista,

Estou fazendo um gráfico simples (datas vs teores de ferro). Eu queria
adicionar os desvios padrão das avaliações dos teores em cada data. Que
foram calculados separadamente.
Mas não estou conseguindo. Talvez alguém tenha experiência com estes
gráficos e possa me ajudar.

Segue o script 


###########################
# SCRIPT PARA GRAFICOS DE ELEMENTOS NO SOLO 
#DATAS DE COLETAS

x<-1:30
x<-factor
x<- c("30/04/2010","13/05/2010","28/05/2010","15/06/2010","30/06/2010",
"15/07/2010", "30/07/2010" ,"15/08/2010",
"30/08/2010","15/09/2010","30/09/2010","15/10/2010","30/10/2010","15/11/2010
","30/11/2010","15/12/2010", "30/12/2010", "15/01/2011", "30/01/2011",
"15/02/2011", "28/02/2011"
,"15/03/2011","30/03/2011","15/04/2011","30/04/2011","15/05/2011","30/05/201
1","15/06/2011","30/06/2011",
      "15/07/2011")


#DADOS DE FE AVALIADOS NO LABORATORIO
y<-1:30
y <- c(0.278, 0.068, 0.064, 0.028, 0.030, 0.024, NA, 0.025, 0.039, 0.077,
0.081, NA, 0.043, NA, 0.343, NA, 0.195, 0.090, 0.027, 0.157, 0.122, 0.069,
0.197,
       0.107, 0.103, 0.154,0.164, 0.051,0.243,0.145)

# DESVIOS PADRAO DOS TEORES DE FE

SDFE10 <- c(0.294, 0.033,
0.058,0.013,0.006,0.006,NA,0.018,0.031,0,0,NA,0.021,NA,0.435,NA,0.057,0.049,
0.005,0.220,0.159,0.029,0.194,0.114,0.109,0.116,0.187,0,0.266,0.121)


#PLOTAR O GRAFICO 
plot(y, ylim=c(-0.1,0.7),  axes=F, ylab="Concentração de Ferro na solução do
solo (mg L-1)", xlab="Dia", pch=4, type="b")
#Adicionar os desvios padrão (SDFE10)
arrows(y,y-SDFE10,y, y+SDFE10, code = 3, angle = 90,length = 0.03)




Wenceslau Geraldes Teixeira

Embrapa Solos (www.cnps.embrapa.br)
Rua Jardim Botânico, 1024
CEP 22460-000
Rio de Janeiro - RJ - Brasil
Tel: 55 21 21794583

Skype: wgt007
E-mail: wenceslau em cnps.embrapa.br / wenceslau.teixeira em pq.cnpq.br

CV: http://lattes.cnpq.br/4136381865367906

________________________________

II Encontro Brasileiro de Física do Solo
2sd Brazilian Soil Physics Meeting
Rio de Janeiro - 06 a 10 de Maio de 2013
http://www.bspm2013.org/



-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121212/ca40e27a/attachment.html>