Im 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]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
