I’m using geoR to interpolate hydraulic conductivities. I constructed the
variogram using the log transform as follows:

dummy.var<-variog(dummy.geo,max.dist=40000,uvec=25,lambda=0)   # Then I
fitted a spherical model using eyefit parameters as starting parameters

dummy.param<-eyefit(dummy.var)

param1<-variofit(vario=dummy.var, ini.cov.pars=dummy.param[[1]]$cov.pars,

cov.model = dummy.param[[1]]$cov.model,nugget=dummy.param[[1]]$nugget)

spatPara<-krige.control(type.krige="ok",obj.model=param1)

# I then interpolate across the map

InterMap<-krige.conv(dummy.geo,loc=loc,krige=spatPara)


Ma question is why are my interpolated values higher than the actual
measurements. I expected the interpolation to honour measurement points at
least. Am I missing something?



Thanks

	[[alternative HTML version deleted]]

