krige function for interpolation - methodology
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>