Adicionar os desvios padrão num gráfico
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/
Alguns cuidados você deve tomar com datas. Representá-las com string não vai reconhecer a ordem cronológica das medidas. É necessário converter para data. Segue um CMR sobre os dados fornecidos. Outras opções existes no pacote gplots, lattice, ggplot2. 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") 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) 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) # colocar num data.frame, eliminar NA, converter para data da <- data.frame(x=as.POSIXlt(x, format="%d/%m/%Y"), y=y, sd=SDFE10) da <- na.omit(da) da$x plot(y~x, data=da, ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})), xlab="Dia", pch=4, type="b", xaxt="n", ylim=c(-0.1,0.7)) with(da, arrows(x, y-sd, x, y+sd, code=3, angle=90,length=0.03)) #axis.POSIXct(side=1, at=da$x, format="%d/%m") axis.POSIXct(side=1, at=seq(min(da$x), max(da$x), by="month"), format="%d/%m") require(gplots) help(plotCI, help_type="html") À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Walmes Zeviani -
Wenceslau