View source: R/compositionsCompatibility.R
logratioVariogram | R Documentation |
Calculation of an empirical logratio variogram (aka variation-variogram)
logratioVariogram(data, ...)
data |
a data container (of class "gmSpatialModel") or a composition (of class "acomp") |
... |
extra arguments for generic functionality |
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").
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.