*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@embrapa.br wenceslau.teixeira@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@gmail.com skype: wgt007