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ê"))