William Thayer wrote: > > I am looking for references that address estimating the mean of an > irregular block using kriging. > [...] > My questions are: [...] 2) > During some testing of Gstat, I have noticed that when I increase the area > of the irregular block beyond approximately one half of the site, Gstat > produces an estimate that is approximately three times the mean of the > observations located within the area. Has anyone encountered this with > Gstat before? FYI: The confusion arose around the definition of a search neighbourhood. Gstat calculates neigbhourhoods relative to the block centre, so a search radius is not defined relative to the block edge. In this case, the (few) data in the search neighbourhood contained one very large outlier. -- Edzer
