library(StatDA)
data(ohorizon)
el=log10(ohorizon[,"Cu"])
pdf("fig-7-1.pdf",width=9,height=4.5)
par(mfcol=c(1,2),mar=c(4,4,2,2))
n=length(el)
plot(sort(el),((1:n)-0.5)/n,pch=3,cex=0.8,xlim=c(0.1,max(el)),
main="",xlab="Cu in O-horizon [mg/kg]",ylab="Probability",cex.lab=1,xaxt="n",cex.lab=1.2)
axis(1,at=log10(alog<-sort(c((10^(-50:50))%*%t(c(10))))),labels=alog)
#abline(h=seq(0,1,by=0.1),lty=3,col=gray(0.5))
#abline(v=log10(alog),lty=3,col=gray(0.5))
abline(v=c(mean(el)-2*sd(el),mean(el)+2*sd(el)),lty=2)
text(mean(el)-2*sd(el),0.5,"MEAN-2*SD",srt=90,adj=c(0.5,-0.3),cex=0.9)
text(mean(el)+2*sd(el),0.3,"MEAN+2*SD",srt=90,adj=c(0.5,-0.3),cex=0.9)
abline(v=c(median(el)-2*mad(el),median(el)+2*mad(el)),lty=4)
text(median(el)-2*mad(el),0.5,"MEDIAN-2*MAD",srt=90,adj=c(0.5,-0.3),cex=0.9)
text(median(el)+2*mad(el),0.3,"MEDIAN+2*MAD",srt=90,adj=c(0.5,-0.3),cex=0.9)
qpplot.das(el,qdist=qnorm,xlab="Cu in O-horizon [mg/kg]",
ylab="Probability [%]", pch=3,cex=0.7, logx=TRUE,
logfinetick=c(0),logfinelab=c(10),cex.lab=1.2,gridy=FALSE,
xlim=c(0.1,max(el)))
#axis(1,at=log10(alog<-sort(c((10^(-50:50))%*%t(c(5,10))))),labels=alog)
#abline(v=log10(alog),lty=3,col=gray(0.5))
abline(v=c(mean(el)-2*sd(el),mean(el)+2*sd(el)),lty=2)
text(mean(el)-2*sd(el),0.5,"MEAN-2*SD",srt=90,adj=c(0.5,-0.3),cex=0.9)
text(mean(el)+2*sd(el),-1.2,"MEAN+2*SD",srt=90,adj=c(0.5,-0.3),cex=0.9)
abline(v=c(median(el)-2*mad(el),median(el)+2*mad(el)),lty=4)
text(median(el)-2*mad(el),0.5,"MEDIAN-2*MAD",srt=90,adj=c(0.5,-0.3),cex=0.9)
text(median(el)+2*mad(el),-1.2,"MEDIAN+2*MAD",srt=90,adj=c(0.5,-0.3),cex=0.9)