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

Edzer Pebesma edzer.pebesma at uni-muenster.de
Fri Sep 10 08:19:14 CEST 2010


In response to this response, I uploaded gstat_0.9-70 to CRAN.

On 09/10/2010 07:52 AM, Edzer Pebesma wrote:
> 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