DIMITRIS KARAKOSTIS
2012-Dec-17 16:06 UTC
[R] How to make Ordinary Kriging using gstat predict?
Hi, I am new in R and trying to implement an algorithm which makes ordinary kriging by using gstat library and the predict method. I use the predict method as following: First, I create an object g: g <- gstat(id="tec", formula=TEC ~ 1, data=data) ## Create gstat object called g And then I use this object in the predict. p <- predict.gstat(g, model=mod, newdata=predGrid, nsim=0) My problem is that instead of Ordinary kriging, when I run the algorithm I get: Inverse distance weighted interpolation. Why is that? What am I missing or doing wrong? Thanks in advance Dimitris [[alternative HTML version deleted]]
> -----Original Message----- > My problem is that instead of Ordinary kriging, when I run > the algorithm I get: Inverse distance weighted interpolation. > Why is that? What am I missing or doing wrong?The gstat manual at http://www.gstat.org/gstat.pdf says on p16 that "When no variograms are specified, inverse distance weighted interpolation is the default action (Fig. 2.1, example [6.3]). When variograms are specified the default prediction method is ordinary kriging Journel and Huijbregts (1978); Cressie (1993) (example [6.4] and example [6.8])." It looks like reading that manual may be useful ... S Ellison ******************************************************************* This email and any attachments are confidential. Any use...{{dropped:8}}