paleodata_scaling: Function to compute scaling of the spectrum following the...

View source: R/Proxytools_tools.R

paleodata_scalingR Documentation

Function to compute scaling of the spectrum following the multi-taper methodology applied in Laepple and Huybers 2014, Rehfeld et al. 2018

Description

Function to compute scaling of the spectrum following the multi-taper methodology applied in Laepple and Huybers 2014, Rehfeld et al. 2018

Usage

paleodata_scaling(
  xin,
  freq.start = NULL,
  freq.end = NULL,
  target = "raw",
  interpolation = TRUE,
  interpolation_type = "spectral",
  interpolation_dates = NULL,
  transformation = FALSE,
  transformation_type = "normalize",
  detrend = TRUE,
  df_log = 0.05,
  bLog = FALSE
)

Arguments

xin

Proxytibble with proxy data in 'zoo::zoo' format, or irregular time series object ('zoo::zoo'), xin can be multivariate

freq.start

Vector containing the start frequencies of the fitting interval(s)

freq.end

Vector containing the start frequencies of the fitting interval(s)

target

either "raw" or "logsmooth", default "raw"

interpolation

Logical: should interpolation be applied?

interpolation_type

Parameters for interpolation (see paleodata_interpolation)

interpolation_dates

Parameters for interpolation (see paleodata_interpolation). If NULL the interpolation is set to the mean temporal resolution of each proxy time series.

transformation

Logical: should transformation be applied?

transformation_type

Parameters for transformation (see paleodata_transformation)

detrend

Logical: should time series be linearly detrend before computing the spectrum

df_log

Width of the log-smoother in log units

bLog

TRUE: average in the log space of the power, FALSE: arithmetic average

Value

List with scaling coefficient, standard deviation and spectra (object description in SpecMTM and LogSmooth functions of PaleoSpec package)

See Also

paleodata_spectrum (from 'PTBoxProxytools') for computation of spectra

SlopeFit (from 'PaleoSpec') for computation of spectral slope

Examples

# Load ice core example data
library(PTBoxProxydata)
mng <- ProxyDataManager()
icecoredata <- load_set(mng,'icecore_testset',zoo_format = 'zoo')
# compute scaling
testout <- icecoredata %>% paleodata_scaling(., target="raw")
print(testout)


paleovar/ptboxproxytools documentation built on June 9, 2025, 1:40 a.m.