This is expected -- tol.vert is ignored because it should be tol.ver (see docs).
The reason why you don't (can't!) get a warning is that it get absorbed by the ... , which just passes on things not known (and hence accepts anything). On 08/11/2010 06:47 AM, Waichler, Scott R wrote: > Hello, > > I have a dataset with coordinates in x, y, z. I would like to generate > directional variograms, including ones for distinct horizontal and vertical > directions, and fit models to them. Apparently geoR does not accept 3D > coordinates, so I am left with (?) gstat. I have trouble understanding and > applying the 5-parameter version of anis in gstat's vgm(), so I thought I > would try generating one direction at a time, as below with a data frame > called "han": > > th.vg1 <- variogram(th ~ 1, han, alpha=0, tol.hor = 90, beta=0, tol.vert=30) > # within 30 deg of horizontal > th.vg2 <- variogram(th ~ 1, han, alpha=0, tol.hor = 90, beta=90, tol.vert=30) > # within 30 deg of vertical > > th.vg1 and th.vg2 are identical; both are omnidirectional. What do I need to > do to to get distinct horizontal and vertical directions? > > Scott Waichler > Pacific Northwest National Laboratory > P.O. Box 999, Richland, WA 99352 > scott.waich...@pnl.gov > 509-372-4423, 509-341-4051 (cell) > > _______________________________________________ > 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
