A explicação do script a seguir está no artigo R na Prática – Parte 27.
setwd("C:\\geoStats\\dados\\jura_goovaerts")
prediction=read.csv("jura_prediction.CSV",sep=";",
header=TRUE)
xp=prediction$Xloc; yp=prediction$Yloc
validation=read.csv("jura_validation.CSV",sep=";",
header=TRUE)
xv=validation$Xloc; yv=validation$Yloc
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_16.jpeg",width=5,height=5,
#units="in",res=300)
plot(xp,yp,pch=21,cex=1.25,bg="red",xlim=c(0,5),
ylim=c(0,6),xlab="LESTE (km)",ylab="NORTE (km)")
points(xv,yv,pch=23,bg="blue",cex=1.25)
points(0,6,pch=21,bg="red",cex=1.25)
points(0,5.75,pch=23,bg="blue",cex=1.25)
text(0.15,6,pos=4,labels="Prediction data")
text(0.15,5.75,pos=4,labels="Validation data")
#dev.off()