https://github.com/cran/gstat
Raw File
Tip revision: 76883c10118b75b3883c3bc21e980beb0bd5783d authored by Edzer Pebesma on 19 February 2013, 15:39:26 UTC
version 1.0-16
Tip revision: 76883c1
vdist.R
library(gstat)

data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
gridded(meuse.grid) = ~x+y

mg = meuse.grid
gridded(mg) = FALSE
mg= mg[1500,]
krige(log(zinc)~1,meuse,mg,vgm(1, "Exp", 300, anis=c(0,0.01)),
	vdist=FALSE, maxdist=1000,nmax=10)
krige(log(zinc)~1,meuse,mg,vgm(1, "Exp", 300, anis=c(0,0.01)),
	vdist=TRUE, maxdist=1000,nmax=10)
back to top