library(StatDA)
data(chorizon)
data(kola.background)
el=chorizon[,"As"]
X=chorizon[,"XCOO"]
Y=chorizon[,"YCOO"]
# true representation of x and y axis of map for plot
xwid=diff(range(X))/12e4
ywid=diff(range(Y))/12e4
pdf("fig-5-5.pdf",width=2*xwid,height=1*ywid)
par(mfrow=c(1,2),mar=c(1.5,1.5,1.5,1.5))
### NEW MAP: EDA symbols; take symbols according to boxplot classes
# generate plot with packground
plot(X,Y,frame.plot=FALSE,xaxt="n",yaxt="n",xlab="",ylab="",type="n")
plotbg(map.col=c("gray","gray","gray","gray"),add.plot=T)
SymbLegend(X,Y,el,type="boxplot",qutiles<-c(0,0.05,0.25,0.75,0.95,1),symbtype="EDA",symbmagn=0.8,
leg.position="topright",leg.title="As [mg/kg]",leg.title.cex=0.8,leg.round=2,leg.wid=4,
leg.just="right",cex.scale=0.75,xf=9e3,logscale=TRUE,accentuate=FALSE)
# Text
text(min(X)+diff(range(X))*3/5,max(Y),"Boxplot class selection",cex=0.7)
text(min(X)+diff(range(X))*3/5,max(Y)-diff(range(Y))/25,"EDA symbolset",cex=0.7)
# scalebar
scalebar(761309,7373050,861309,7363050,shifttext=-0.5,shiftkm=37e3,sizetext=0.8)
# North arrow
Northarrow(362602,7818750,362602,7878750,362602,7838750,Alength=0.15,Aangle=15,Alwd=1.3,Tcex=1.6)
### NEW MAP: EDA accentuated symbols; take symbols according to boxplot classes
# generate plot with packground
plot(X,Y,frame.plot=FALSE,xaxt="n",yaxt="n",xlab="",ylab="",type="n")
plotbg(map.col=c("gray","gray","gray","gray"),add.plot=T)
SymbLegend(X,Y,el,type="boxplot",qutiles<-c(0,0.05,0.25,0.75,0.95,1),symbtype="EDA",symbmagn=0.8,
leg.position="topright",leg.title="As [mg/kg]",leg.title.cex=0.8,leg.round=2,leg.wid=4,
leg.just="right",cex.scale=0.75,xf=9e3,logscale=TRUE,accentuate=TRUE)
# Text
text(min(X)+diff(range(X))*5/8,max(Y),"Boxplot class selection",cex=0.7)
text(min(X)+diff(range(X))*5/8,max(Y)-diff(range(Y))/25,"EDA symbolset",cex=0.7)
text(min(X)+diff(range(X))*5/8,max(Y)-2*diff(range(Y))/25,"(accentuated)",cex=0.7)
# scalebar
scalebar(761309,7373050,861309,7363050,shifttext=-0.5,shiftkm=37e3,sizetext=0.8)
# North arrow
Northarrow(362602,7818750,362602,7878750,362602,7838750,Alength=0.15,Aangle=15,Alwd=1.3,Tcex=1.6)
dev.off()