plot.logratioVariogramAnisotropy {gmGeostats}

Plot variogram lines of empirical directional logratio variograms

Description

Plots an "logratioVariogramAnisotropy" object in a series of panels, with each direction represented as a broken line.

Usage

## S3 method for class 'logratioVariogramAnisotropy'
plot(
  x,
  azimuths = colnames(x),
  col = rev(rainbow(length(azimuths))),
  type = "o",
  V = NULL,
  lty = 1,
  pch = 1:length(azimuths),
  model = NULL,
  figsp = 0,
  closeplot = TRUE,
  ...
)

Arguments

x

logratio variogram with anisotropy, i.e. object of class c("logratioVariogramAnisotropy", "logratioVariogram")

azimuths

which directions do you want to plot? default: all directions available

col

colors to be used for plotting

type

type of representation, see graphics::plot()

V

optionally, a matrix of logcontrasts, or else one of the following strings: "alr", "ilr" or "clr"; to produce a plot of the empirical variogram in the corresponding representation; default to variation-variograms

lty

style of the lines, potentially different for each directions

pch

symbols for the points, potentially different for each directions

model

eventually, variogram model to plot on top of the empirical variogram

figsp

spacing between the several panels, if desired

closeplot

boolean, should the plotting par() be returned to the starting values? (defaults to TRUE; see plot.gmCgram() for details)

...

additional graphical arguments, to be passed to the underlying graphics::matplot() function

Value

Nothing. The function is called to create a plot.

Examples

data("jura", package="gstat")
X = jura.pred[,1:2]
Zc = compositions::acomp(jura.pred[,7:9])
vg = logratioVariogram(data=Zc, loc=X, azimuth=c(0:3)*45, 
                         azimuth.tol=22.5)
plot(vg)