library(StatDA)
data(bhorizon)
data(chorizon)
data(moss)
data(ohorizon)
# only samples from same locations in ALL layers:
c.ind=NULL
b.ind=NULL
m.ind=NULL
h.ind=NULL
c.id=chorizon[,1]
b.id=bhorizon[,1]
m.id=moss[,1]
h.id=ohorizon[,1]
for (i in 1:1000){
if (sum(c.id==i)+sum(b.id==i)+sum(m.id==i)+sum(h.id==i)==4){
c.ind=c(c.ind,which(c.id==i))
b.ind=c(b.ind,which(b.id==i))
m.ind=c(m.ind,which(m.id==i))
h.ind=c(h.ind,which(h.id==i))
}
}
# Cu data
c.Cu=chorizon[c.ind,"Cu"]
b.Cu=bhorizon[b.ind,"Cu"]
pdf("fig-9-5.pdf",width=9,height=4.5)
par(mfrow=c(1,2),mar=c(4,4,2,2))
boxplot(log10(c.Cu)-log10(b.Cu),notch=T,horizontal=T,pch=3,cex=0.8,
xlab="Difference of Cu in B- and C-horizon",cex.lab=1.2)
abline(v=0,lty=2)
boxplot(log10(c.Cu),log10(b.Cu),horizontal=T,xaxt="n",pch=3,cex=0.8,notch=T,
xlab="Cu [mg/kg]",names=c("C-horizon","B-horizon"),cex.lab=1.2)
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
abline(a=0.5,b=1,lwd=2,lty=2,col=grey(0.5))
abline(a=2.5,b=-1,lwd=2,lty=2,col=grey(0.5))
dev.off()