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

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)

[Package gmGeostats version 0.11.3 Index]