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/2011","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@cnps.embrapa.br / wenceslau.teixeira@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