3.4.4 The Cumulative Probability Plot (CP-plot)

# Fig. 3.10.: QP plots for Au, log(Au) with line library(StatDA) data(chorizon) # need data Au=chorizon$Au n=length(Au) pdf("fig-3-10.pdf",width=9,height=4.5) par(mfcol=c(1,2),mar=c(4,4,2,2)) qpplot.das(Au,qdist=qlnorm,xlab="Au [ g/kg]", ylab="Probabilities of lognormal distribution", pch=3,cex=0.7) mtext(expression(mu),side=1,line=3.05,at=72.5,cex=1.2) qpplot.das(log10(Au),qdist=qnorm,xlab="Au [ g/kg]", ylab="Probabilities of standard normal distribution", pch=3,cex=0.7, logx=TRUE, logfinetick=c(5,10),logfinelab=c(10)) mtext(expression(mu),side=1,line=3.05,at=0.40,cex=1.2) dev.off()