warning: Parameter 2 to gmap_gmap() expected to be a reference, value given in /home/spatiala/public_html/book/includes/module.inc on line 497.

Fig. 5.13: Four simulations of liming requirements (indicator variable) using ordinary kriging.

Tags:
Fig. 5.13: Four simulations of liming requirements (indicator variable) using ordinary kriging.
# geostatistical simulations (lime)
lime.ovgm <- fit.variogram(variogram(I(lime == 1) ~ 1, meuse), vgm(1, "Sph", 800, 1))
plot(variogram(I(lime == 1) ~ 1, meuse), lime.ovgm)
lime.sim <- krige(I(lime == 1) ~ 1, meuse, meuse.grid, lime.ovgm, nsim=4, indicators = TRUE, nmax=40)
spplot(lime.sim, col.regions=grey(c(rev(seq(0,0.9,0.05)), 0)), main="OK simulations of lime")