spec.ar | R Documentation |
Fits an AR model to x
(or uses the existing fit) and computes
(and by default plots) the spectral density of the fitted model.
spec.ar(x, n.freq, order = NULL, plot = TRUE, na.action = na.fail, method = "yule-walker", ...)
x |
A univariate (not yet:or multivariate) time series or the
result of a fit by |
n.freq |
The number of points at which to plot. |
order |
The order of the AR model to be fitted. If omitted, the order is chosen by AIC. |
plot |
Plot the periodogram? |
na.action |
|
method |
|
... |
Graphical arguments passed to |
An object of class "spec"
.
The result is returned invisibly if plot
is true.
Some authors, for example Thomson (1990), warn strongly that AR spectra can be misleading.
The multivariate case is not yet implemented.
Thompson, D.J. (1990). Time series analysis of Holocene climate data. Philosophical Transactions of the Royal Society of London Series A, 330, 601–616. \Sexpr[results=rd,stage=build]{tools:::Rd_expr_doi("10.1098/rsta.1990.0041")}. https://www.jstor.org/stable/53609.
Venables, W.N. and Ripley, B.D. (2002) Modern Applied Statistics with S. Fourth edition. Springer. (Especially page 402.)
ar
, spectrum
.
require(graphics) spec.ar(lh) spec.ar(ldeaths) spec.ar(ldeaths, method = "burg") spec.ar(log(lynx)) spec.ar(log(lynx), method = "burg", add = TRUE, col = "purple") spec.ar(log(lynx), method = "mle", add = TRUE, col = "forest green") spec.ar(log(lynx), method = "ols", add = TRUE, col = "blue")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.