PlotArraySpectra: Plot proxy array spectra (Figure 1 in Münch and Laepple,...

In EarthSystemDiagnostics/proxysnr: Separate Signal and Noise in Climate Proxy Records

Description

Plot the spectral estimates from a spatial proxy core array.

Usage

PlotArraySpectra(
  spec,
  f.cutoff = NA,
  xlim = c(100, 2),
  ylim = c(0.005, 50),
  col = c("darkgrey", "black", "burlywood4"),
  col.sn = c("dodgerblue4", "firebrick4"),
  alpha.sn = 0.2,
  plt.ann = "default",
  xlab = NULL,
  ylab = NULL,
  xtm = NULL,
  ytm = NULL,
  xtl = NULL,
  ytl = NULL
)

Arguments

spec	output from ArraySpectra.
f.cutoff	optional frequency value to draw a vertical line symbolising a cutoff frequency.
xlim	the x limits (x1, x2) of the plot.
ylim	the y limits (y1, y2) of the plot.
col	a three-element vector with the colours for the individual spectra (col[1]), the mean (col[2]) and the stack spectrum (col[3]).
col.sn	a two-element vector with the colours for the approximate signal (col[1]) and noise shading (col[2]).
alpha.sn	opacity factor for the colours in col.sn within [0,1].
plt.ann	if "default" use axis annotation as in Fig. 1 of Münch and Laepple (2018). Since no other fixed annotation scheme is implemented, setting plt.ann to a different value will result in an error.
xlab	if not NULL use this specific x axis label to override the default setting.
ylab	if not NULL use this specific y axis label to override the default setting.
xtm	if not NULL use this vector of specific x axis tick mark positions to override the default setting.
ytm	if not NULL use this vector of specific y axis tick mark positions to override the default setting.
xtl	x axis tick mark labels. If xtm is not NULL, xtl is used for the labels parameter of the axis function: if xtl is set to NULL, the labels are determined automatically, else set it explicitly by specifying a vector of labels of the same length as xtm; see also ?axis.
ytl	equivalent to xtl for the y axis tick mark labels.

Author(s)

Thomas Münch

References

Münch, T. and Laepple, T.: What climate signal is contained in decadal- to centennial-scale isotope variations from Antarctic ice cores? Clim. Past, 14, 2053–2070, https://doi.org/10.5194/cp-14-2053-2018, 2018.

See Also

ArraySpectra

Examples

# Plot Figure 1 in Münch and Laepple (2018) (DML1 oxygen isotope data set):

PlotArraySpectra(ArraySpectra(dml$dml1, df.log = 0.12))

EarthSystemDiagnostics/proxysnr documentation built on Oct. 2, 2021, 3:03 p.m.