logratioVariogram {gmGeostats} | R Documentation |
Empirical logratio variogram calculation
Description
Calculation of an empirical logratio variogram (aka variation-variogram)
Usage
logratioVariogram(data, ...)
Arguments
data |
a data container (of class "gmSpatialModel") or a composition (of class "acomp") |
... |
extra arguments for generic functionality |
Value
The output of this function depends on the number of azimuth
values provided:
if it is one single value (or if you explicitly call logratioVariogram.default
) the result is
a list of class "logratioVariogram" with the following elements
vg: A nbins x D x D array containing the logratio variograms
h: A nbins x D x D array containing the mean distance the value is computed on.
n: A nbins x D x D array containing the number of nonmissing pairs used for the corresponding value. If
azimuth
is a vector of directions, then the result is a matrix-like list of "logratioVariogram" objects. Each element of the mother list (or column of the matrix) is the variogram of one direction. The output has a compound class c("logratioVariogramAnisotropy", "logratioVariogram").
Examples
data("jura", package="gstat")
X = jura.pred[,1:2]
Zc = compositions::acomp(jura.pred[,7:13])
vg = logratioVariogram(data=Zc, loc=X)
class(vg)
summary(vg)
vg = logratioVariogram(data=Zc, loc=X, azimuth=c(0,90))
class(vg)
summary(vg)