library(StatDA)
data(chorizon)
pdf("fig-9-2.pdf",width=9,height=4.5)
par(mfrow=c(1,2),mar=c(4,4,2,2))
Co=log10(chorizon[,"Co"])
lit=chorizon[,"LITO"]
lit1=rep(NA,length(lit))
lit1[lit==9] <- 1
lit1[lit==20] <- 2
lit1[lit==32] <- 3
lit1[lit==51] <- 4
lit1[lit==52] <- 5
G1=density(Co[lit==9])
G2=density(Co[lit==20])
G3=density(Co[lit==32])
G4=density(Co[lit==51])
G5=density(Co[lit==52])
Gall=density(Co[!is.na(lit1)],width=.4)
plwd=c(1,1,1.6,1,1.5)
plot(Gall, main="",xlab="Co in C-horizon [mg/kg]",ylab="Density",
xaxt="n",cex.lab=1.2)
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
plot(0,0,xlim=c(min(G1$x,G2$x,G3$x,G4$x,G5$x),max(G1$x,G2$x,G3$x,G4$x,G5$x)),
ylim=c(min(G1$y,G2$y,G3$y,G4$y,G5$y),max(G1$y,G2$y,G3$y,G4$y,G5$y)),
main="",xlab="Co in C-horizon [mg/kg]",ylab="Density",
xaxt="n",cex.lab=1.2)
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
lines(G1,lty=1,lwd=plwd[1])
lines(G2,lty=2,lwd=plwd[2])
lines(G3,lty=3,lwd=plwd[3])
lines(G4,lty=4,lwd=plwd[4])
lines(G5,lty=5,lwd=plwd[5])
legend("topleft",legend=c(9,20,32,51,52),lty=c(1,2,3,4,5),lwd=plwd,col=1)
dev.off()