plot.logratioVariogramAnisotropy: Plot variogram lines of empirical directional logratio...

In gmGeostats: Geostatistics for Compositional Analysis

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) 

gmGeostats documentation built on April 18, 2023, 5:08 p.m.