10.2.3 Log(arithmic) Transformation

# Fig. 10.3.: Cu against Cr in C-horizon, also log-transformed library(StatDA) data(chorizon) Cu=chorizon[,"Cu"] Cr=chorizon[,"Cr"] pdf("fig-10-3.pdf",width=8,height=4) par(mfrow=c(1,2),mar=c(4,4,2,2)) plot(Cu,Cr,main="",xlab="Cu in C-horizon [mg/kg]",ylab="Cr in C-horizon [mg/kg]",cex.lab=1.2,pch=3,cex=0.6) plot(log10(Cu),log10(Cr),main="",xlab="Cu in C-horizon [mg/kg]",ylab="Cr in C-horizon [mg/kg]", cex.lab=1.2,xaxt="n",yaxt="n",pch=3,cex=0.6) 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()