7.3.2 The Concentration-Area Plot

# Fig. 7.3: Concentration are plot for Cu in O-Horizon library(StatDA) data(ohorizon) data(kola.background) data(bordersKola) Cu=ohorizon[,"Cu"] X=ohorizon[,"XCOO"] Y=ohorizon[,"YCOO"] # true representation of x and y axis of map for plot xwid=diff(range(X))/12e4 ywid=diff(range(Y))/12e4 pdf("fig-7-3.pdf",width=2*xwid,height=2*ywid) par(mfrow=c(2,2),mar=c(1.5,1.5,1.5,1.5)) concareaExampleKola(X,Y,Cu,log=T,zname="Cu in O-horizon [mg/kg]", x.logfinetick=c(2,5,10),y.logfinetick=c(10)) dev.off()

# Fig. 7.4: Concentration are plot for Cu in O-Horizon (upwards) library(StatDA) data(ohorizon) data(kola.background) Cu=ohorizon[,"Cu"] X=ohorizon[,"XCOO"] Y=ohorizon[,"YCOO"] # true representation of x and y axis of map for plot xwid=diff(range(X))/12e4 ywid=diff(range(Y))/12e4 pdf("fig-7-4.pdf",width=9,height=4.5) par(mfrow=c(1,2),mar=c(4,4,2,2)) concarea(X,Y,Cu,log=TRUE,zname="Cu in O-horizon [mg/kg]",borders="bordersKola", ifrev=FALSE, x.logfinetick=c(2,5,10),y.logfinetick=c(10)) dev.off()