library(StatDA)
data(chorizon)
pdf("fig-11-2.pdf",width=9,height=4.5)
par(mfrow=c(1,2),mar=c(4,4,2,2))
Be=chorizon[,"Be"]
Sr=chorizon[,"Sr"]
rBe=rank(Be)
rSr=rank(Sr)
plot(rBe, rSr, main=paste("Spearman correlation =",round(cor(rBe, rSr),2)),
xlab="Rank of Be in C-horizon",ylab="Rank of Sr in C-horizon",
cex.lab=1.2, pch=3, cex=0.7)
La=chorizon[,"La_INAA"]
Ce=chorizon[,"Ce_INAA"]
rLa=rank(La)
rCe=rank(Ce)
plot(rLa, rCe, main=paste("Spearman correlation =",round(cor(rLa, rCe),2)),
xlab="Rank of La in C-horizon",ylab="Rank of Ce in C-horizon",
cex.lab=1.2, pch=3, cex=0.7)
dev.off()