https://github.com/cran/gstat
Raw File
Tip revision: e0348cd08bd1aa1c5af65a89bd7ffcafe96f8fdb authored by Edzer J. Pebesma on 25 January 2006, 20:47:01 UTC
version 0.9-25
Tip revision: e0348cd
sim.R
library(gstat)
data(meuse)
set.seed(158229572)
new.locs <- data.frame(x = c(181170, 180310, 180205, 178673, 178770, 178270),
	y = c(333250, 332189, 331707, 330066, 330675, 331075))
krige(zinc ~ 1, ~ x + y, meuse, newdata = new.locs, 
		model = vgm(1.34e5, "Sph", 800, nug = 2.42e4), 
		block = c(40,40), nmax = 40, nsim = 10)
back to top