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",
header=TRUE,sep=";")
validation<-read.csv("jura_validation.CSV",
header=TRUE,sep=";")
x=prediction$Cr; y=validation$Cr
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_17.jpeg",width=5,height=5,
#units="in",res=300)
qqplot(x,y,cex=1.00,pch=21,bg="red",
xlab="Cr-Prediction",
ylab="Cr-Validation",col="red")
abline(0,1,lwd=3,col="blue")
#dev.off()