library(StatDA)
data(chorizon)
pdf("fig-11-3.pdf",width=9,height=4.5)
par(mfrow=c(1,2),mar=c(4,4,2,2))
Be=chorizon[,"Be"]
Sr=chorizon[,"Sr"]
RobCor.plot(log10(Be),log10(Sr),
xlab="Be in C-horizon [mg/kg]",ylab="Sr in C-horizon [mg/kg]",
cex.lab=1.2, pch=3, cex=0.7, xaxt="n", yaxt="n",
colC=1,colR=1,ltyC=2,ltyR=1)
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)
legend("topleft",legend=c("Pearson","Robust"),lty=c(2,1))
La=chorizon[,"La_INAA"]
Ce=chorizon[,"Ce_INAA"]
RobCor.plot(log10(La),log10(Ce),
xlab="La in C-horizon [mg/kg]",ylab="Ce in C-horizon [mg/kg]",
cex.lab=1.2, pch=3, cex=0.7, xaxt="n", yaxt="n",
colC=1,colR=1,ltyC=2,ltyR=1)
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)
legend("topleft",legend=c("Pearson","Robust"),lty=c(2,1))
dev.off()