10.2.4 Box-Cox Transformation

# Fig. 10.4.: Fe_XRF in C-horizon, log10 and Box-Cox transformation library(StatDA) data(chorizon) Fe=chorizon[,"Fe_XRF"]/10000 pdf("fig-10-4.pdf",width=8,height=8) par(mfrow=c(3,1),mar=c(5,5,2,2)) plot(density(Fe),main="",xlab="Fe_XRF [wt.-%]",cex.lab=1.8, cex.axis=1.6) plot(density(log10(Fe)),main="",xlab="Fe_XRF [wt.-%]",cex.lab=1.8,xaxt="n", cex.axis=1.6) axis(1,at=log10(a<-sort(c((10^(-50:50))%*%t(c(2,5,10))))),labels=a, cex.axis=1.6) plot(density(Fe^0.31),main="",xlab="Fe_XRF [wt.-%]",cex.lab=1.8,xaxt="n", cex.axis=1.6) axis(1,at=(a<-sort(c((10^(-50:50))%*%t(c(2,5,10)))))^0.31,labels=a, cex.axis=1.6) dev.off()