PlotSNR: Plot proxy SNR (Figure 3 in Münch and Laepple, 2018)

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

Description

Plot the timescale dependence of proxy signal-to-noise ratios.

Usage

PlotSNR(
  spec,
  f.cut = FALSE,
  names = NULL,
  col = 1:length(spec),
  xlim = c(500, 2),
  ylim = c(0.05, 5),
  plt.ann = "default",
  xlab = NULL,
  ylab = NULL,
  xtm = NULL,
  ytm = NULL,
  xtl = NULL,
  ytl = NULL
)

Arguments

spec	a (named) list of signal-to-noise ratio data sets: each data set itself should be list containing at least a named element snr which is an object of class "spec" providing signal-to-noise ratios as a function of frequency. For Figure 3 in Münch and Laepple (2018) set spec to the output from PublicationSNR.
f.cut	Shall the spectra be cut at the cutoff frequency constrained by the diffusion correction strength? Defaults to FALSE.
names	an optional character vector of names of the proxy data sets. If NULL, the names of spec are used or, if not present, default names.
col	a numeric or character vector of colors to use for the plotting with length recycled to match length(spec).
xlim	the x limits (x1, x2) of the plot.
ylim	the y limits (y1, y2) of the plot.
plt.ann	if "default" use axis annotation as in Fig. 3 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

PublicationSNR

Examples

# Plot Figure 3 in Münch and Laepple (2018)
# (DML and WAIS oxygen isotope data sets):

# Load main spectral results
DWS <- WrapSpectralResults(
       dml1 = dml$dml1, dml2 = dml$dml2, wais = wais,
       diffusion = diffusion.tf,
       time.uncertainty = time.uncertainty.tf,
       df.log = c(0.15, 0.15, 0.1))

# Calculate the final signal-to-noise ratio spectra
SNR <- proxysnr:::PublicationSNR(DWS)

# Plot it
PlotSNR(SNR, f.cut = TRUE,
        names = c("DML", "WAIS"),
        col = c("black", "dodgerblue4"))

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