[R-br] Adicionar um intervalo (quebra - axis.break) no eixo y em gráficos no R

[Wenceslau Teixeira]

*Caros colegas da lista R,*

* *

*Envio este post para saber se alguém pode me auxiliar ou dar algumas dicas
de bibliografia sobre como colocar uma quebra (intervalo) no eixo y em
gráficos.*

* *

*Segue o código reproduzível abaixo cujo intervalo ultrapassa o eixo y e
não gostaria de aumentar os limites de Maximo no eixo y  pois a maioria dos
dados acaba ficando restrita num pequeno intervalo que dificulta a
visualização. Eu encontrei o comando (**axis.break) mas não consegui
implementar no grafico.
*

* *

*Sds estatísticas e antecipadamente agradeço alguma orientação.*

* *

*Wenceslau Teixeira *

*Embrapa Solos*

* *

*#X SÃO AS DATAS*

* *

*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*

*y <-
c(0.067,0.094,0.039,0.048,0.123,0.080,0.036,0.179,0.069,0.601,NA,NA,0.057,NA,NA,0.028,0.273,0.052,0.045,0.476,0.154,0.314,0.120,0.100,0.079,1.208,8.388,3.307,0.774,0.244)
*

* *

*#SDXX E O DESVIO PADRAO DO ELEMENTO FE A 10 CM QUE SERA FEITO O GRAFICO*

* *

*SDFE10
<-c(0.018,0.028,0.013,0.047,0.132,0.056,0.011,0.177,0.008,0.000,NA,NA,0.028,NA,NA,0.000,0.387,0.030,0.000,0.622,0.037,0.157,0.054,0.030,0.031,1.907,9.203,4.629,1.041,0.189)
*

* *

*# colocar num data.frame, eliminar NA, converter para data*

*da <- data.frame(x=as.POSIXlt(x, format="%d/%m/%Y"), y=y,cex.axis=0.5,
sd=SDFE10)*

*#da <- na.omit(da)*

*da$x*

* *

*#fazer o grafico *

* *

*plot(y~x,tck=0.01,
data=da,ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),tck=0.01,xlab="",
pch=24,cex.axis=0.5, type="p", xaxt="n",ylim=c(-0.2,0.8))*

*with(da, arrows(x, y-sd, x, y+sd,cex.axis=0.5, code=3,
angle=90,length=0.03))*

* *

*axis.POSIXct(side=1, at=da$x, las=2, format="%d/%m/%Y",cex.axis=0.5)*

* *

*#adicionado os outros dados da prof50 cm *

* *

*par(new=TRUE)*

* *

*#Y SÃO OS DADOS DO FE50CM*

* *

*y1 <-
c(0.097,0.029,0.026,0.053,0.056,0.014,0.015,0.122,NA,0.140,0.160,NA,0.088,NA,NA,NA,NA,NA,0.017,0.070,0.052,0.173,0.079,0.037,0.050,0.420,0.280,0.030,0.057,0.037)
*

* *

*#SDXX E O DESVIO PADRAO FERRO A 50 CM  FEITO O GRAFICO *

* *

*SDFE50 <-
c(0.078,0.000,0.012,0.057,0.011,0.000,0.000,0.000,NA,0.000,0.000,,0.000,NA,NA,NA,NA,NA,0.000,0.027,0.029,0.062,0.019,0.006,0.000,0.114,0.227,0.000,0.000,0.000)
*

*# break axis*

*#axis.break(2, 22, style="zigzag")*

* *

* *

*#fazer o segundo grafico *

* *

*plot(y1~x, data=da,cex.axis=0.5,*

*     ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),xlab="",tck=0.01,
pch=19, type="p", xaxt="n",ylim=c(-0.2,0.8))*

*with(da, arrows(x, y1-sd, x, y1+sd, code=3, angle=90,length = 0.03))*

* *

*#Adicionando uma legenda*

* *

*legend("topright", c("10 cm", "50 cm"),  pch = c(24,19))*

* *

*#informação sobre O tratamento *

*legend("topleft", c("Tratamento - TENS Dendê"))*

* *


-- 
Wenceslau Geraldes Teixeira
Pesquisador Embrapa


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.teixeira em embrapa.br
        wenceslau.teixeira em pq.cnpq.br

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

II ENCONTRO BRASILEIRO DE FÍSICA DO SOLO - 2013
2nd Brazilian Soil Physics Meeting
Rio de Janeiro, 06 a 10 de Maio de 2013

http://www.bspm2013.org/





-- 
Wenceslau Teixeira

wgt007 em gmail.com
skype: wgt007
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130111/0f7470ed/attachment.html>