IntegrateErrorSpectra: Integrate Error Spectra

In EarthSystemDiagnostics/psem: Proxy Spectral Error Model

Integrate Error Spectra

Description

Integrates error spectra to get the variance(s) and optionally filters to simulate the case where a running mean has been passed over a proxy record.

Usage

IntegrateErrorSpectra(
  pes,
  method = "sincfilter",
  tau_smooth = NULL,
  thin.tau_smooth = TRUE
)

Arguments

pes	Object of class proxy.error.spec, e.g. output from ProxyErrorSpectrum
method	sinc filter is the only currently supported method
tau_smooth	temporal resolution to which the proxy timeseries is smoothed. If NULL, variance is returned at all odd multiples of delta_t
thin.tau_smooth	If TRUE and the vector of tau_smooth values is longer than 1000, these are thinned to 1000 values equally spaced on a log-scale

Value

timescale dependent variance object (dataframe)

Examples

spec.pars <- GetSpecPars("Mg_Ca", tau_p = 1 / 12, phi_c = 0, seas.amp = 4, T = 100 * 101)
spec.obj <- do.call(ProxyErrorSpectrum, spec.pars)
PlotSpecError(spec.obj, show.low.power.panel = FALSE)
var.obj <- IntegrateErrorSpectra(spec.obj)
PlotTSDVariance(var.obj)

var.obj <- IntegrateErrorSpectra(spec.obj, tau_smooth = 1000)

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