zijderveld: Draws a Zijderveld plot

Description

Draws a Zijderveld plot: it projects 3D points (having declination, inclination and intensity) in 2D, horizontally and vertically.

Usage

zijderveld(
  dec,
  inc,
  int,
  xh = "WE",
  xv = xh,
  centre = F,
  xlim = NA,
  ylim = NA,
  unit = NA,
  xlab = "",
  ylab = "",
  labels = NA,
  nlabels = 1,
  h = list(pch = 19),
  v = list(pch = 21, bg = "white"),
  f = list(pch = 21, bg = "white", cex = 1.5),
  t = list(pos = 3, offset = 0.5),
  l = list(),
  anchored = T,
  style = "branches",
  tcl = 0.2,
  orientation = TRUE,
  scientific = NA,
  decimals = 10,
  add = FALSE
)

Arguments

dec: declination of the data; it is the angle from the north taken on an horizontal plane. It is measured clockwise from North and ranges from 0 to 360° (Tauxe 2010). Values outside this range are corrected by incfix().
inc: inclination of the data; it is the angle from the horizontal, is positive downward, and ranges from +90° for straight down to -90° for straight up (Tauxe, 2010). Values outside this range are corrected by incfix().
int: intensity of the data.
xh: orientation of the x axis for the horizontal points: can be 'SN' or 'WE'.
xv: orientation of the x axis for the horizontal points: can be 'SN', 'WE' or 'modified' (for the latter the horizontal projection of the vector given by the square root of the addition of the squared horizontal components).
centre: logical, whether the [0,0] point should be in the centre of the plot. Is ignored if xlim and/or ylim are defined.
xlim, ylim: the x and y minimal limits. The actual limits can change to keep a x/y ratio of 1.
unit: the tick interval.
xlab, ylab: the titles for the axes.
labels: a character vector of labels to add to each point.
nlabels: the number of labels to skip (for clarity).
h, v, f, t, l: list of graphical parameters to feed the graphical functions: h, v and f are fed to points() for the horizontal, vertical and first points respectively; t is fed to the text() for the labels and l is fed to lines() for the lines joining each horizontal and vertical points. See ?points, ?text and ?lines help page for the possible arguments. See ?merge_list for further information.
anchored: logical, whether the lines should be anchored to the [0,0] point.
style: the style of the plot: 'branches', 'box0', 'box1', or 'box2'. The boxes are advised when zooming using xlim and/or ylim.
tcl: The length of tick marks (see par() help page).
orientation: logical, whether to add captions indicating the orientation of the plot.
scientific: logical or NA, whether have scientific notation (e.g. -1.0E-06) or not (e.g. 0.00015). If NA, R will be left only judge.
decimals: the number of decimals if scientific is T or F. Having not enough decimals can lead to override the unit parameter, but the tick labels will be correctly aligned.
add: logical, whether to add the plot to an existing plot.

Details

By default horizontal projection is made of black points, vertical of white points.

References

Tauxe, L., 2010. Essentials of Paleomagnetism. University of California Press.

Examples

Run this code

zd <- zeq_example

ori <- par()$mfrow

par(mfrow = c(1,2))

zijderveld(dec = zd$Dec, inc = zd$Inc, int = zd$Int,
           xh = "WE", unit = 10^-5)

zijderveld(dec = zd$Dec, inc = zd$Inc, int = zd$Int,
           style = "box1", scientific = FALSE, decimals = 5,
           labels = zd$Treat, nlabels = 2)

par(mfrow = ori)