8.8 Data Subsets in the Scatterplot Matrix

# Fig. 8.13.: draftman's plot with subgroups library(StatDA) data(chorizon) pdf("fig-8-13.pdf",width=9,height=9) par(mfrow=c(1,1),mar=c(0.1,0.1,0.1,0.1)) x=chorizon[,c("Ca","Cu","Mg","Na","P","Sr","Zn")] lit=chorizon[,"LITO"] litolog=rep(NA, length(lit)) litolog[lit==10] <- 1 litolog[lit==52] <- 2 litolog[lit==81 | lit==82 | lit==83] <- 3 litolog[lit==7] <- 4 x.sub <- log10(x[!is.na(litolog),]) litolog <- litolog[!is.na(litolog)] litolog <- factor(litolog, labels=c(1:4)) scatterplot.matrix(~Ca+Cu+Mg+Na+P+Sr+Zn | litolog, reg.line="none", smooth=F, span=0.5, diagonal= "none", by.groups=TRUE, data=x.sub, pch=c(3,4,1,22), col=rep(gray(0.1),length(litolog)), xaxt="n",yaxt="n",legend.plot=F,cex.labels=1.9,cex=1.2) mtext("Caledonian Sediments", side=1, line=-2,outer=F, at=0.03,adj=0) points(0.01,-0.025,pch=3,cex=1.3,col=1) mtext("Basalts", side=1, line=-2,outer=F, at=0.34,adj=0) points(0.32,-0.025,pch=4,cex=1.3,col=1) mtext("Alkaline Rocks", side=1, line=-2,outer=F, at=0.50,adj=0) points(0.48,-0.025,pch=1,cex=1.3,col=1) mtext("Granites", side=1, line=-2,outer=F, at=0.74,adj=0) points(0.72,-0.025,pch=22,cex=1.3,col=1) dev.off()

# Fig. 8.14.: draftman's plot with subgroups and lines library(StatDA) data(chorizon) pdf("fig-8-14.pdf",width=9,height=9) par(mfrow=c(1,1),mar=c(0.1,0.1,0.1,0.1)) x=chorizon[,c("Ca","Cu","Mg","Na","P","Sr","Zn")] lit=chorizon[,"LITO"] litolog=rep(NA, length(lit)) litolog[lit==10] <- 1 litolog[lit==52] <- 2 litolog[lit==81 | lit==82 | lit==83] <- 3 litolog[lit==7] <- 4 x.sub <- log10(x[!is.na(litolog),]) litolog <- litolog[!is.na(litolog)] litolog <- factor(litolog, labels=c(1:4)) scatterplot.matrix(~Ca+Cu+Mg+Na+P+Sr+Zn | litolog, reg.line=lm, smooth=F, span=0.5, diagonal= "none", by.groups=TRUE, data=x.sub, pch=c(3,4,1,22), col=rep(gray(0.1),length(litolog)), xaxt="n",yaxt="n",legend.plot=F,cex.labels=1.9,cex=1.2) mtext("Caledonian Sediments", side=1, line=-2,outer=F, at=0.03,adj=0) points(0.01,-0.025,pch=3,cex=1.3,col=1) mtext("Basalts", side=1, line=-2,outer=F, at=0.34,adj=0) points(0.32,-0.025,pch=4,cex=1.3,col=1) mtext("Alkaline Rocks", side=1, line=-2,outer=F, at=0.50,adj=0) points(0.48,-0.025,pch=1,cex=1.3,col=1) mtext("Granites", side=1, line=-2,outer=F, at=0.74,adj=0) points(0.72,-0.025,pch=22,cex=1.3,col=1) dev.off()