11.7 Correlation Analysis and Inhomogeneous Data

# Fig.11.5.: Plot correlation matrix with subsets library(StatDA) data(chorizon) 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")) pdf("fig-11-5.pdf",width=5,height=5) par(mfrow=c(1,1),mar=c(0.1,0.1,0.1,0.1)) CorGroups(log10(x), grouping=litolog, labels1=dimnames(x)[[2]],labels2=dimnames(x)[[2]], legend=c("Caledonian Sediments","Basalts","Alkaline Rocks","Granites"),ndigits=2) dev.off()