mtm | R Documentation |
Multitaper method (MTM) spectral analysis (Thomson, 1982)
mtm(dat,tbw=3,ntap=NULL,padfac=5,demean=T,detrend=F,siglevel=0.9,ar1=T,output=0,
CLpwr=T,xmin,xmax,pl=1,sigID=T,genplot=T,verbose=T)
dat |
Stratigraphic series for MTM spectral analysis. First column should be location (e.g., depth), second column should be data value. |
tbw |
MTM time-bandwidth product. |
ntap |
Number of DPSS tapers to use. By default, this is set to (2*tbw)-1. |
padfac |
Pad with zeros to (padfac*npts) points, where npts is the original number of data points. |
demean |
Remove mean from data series? (T or F) |
detrend |
Remove linear trend from data series? (T or F) |
siglevel |
Significance level for peak identification. (0-1) |
ar1 |
Estimate conventional AR(1) noise spectrum and confidence levels? (T or F) |
CLpwr |
Plot AR(1) noise confidence levels on power spectrum? (T or F) |
output |
What should be returned as a data frame? (0=nothing; 1= power spectrum + harmonic CL + AR1 CL + AR1 fit + 90%-99% AR1 power levels (ar1 must be set to TRUE to output AR model results); 2=significant peak frequencies; 3=significant peak frequencies + harmonic CL; 4=internal variables from spec.mtm). Option 4 is intended for expert users, and should generally be avoided. |
xmin |
Smallest frequency for plotting. |
xmax |
Largest frequency for plotting. |
pl |
Power spectrum plotting: (1) linear frequency-log spectral power, (2) linear frequency-linear spectral power (3) log frequency-log spectral power, (4) log frequency-linear spectral power |
sigID |
Identify significant frequencies on power and probabilty plots? (T or F) |
genplot |
Generate summary plots? (T or F) |
verbose |
Verbose output? (T or F) |
If ar1=T, candidiate astronomical cycles are identified via isolation of those frequencies that achieve the required (e.g., 90 percent) "red noise" confidence level and MTM harmonic F-test confidence level. Allowance is made for the smoothing inherent in the MTM power spectral estimate as compared to the MTM harmonic spectrum. That is, an F-test peak is reported if it achieves the required MTM harmonic confidence level, while also achieving the required red noise confidence level within +/- half the power spectrum bandwidth resolution. One additional criterion is included to further reduce the false positive rate, a requirement that significant F-tests must occur on a local power spectrum high, which is parameterized as occurring above the local red noise background estimate. See Meyers (2012) for futher information.
S.R. Meyers, 2012, Seeing Red in Cyclic Stratigraphy: Spectral Noise Estimation for Astrochronology: Paleoceanography, 27, PA3228, doi:10.1029/2012PA002307.
Rahim, K.J. and Burr W.S. and Thomson, D.J., 2014, Appendix A: Multitaper R package in "Applications of Multitaper Spectral Analysis to Nonstationary Data", PhD diss., Queen's Univieristy, pp 149-183. http://hdl.handle.net/1974/12584
Thomson, D. J., 1982, Spectrum estimation and harmonic analysis, Proc. IEEE, 70, 1055-1096, doi:10.1109/PROC.1982.12433.
eha
, lowspec
, mtmAR
, mtmML96
, periodogram
, and spec.mtm
# generate example series with periods of 400 ka, 100 ka, 40 ka and 20 ka
ex = cycles(freqs=c(1/400,1/100,1/40,1/20),start=1,end=1000,dt=5)
# add AR1 noise
noise = ar1(npts=200,dt=5,sd=.5)
ex[2] = ex[2] + noise[2]
# MTM spectral analysis, with conventional AR1 noise test
pl(1,title="mtm")
mtm(ex,ar1=TRUE)
# compare to ML96 analysis
pl(1, title="mtmML96")
mtmML96(ex)
# compare to analysis with LOWSPEC
pl(1, title="lowspec")
lowspec(ex)
# compare to amplitudes from eha
pl(1,title="eha")
eha(ex,tbw=3,win=1000,pad=1000)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.