PlotTSDVariance: Plot Time-scale Dependent Variance Object.

In EarthSystemDiagnostics/psem: Proxy Spectral Error Model

Plot Time-scale Dependent Variance Object.

Description

Plot Time-scale Dependent Variance Object.

Usage

PlotTSDVariance(
  var.obj,
  components = c("Bioturbation", "Seasonal.bias", "Seasonal.bias.unc.",
    "Calibration.unc.", "Aliasing.seasonal", "Aliasing.stochastic",
    "Individual.variation", "Meas.error"),
  include = TRUE,
  include.constant.errors = FALSE,
  units = c("degC", "d18O"),
  thin.var.obj = TRUE
)

Arguments

var.obj	Object of class proxy.error.var, e.g. output from IntegrateErrorSpectra
components	Which error components to include or exclude
include	Include or exclude components named in components
include.constant.errors	Include the constant (nu = 0) error components
units	Are the units degC or d18O, used for axis labelling
thin.var.obj	If TRUE and the variance object contains more than 1000 frequencies these are thinned to 1000 frequencies equally spaced on a log-scale

Value

ggplot object

Examples

spec.pars <- GetSpecPars("Mg_Ca", delta_phi_c = pi)
spec.obj <- do.call(ProxyErrorSpectrum, spec.pars)

var.obj <- IntegrateErrorSpectra(spec.obj)
PlotTSDVariance(var.obj, units = "degC")
PlotTSDVariance(var.obj, units = "d18O")
PlotTSDVariance(var.obj, include = FALSE, c("Seasonal.bias", "Calibration.unc."))
PlotTSDVariance(var.obj, include.constant.error = TRUE)

EarthSystemDiagnostics/psem documentation built on March 3, 2024, 10:24 p.m.