https://github.com/cran/gstat
Raw File
Tip revision: bd635dd8d73a2ff20b2ca8c7bd579ad62d0e19b0 authored by Edzer Pebesma on 15 March 2011, 17:40:34 UTC
version 0.9-79
Tip revision: bd635dd
rep.R
library(gstat)
data(meuse)
data(meuse.grid)
coordinates(meuse) = ~x + y
coordinates(meuse.grid) = ~x + y

# Variogram log Zn

lzn.vgm = variogram(log(zinc) ~ 1, meuse)
lzn.fit = fit.variogram(lzn.vgm, model = vgm(1, "Sph", 900, 1))

#Conditional simulation

nsim = 100
lzn.sim = krige(log(zinc) ~ 1, meuse, meuse.grid, model = lzn.fit,
	nmax = 30, nsim = nsim)

# Variogram of all relizations

m = out = list()
for (i in 1:nsim) {
	s = paste("sim", i, sep="")
	f = as.formula(paste(s, "~1"))
	v = variogram(f, lzn.sim)
	v$id = rep(s, nrow(v))
	out[[s]] = v
	m[[s]] = fit.variogram(v, lzn.fit)
}
plot(do.call(rbind, out), m, layout=c(10,10), skip = FALSE, 
	scales = list(y = list(relation = "same")))
back to top