
415Geostatistics
12.6 COMPUTER SESSION: GEOSTATISTICS, KRIGING
12.6.1 orDinary kriGinG
Use the lab12/example-1x1.txt le, which is in geoEAS format. It contains x, y coordinates and
marks labeled “z”
xyz marked point pattern
3
x
y
z
0.4785 0.4668 103
0.7734 0.5599 110
0.934 0.8124 113
0.3497 0.5704 70
0.9335 0.4743 97
0.4852 0.9415 114
0.2065 0.3103 122
… etc …
First, read the le, convert to point pattern, and plot the point pattern marking by quartiles using
the script
xy <- scan.geoeas.ppp(“lab12/example-1x1.txt”)
xy.ppp <- ppp(xy$x, xy$y, marks=xy$z)
xyz <- point(xy)
plot.point.bw(xyz,v=‘z’,legend.pos=2,pch=c(21:24),main=“”,cex=0.7)
mtext(side=1,line=2,“x”) ...