[R-sig-Geo] krige function for interpolation - methodology

Edzer Pebesma edzer.pebesma at uni-muenster.de
Fri Mar 7 11:32:55 CET 2014

Please do read documentation before you post. This list is not meant to
remind you it is there.

Look up parameter maxdist in ?krige.

On 03/07/2014 12:14 AM, Moshood Agba Bakare wrote:
> Hi,
> I used the krige function of gstat to obtain predicted value for the grid
> points using block ordinary kriging. I defined my grid  spacing of 100m
> with block size of 10m x 10m using 100 nearest obervation in the local
> neighborhood of unsampled locations for prediction.
> I do not understand how this krige function will capture the number of
> nearest observations in the local neighborhood for prediction.
>  Could anyone please explain how it will capture the 100 obervations for
> the prediction? does this function define the search radius to capture
> nearest observation for prediction? Even though I define my nmax to be 100
> (nmax=100), is it possible that number of observations less than hundred
> can be used if the unsampled location is within sparse region of the field?
> prd_ok <- krige(id="yield",yield ~ 1, canmod.sp, newdata = grd,
> model=exp.mod,
> block=c(10,10),nmax=100)
> Thanks
> Moshood

Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251
83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 490 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20140307/f78d1a1d/attachment.bin>

More information about the R-sig-Geo mailing list