library(StatDA)
data(chorizon)
# need Au data:
Au=chorizon$Au
n=length(Au)
pdf("fig-3-7.pdf",width=8,height=10)
par(mfcol=c(3,2),mar=c(4,4,2,2))
###################
edaplot(Au,H.freq=F,box=F,H.breaks=100,S.pch=3,S.cex=0.5,D.lwd=1.5,P.log=F,
P.main="",P.xlab="Au [ g/kg]",P.ylab="Density")
mtext(expression(mu),side=1,line=3.05,at=72.9,cex=0.9)
###################
edaplot(log10(Au),H.freq=F,box=F,H.breaks=15,S.pch=3,S.cex=0.5,D.lwd=1.5,
P.log=T,P.logfine=c(5,10),P.main="",P.xlab="Au [ g/kg]",
P.ylab="Density")
mtext(expression(mu),side=1,line=3.05,at=0.40,cex=0.9)
###################
edaplot(Au,H.freq=F,box=F,H.breaks=290,S.pch=3,S.cex=0.5,D.lwd=1.5,P.log=F,
P.main="",P.xlab="Au [ g/kg]",P.ylab="Density",P.xlim=c(0,6))
mtext(expression(mu),side=1,line=3.05,at=2.95,cex=0.9)
text(2.5,0.6,paste("Maximum =",round(max(Au),1)," g/kg"),pos=4,cex=1.3)
text(4.82,0.585,expression(mu),pos=4,cex=1.4)
#####
plot(sort(Au),((1:n)-0.5)/n,pch=3,cex=0.8,
main="",xlab="Au [ g/kg]",ylab="Probability",cex.lab=1.4)
abline(h=seq(0,1,by=0.1),lty=3,col=gray(0.5))
abline(v=seq(0,150,by=20),lty=3,col=gray(0.5))
mtext(expression(mu),side=1,line=3.05,at=72.9,cex=0.9)
#####
plot(sort(log10(Au)),((1:n)-0.5)/n,pch=3,cex=0.8,
main="",xlab="Au [ g/kg]",ylab="Probability",xaxt="n",cex.lab=1.4)
mtext(expression(mu),side=1,line=3.05,at=0.40,cex=0.9)
axis(1,at=log10(alog<-sort(c((10^(-50:50))%*%t(c(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))
#####
plot(sort(Au),((1:n)-0.5)/n,pch=3,cex=0.8,
main="",xlab="Au [ g/kg]",ylab="Probability",xlim=c(0,6),cex.lab=1.4)
abline(h=seq(0,1,by=0.1),lty=3,col=gray(0.5))
abline(v=seq(0,6),lty=3,col=gray(0.5))
mtext(expression(mu),side=1,line=3.05,at=2.95,cex=0.9)
text(2.5,0.35,paste("Maximum =",round(max(Au),1)," g/kg"),pos=4,cex=1.3)
text(4.82,0.335,expression(mu),pos=4,cex=1.4)
dev.off()