back with cutoff/variogram question
Hi! This is from the gstat variogram. I tried to track this down via the vario.c function, but I'm having some trouble with it.
From: b.rowlingson at gmail.com [b.rowlingson at gmail.com] on behalf of Barry Rowlingson [b.rowlingson at lancaster.ac.uk]
Sent: Sunday, November 24, 2013 3:08 AM
To: Hodgess, Erin
Cc: r-sig-geo at r-project.org
Subject: Re: [R-sig-Geo] back with cutoff/variogram question
Sent: Sunday, November 24, 2013 3:08 AM
To: Hodgess, Erin
Cc: r-sig-geo at r-project.org
Subject: Re: [R-sig-Geo] back with cutoff/variogram question
On Sun, Nov 24, 2013 at 5:46 AM, Hodgess, Erin <HodgessE at uhd.edu> wrote: > Hello again. > > I'm back with the variogram question. Here is the info from the help file: > > cutoff: spatial separation distance up to which point pairs are > included in semivariance estimates; as a default, the length > of the diagonal of the box spanning the data is divided by > three. > > By the box spanning the data, does that mean the bounding box, please? If so, wouldn't that just be 2, please? There's quite a few 'variogram' functions in various R packages so it might help if you told us which one! I should think that "the box spanning the data" is indeed the bounding box. I don't see what you think would be '2' in this case though, unless the diagonal length of the bounding box happened to be 6. For data that included points at (0,0) and (1,1) the diagonal would have length sqrt(2) and so the cutoff default should be sqrt(2)/3. Note you have to actually have points at those coordinates and not just generate 1000 random points in the unit square. A few simple experiments with simulated point data sets with various bounding boxes should confirm this, or you could just look at the code? Barry