6.7 Scatterplot Matrix

# Fig. 6.13.: draftman's plot library(StatDA) data(chorizon) pdf("fig-6-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 <- x[!is.na(litolog),] litolog <- litolog[!is.na(litolog)] litolog <- factor(litolog, labels=c("AB","PG","AR","LPS")) scatterplot.matrix(x.sub, diagonal="none", smooth=F, reg.line=F, robust=F, ellipse=F, transform=F, groups=F, by.groups=F, pch=3, col=rep(1,length(litolog)), lwd=1, cex=0.3, cex.labels=1.9, legend.plot=F,xaxt="n",yaxt="n") dev.off()

# Fig. 6.14.: draftman's plot library(StatDA) data(chorizon) pdf("fig-6-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("AB","PG","AR","LPS")) scatterplot.matrix(x.sub, diagonal="none", smooth=F, reg.line=F, robust=F, ellipse=F, transform=F, groups=F, by.groups=F, pch=3, col=rep(1,length(litolog)), lwd=1, cex=0.3, cex.labels=1.9, legend.plot=F,xaxt="n",yaxt="n") dev.off()