[R-sig-Geo] How to make directional h-scatterplots

Edzer Pebesma edzer.pebesma at uni-muenster.de
Fri Sep 10 07:52:48 CEST 2010


you can pass the alpha, beta, tol.hor and tol.ver parameters to hscat,
who will pass it on to variogram(). The problem is that the alpha
parameter, when set, is also passed to xyplot() and messes up the plot
there - I will choose a different name for that one.

On 09/10/2010 01:59 AM, Waichler, Scott R wrote:
> What is the best way to make directional h-scatterplots (lagged scatterplots)?  There is hscat() in gstat, but I don't see a way to limit point pair selection to certain directions.  I have a 3D dataset based on logs from 30 different wells, with vertical spacing between points and total extent both smaller than horizontal spacing and extent. 
> 
> Here is what I have figured out so far.  For the horizontal direction, I can select points in a given x-y plane before calling hscat().  For the vertical direction, I can do a de facto selection by limiting the maximum lag distance to something smaller than the minimum distance between my wells.  
> 
> My main interest for this is in 2 directions, horizontal and vertical, but it would be nice to specify angles as in the anis argument for vgm().
> 
> Thanks,
> Scott Waichler
> Pacific Northwest National Laboratory
> scott.waichler at pnl.gov
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> 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.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de



More information about the R-sig-Geo mailing list