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@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 http://www.bspm2013.org/