library(StatDA)
data(chorizon)
Cu=chorizon[,"Cu"]
pdf("fig-9-3.pdf",width=9,height=4.5)
par(mfrow=c(1,2),mar=c(4,4,2,2))
plot(d <- density(log10(Cu)),main="",xlab="Cu in C-horizon [mg/kg]",xaxt="n",cex.lab=1.2)
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
abline(v=median(log10(Cu)))
lines(d$x,dnorm(d$x,log10(18),sd(log10(Cu))),lty=2)
abline(v=log10(18),lty=2)
legend(log10(28),1.32,legend=c("Empirical","Hypothetical"),lty=c(1,2),cex=0.85)
coun=factor(chorizon[,"COUN"],levels=c("FIN","NOR","RUS"),labels=c("Finland","Norway","Russia"))
dFIN=density(log10(Cu[coun=="Finland"]))
dNOR=density(log10(Cu[coun=="Norway"]))
dRUS=density(log10(Cu[coun=="Russia"]))
xl=c(min(dFIN$x,dNOR$x,dRUS$x),max(dFIN$x,dNOR$x,dRUS$x))
yl=c(0,max(dFIN$y,dNOR$y,dRUS$y))
plot(dFIN,main="",xlab="Cu in C-horizon [mg/kg]",xaxt="n",
xlim=xl,ylim=yl,cex.lab=1.2)
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
lines(dNOR,lty=4)
lines(dRUS,lty=5)
abline(v=log10(18),lty=2)
legend(log10(47),1.52,legend=c("Finland","Norway","Russia"),lty=c(1,4,5),cex=0.85)
dev.off()