Reply to comment

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")

Reply

CAPTCHA
This question is for testing whether you are a human visitor and to prevent automated spam submissions.
Question text provided by textcaptcha.com