ArraySpectra: Spectral estimates of core array

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

Description

ArraySpectra calculates all relevant spectral estimates from a given core array of N proxy records. The spectral estimates can be smoothed in logarithmic space.

Usage

ArraySpectra(cores, res = 1, neff = length(cores), df.log = NULL, ...)

Arguments

cores	a list of the proxy data from the core array. Each component is expected to be a numeric vector of the proxy values of a common length.
res	the sampling (e.g., temporal) resolution of the proxy data; determines the frequency axis of the spectral estimates.
neff	the effective number of records (e.g. to account for an expected spatial correlation of the local noise). Per default, no spatial correlation is assumed and neff is set to the number of proxy records (the length of cores).
df.log	width of the Gaussian kernel in logarithmic frequency units to smooth the spectral estimates; NULL (the default) suppresses smoothing.
...	additional parameters which are passed to the spectral estimation function SpecMTM.

Details

The spectral estimates are calculated using Thomson’s multitaper method with three windows with linear detrending before analysis (see SpecMTM). Each spectral result is returned as an object of class "spec" with the minimum elements freq and spec.

Value

A list of the following components:

N:

the number of (effective) proxy records of the core array.

single:

a list of length N with the spectra of each individual proxy record.

mean:

the mean spectrum across all individual spectra.

stack:

the spectrum of the average proxy record in the time domain ("stacked record").

Author(s)

Thomas Münch

See Also

SpecMTM

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