R: Compute and plot the multitaper spectrum estimate

plot.mtm {multitaper} | R Documentation

## 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)

[Package multitaper version 1.0-13 Index]