library(StatDA)
data(chorizon)
Cu=chorizon[,"Cu"]
pdf("fig-9-4.pdf",width=6,height=4.5)
par(mfrow=c(1,1),mar=c(4,4,2,2))
c2=chorizon[,"COUN"]
c2[c2=="NOR"] <- "FIN"
coun2=factor(c2,levels=c("FIN","RUS"),labels=c("Finland and Norway", "Russia"))
boxplot(log10(Cu)~coun2,horizontal=T,xaxt="n",pch=3,cex=0.8,notch=T,
xlab="Cu in C-horizon [mg/kg]",cex.lab=1.2)
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
dev.off()