library(StatDA)
data(chorizon)
Be=chorizon[,"Be"]
Sr=chorizon[,"Sr"]
pdf("fig-11-1.pdf",width=9,height=4.5)
par(mfrow=c(1,2),mar=c(4,4,2,2))
plot(Be, Sr, main=paste("Pearson correlation =",round(cor(Be, Sr),2)),
xlab="Be in C-horizon [mg/kg]",ylab="Sr in C-horizon [mg/kg]",
cex.lab=1.2, pch=3, cex=0.7)
plot(log10(Be), log10(Sr), main=paste("Pearson correlation =",
round(cor(log10(Be), log10(Sr)),2)),
xlab="Be in C-horizon [mg/kg]",ylab="Sr in C-horizon [mg/kg]",
cex.lab=1.2, pch=3, cex=0.7, xaxt="n", yaxt="n")
axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
axis(2,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a)
dev.off()