Multitaper Method (MTM) Spectral Analysis with Power Law (1/f) fit
mtmPL(dat,tbw=3,ntap=NULL,padfac=5,demean=T,detrend=F,siglevel=0.9,flow=NULL,fhigh=NULL,
output=0,CLpwr=T,xmin=0,xmax=Nyq,pl=1,sigID=F,genplot=T,verbose=T)Stratigraphic series for MTM spectral analysis. First column should be location (e.g., depth), second column should be data value.
MTM time-bandwidth product.
Number of DPSS tapers to use. By default, this is set to (2*tbw)-1.
Pad with zeros to (padfac*npts) points, where npts is the original number of data points.
Remove mean from data series? (T or F)
Remove linear trend from data series? (T or F)
Significance level for peak identification.
lowest frequency to include in 1/f fit
highest frequency to include in 1/f fit
What should be returned as a data frame? (0=nothing; 1=spectrum + CLs + AR1 fit; 2=sig peak freqs; 3=sig peak freqs + prob; 4=all)
Plot AR(1) noise confidence levels on power spectrum? (T or F)
Smallest frequency for plotting.
Largest frequency for plotting.
Plot logarithm of spectral power (1) or linear spectral power (2)?
Identify signficant frequencies on power and probabilty plots? (T or F)
Generate summary plots? (T or F)
Verbose output? (T or F)
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.
spec.mtm, lowspec, and periodogram