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.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