library(StatDA)
data(chorizon)
As <- chorizon[,"As"]
n <- length(As)
pdf("fig-3-12.pdf",width=8,height=10)
par(mfcol=c(3,2),mar=c(4,4,2,2))
##### ECDF original
plot(sort(As),((1:n)-0.5)/n,pch=3,cex=0.8,
main="",xlab="As in C-horizon [mg/kg]",ylab="Probability",cex.lab=1.2)
abline(h=seq(0,1,by=0.1),lty=3,col=gray(0.5))
abline(v=seq(0,35,by=5),lty=3,col=gray(0.5))
##### QP original
qpplot.das(As,qdist=qnorm,xlab="As in C-horizon [mg/kg]", cex.lab=1.2,
ylab="Probabilities of standard normal distribution", pch=3,cex=0.7, logx=FALSE)
##### PP original
ppplot.das(As,pdist=pnorm,xlab="Probability of As in C-horizon", cex.lab=1.2,
ylab="Probabilities of standard normal distribution", pch=3,cex=0.7)
##### ECDF log
plot(sort(log10(As)),((1:n)-0.5)/n,pch=3,cex=0.8,
main="",xlab="As in C-horizon [mg/kg]",ylab="Probability",xaxt="n",cex.lab=1.2)
axis(1,at=log10(alog<-sort(c((10^(-50:50))%*%t(c(2,5,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))
##### QP log
qpplot.das(log10(As),qdist=qnorm,xlab="As in C-horizon [mg/kg]", cex.lab=1.2,
ylab="Probabilities of standard normal distribution", pch=3,cex=0.7, logx=TRUE,
logfinetick=c(2,5,10),logfinelab=c(2,5,10))
##### PP log
ppplot.das(log10(As),pdist=pnorm,xlab="Probability of log10(As)",cex.lab=1.2,
ylab="Probabilities of standard normal distribution", pch=3,cex=0.7)
dev.off()