10.2.2 Power Transformation

# Fig. 10.2.: Density traces of Ca in C-horizon with power transformation
library(StatDA)
data(chorizon)
Cr=chorizon[,"Cr"]

pdf("fig-10-2.pdf",width=8,height=8)
par(mfrow=c(3,1),mar=c(5,5,2,2))

plot(density(Cr),main="",xlab="Cr [mg/kg]",cex.lab=1.8,
cex.axis=1.6)

plot(density(Cr^(1/4)),main="",xlab="Cr power transformed (1/4)",cex.lab=1.8,
cex.axis=1.6)

plot(density(Cr^(-1/4)),main="",xlab="Cr power transformed (-1/4)",cex.lab=1.8,
cex.axis=1.6)

dev.off()