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.waich...@pnl.gov > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@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.pebe...@wwu.de _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
