cross_spectrum: Calculate the cross-spectrum of two timeseries

In kitagawa: Spectral Response of Water Wells to Harmonic Strain and Pressure Signals

Description

Calculate the cross-spectrum of two timeseries

Usage

cross_spectrum(x, ...)

## S3 method for class 'mts'
cross_spectrum(x, ...)

## Default S3 method:
cross_spectrum(
  x,
  y,
  k = 10,
  samp = 1,
  q,
  adaptive = FALSE,
  verbose = FALSE,
  ...
)

Arguments

x	numeric; timeseries
...	additional arguments to pspectrum
y	numeric; timeseries. if missing, assumed to be column no. 2 in x
k	integer; the number of sine tapers, unless this is NULL; in the latter case a Welch-based spectrum is calculated rather than a multitaper spectrum. There are distinct advantages and disadvantages to either of these.
samp	numeric; the sampling rate (e.g., deltat) of the data; must be the same for x and y
q	numeric; the probability quantile [0,1] to calculate coherence significance levels; if missing, a pre-specified sequence is included. This is will be ignored for Welch-based spectra (see k).
adaptive	logical; should adaptive multitaper estimation be used?
verbose	logical; should messages be printed?

Examples

require(stats)
require(psd)

n <- 1000
ramp <- seq_len(n)
parab <- ramp^2

set.seed(1255)
X <- ts(rnorm(n) + ramp/2)
Y <- ts(rnorm(n) + ramp/10 + parab/100)

# Calculate the multitaper cross spectrum
csd <- cross_spectrum(X, Y, k=20)

kitagawa documentation built on July 2, 2020, 1:47 a.m.