plot.mtm {multitaper} | R Documentation |
Plots the multitaper spectral estimate and the multitaper F-test.
## S3 method for class 'mtm' plot(x, jackknife = FALSE, Ftest = FALSE, ftbase = 1.01, siglines = NULL, ...)
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'). |
The value log can be set to “yes” (default), “no”, or “dB” as in the function plot.spec.
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.
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)