plot.mtm: Compute and plot the multitaper spectrum estimate

Description

Plots the multitaper spectral estimate and the multitaper F-test.

Usage

# S3 method for mtm
plot(x, jackknife = FALSE, Ftest = FALSE, ftbase = 1.01, siglines = NULL, ...)

Arguments

x: An object of the class mtm generated by spec.mtm.
jackknife: Boolean variable indicating if jackknife confidence intervals should be plotted, only applies if Ftest=FALSE.
Ftest: Boolean variable indicating if the multitaper harmonic F-test should be plotted instead of the spectrum.
ftbase: Lowest value to be plotted when the F-test is plotted. When Ftest = TRUE, max(ftestvalue, ftbase) is plotted.
siglines: Vector of significance values (as probabilities, 0.0 to 1.0) to plot as horizontal significance lines on the F-test plot.
...: Arguments to be passed to methods, such as graphical parameters (see 'par').

Details

The value log can be set to “yes” (default), “no”, or “dB” as in the function plot.spec.

References

Thomson, D.J (1982) Spectrum estimation and harmonic analysis. Proceedings of the IEEE Volume 70, number 9, pp. 1055--1096.

Percival, D.B. and Walden, A.T. (1993) Spectral analysis for physical applications Cambridge University Press.

Examples

Run this code

data(willamette)
resSpec <- spec.mtm(willamette, nw=4.0, k=8, Ftest=TRUE, plot=FALSE, deltat=1.0, dtUnits="month")
plot(resSpec)
plot(resSpec, Ftest=TRUE)
plot(resSpec, Ftest=TRUE, siglines=c(0.90, 0.99))

# with jackknife estimate
resSpec2 <- spec.mtm(willamette, nw=4.0, k=8, Ftest=TRUE, jackknife=TRUE, plot=FALSE,
                     deltat=1.0, dtUnits="month")
plot(resSpec2,jackknife=TRUE)