swh:1:snp:33a53053e50f7abe7d281cc0c803be827debf4a3
Tip revision: 1402a3d4a695b8bf8dc225e4d787d46a21fb6ff6 authored by Edzer Pebesma on 20 September 2010, 07:41:03 UTC
version 0.9-71
version 0.9-71
Tip revision: 1402a3d
blockkr.R
library(gstat)
data(meuse)
coordinates(meuse) = c("x", "y")
new.locs <- SpatialPoints(data.frame(
x = c(181170, 180310, 180205, 178673, 178770, 178270),
y = c(333250, 332189, 331707, 330066, 330675, 331075)))
krige(zinc ~ 1, meuse, new.locs, vgm(1.34e5, "Sph", 800, nug = 2.42e4),
block = c(40,40), nmax = 40)
new.locs <- SpatialPoints(data.frame(x = c(181170), y = c(333250)))
disc = c(-15,-5,5,15)
block.irreg <- data.frame(expand.grid(x = disc, y = disc))
block.irreg
# first disable default Gaussian quadrature used for block integration, by
# setting nblockdiscr explicitly:
krige(zinc ~ 1, meuse, new.locs, model = vgm(1.34e5, "Sph", 800, nug = 2.42e4),
block = c(40,40), nmax = 40, set = list(nblockdiscr=4))
# now pass the same block discretization as block.irreg:
krige(zinc ~ 1, meuse, new.locs, vgm(1.34e5, "Sph", 800, nug = 2.42e4),
block = block.irreg, nmax = 40)
# check weights argument:
block.irreg <- data.frame(expand.grid(x = disc, y = disc), weights = rep(1/16, 16))
krige(zinc ~ 1, meuse, new.locs, vgm(1.34e5, "Sph", 800, nug = 2.42e4),
block = block.irreg, nmax = 40)