IntegrateErrorSpectra | R Documentation |
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.
IntegrateErrorSpectra(
pes,
method = "sincfilter",
tau_smooth = NULL,
thin.tau_smooth = TRUE
)
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 |
timescale dependent variance object (dataframe)
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.