plot.mtm {multitaper}

Compute and plot the multitaper spectrum estimate

Description

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

Usage

## S3 method for class '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.

See Also

spec.mtm and plot.spec

Examples

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)