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: [email protected] / [email protected]
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/
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
