# Question on lags in h-Scattergrams Classic List Threaded 2 messages Reply | Threaded
Open this post in threaded view
|

## Question on lags in h-Scattergrams

 My question is not R specific, but about h-scattergrams as described in Pierre Goovaerts's "Geostatistics for Natural Resource Evaluation." No matter how many times I read these sections I still cannot grok the relationship between delta-h (which he sets at 100m in the discussions beginning in Section 2.3.2 on page 25) and the examples in following sections.    For example, on page 27, Figure 2.8 shows Cd concentrations at 6 classes of distances, 'h' with distance tolerance 100m (the delta-h). Where does the 100m distance show up in the plots with h=43m through h=1024m? This is what confuses me.    Learning spatial statistics on my own (all my stat classes in grad school were biometry focused) leaves me lacking a local resource to ask. TIA, Rich _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

## Re: Question on lags in h-Scattergrams

 On 04/01/2011 10:22 PM, Rich Shepard wrote: >   My question is not R specific, but about h-scattergrams as described in > Pierre Goovaerts's "Geostatistics for Natural Resource Evaluation." No > matter how many times I read these sections I still cannot grok the > relationship between delta-h (which he sets at 100m in the discussions > beginning in Section 2.3.2 on page 25) and the examples in following > sections. > >   For example, on page 27, Figure 2.8 shows Cd concentrations at 6 classes > of distances, 'h' with distance tolerance 100m (the delta-h). Where does > the > 100m distance show up in the plots with h=43m through h=1024m? This is what > confuses me. It seems to me that in figure 2.8 the average distance of all point pairs in a particular lag class is shown, but the exact class boundaries are not given (delta h = 100, but the h values in h +/- delta h are not). Looking at the average values, I would assume lag intervals 0-100, 100-300, 300-500 etc. If I try to reproduce this, as in library(gstat) data(jura) example(jura) hscat(Cd~1, jura.pred, c(0,100,300,500,700,900,1100)/1000, variogram.alpha=90, tol.hor=22.5, asp="iso") I'm getting quite similar, but not identical results, so some of my assumptions must be false. Maybe Pierre can help us out. Note that the intervals in this plot are in km. With best regards, >   Learning spatial statistics on my own (all my stat classes in grad school > were biometry focused) leaves me lacking a local resource to ask. > > TIA, > > Rich > > _______________________________________________ > R-sig-Geo mailing list > [hidden email] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo-- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.dehttp://www.52north.org/geostatistics      [hidden email] _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo