gmGeostats (version 0.10-6)

logratioVariogram: 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

# NOT RUN {
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)
# }

Run the code above in your browser using DataCamp Workspace