library(StatDA)
data(chorizon)
Cr=chorizon[,"Cr"]
Fe=chorizon[,"Fe_XRF"]/10000
pdf("fig-10-6.pdf",width=8,height=8)
par(mfrow=c(3,1),mar=c(5,5,2,2))
densCr <- density(log10(Cr))
densFe <- density(log10(Fe))
plot(densCr,main="",xlab="log10-transformation",cex.lab=1.8, cex.axis=1.6,
xlim=c(min(densCr$x,densFe$x),max(densCr$x,densFe$x)),
ylim=c(0,max(densCr$y,densFe$y)))
lines(densFe,lty=2)
text(0.15,1.5,"Fe_XRF",cex=1.8)
text(1.7,1.2,"Cr",cex=1.8)
densCr <- density(scale(log10(Cr),T,F))
densFe <- density(scale(log10(Fe),T,F))
plot(densCr,main="",xlab="log10-transformation and centring",cex.lab=1.8, cex.axis=1.6,
xlim=c(min(densCr$x,densFe$x),max(densCr$x,densFe$x)),
ylim=c(0,max(densCr$y,densFe$y)))
lines(densFe,lty=2)
text(-0.3,1.9,"Fe_XRF",cex=1.8)
text(0.5,0.6,"Cr",cex=1.8)
densCr <- density(scale(log10(Cr),T,T))
densFe <- density(scale(log10(Fe),T,T))
plot(densCr,main="",xlab="log10-transformation, centring and scaling",
cex.lab=1.8, cex.axis=1.6,
xlim=c(min(densCr$x,densFe$x),max(densCr$x,densFe$x)),
ylim=c(0,max(densCr$y,densFe$y)))
lines(densFe,lty=2)
text(1.8,0.25,"Fe_XRF",cex=1.8)
text(-0.6,0.42,"Cr",cex=1.8)
dev.off()