Message-ID: <4C63057D.2080301@uni-muenster.de>
Date: 2010-08-11T20:18:05Z
From: Edzer Pebesma
Subject: Directional variograms with 3D data
In-Reply-To: <A632B3F93D861843A3C966F3962C3F8401080DD39972@EMAIL05.pnl.gov>
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.waichler at pnl.gov
> 509-372-4423, 509-341-4051 (cell)
>
> _______________________________________________
> 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