testBackground: Evaluate power spectrum false positive rates via Monte Carlo...

In astrochron: A Computational Tool for Astrochronology

Evaluate power spectrum false positive rates via Monte Carlo simulation

Description

This is a simulation tool to evaluate power spectrum false positive rates, the frequency distribution of the false positives, and the behavior of numerous "multiple correction" procedures, for a range of background estimation approaches that are implemented in Astrochron. The tool can be used to conduct surrogate analyses, alongside analysis of real data, to better understand the suitability of particular background estimation approaches. The resulting simulations are similar to those presented in Figure 3 of Meyers (2012) and Crampton et al. (PNAS).

Usage

testBackground(npts=1001,dt=5,noiseType="ar1",coeff=NULL,method="periodogramAR1",
               opt=NULL,demean=T,detrend=F,low=0,tbw=3,multi=F,iter=2000,output=F,
               genplot=F,verbose=T)

Arguments

npts	Number of points in simulated stratigraphic series (surrogates).
dt	Sampling interval for surrogates.
noiseType	Select "ar1" for AR1 noise surrogates, or "pwrLaw" for Power Law noise surrogates
coeff	AR1 coefficient (rho) or Power Law coefficient (beta) for surrogates.
method	Background estimation method: (1) "mtmAR1" (function mtm), (2) "mtmML96" (function mtmML96), (3) "lowspec" (function lowspec), (4) "mtmPL" (function mtmPL), (5) "periodogramPL" (function periodogram), (6) "periodogramAR1" (function periodogram)
opt	Method specific options. For mtmML96, this is medsmooth (see function mtmML96); for lowspec this is lowspan (see function lowspec); for periodogram this is percent cosine taper (see function cosTaper).
demean	Remove mean value from simulated surrogates? (T or F; this option does not apply to lowspec)
detrend	Remove linear trend from simulated surrogates? (T or F)
low	Remove long-term trend using a LOWESS smoother? Choose a value ranging from 0-1 (see function noLow). 0 = no long-term trend removal.
tbw	MTM time-bandwidth product. This option is ignored for methods 5 and 6.
multi	Evaluate a range of multiple-comparison tests too? (T or F)
iter	Number of iterations (surrogate series) for Monte Carlo simulation.
output	Output data frame? (T or F)
genplot	Generate summary plots? (T or F)
verbose	Verbose output? (T or F)

Details

The Monte Carlo simulations can utilize AR1 or Power Law noise surrogates. Background estimation approaches include conventional AR1, ML96, LOWSPEC and Power Law. The function also allows evaluation of common data detrending approaches (linear trend removal, LOWESS trend removal).

Note that MTM-ML96 conducts the Mann and Lees (1996; ML96) "robust red noise" analysis, with an improved median smoothing approach. The original Mann and Lees (1996) approach applies a truncation of the median smoothing window to include fewer frequencies near the edges of the spectrum; while truncation is required, its implementation in the original method often results in an "edge effect" that can produce excess false positive rates at low frequencies, commonly within the eccentricity-band (Meyers, 2012). To help address this issue, an alternative median smoothing approach is applied that implements Tukey's robust end-point rule and symmetrical medians (see the function mtmML96 for more details). This version of the ML96 algorithm was first implemented in Patterson et al. (2014).

See function multiTest for more information on the multiple comparison tests evaluated.

References

W.S. Cleveland, 1979, Locally weighted regression and smoothing scatterplots: Journal of the American Statistical Association, v. 74, p. 829-836.

J.S. Campton, S.R. Meyers, R.A. Cooper, P.M Sadler, M. Foote, D. Harte, 2018, Pacing of Paleozoic macroevolutionary rates by Milankovitch grand cycles: Proceedings of the National Academy of Sciences, doi:10.1073/pnas.1714342115.

M.E. Mann, and J.M. Lees, 1996, Robust estimation of background noise and signal detection in climatic time series, Clim. Change, 33, 409-445.

S.R. Meyers, 2012, Seeing Red in Cyclic Stratigraphy: Spectral Noise Estimation for Astrochronology: Paleoceanography, 27, PA3228, doi:10.1029/2012PA002307.

M.O. Patterson, R. McKay, T. Naish, C. Escutia, F.J. Jimenez-Espejo, M.E. Raymo, M.E., S.R. Meyers, L. Tauxe, H. Brinkhuis, and IODP Expedition 318 Scientists,2014, Response of the East Antarctic Ice Sheet to orbital forcing during the Pliocene and Early Pleistocene, Nature Geoscience, v. 7, p. 841-847.

D.J. Thomson, 1982, Spectrum estimation and harmonic analysis: IEEE Proceedings, v. 70, p. 1055-1096.

See Also

confAdjust,multiTest,lowspec, mtm, mtmML96, mtmPL, and periodogram

Examples

 ## Not run: 
# evaluate false positive rate for MTM-AR1 using AR1 surrogates
testBackground(noiseType="ar1",method="mtmAR1")

# evaluate false positive rate for MTM-AR1 using Power Law surrogates
testBackground(noiseType="pwrLaw",method="mtmAR1")
 
## End(Not run)

astrochron documentation built on Sept. 30, 2024, 9:14 a.m.