\HeaderA{spec.ar}{Estimate Spectral Density of a Time Series from AR Fit}{spec.ar}
\keyword{ts}{spec.ar}
\begin{Description}\relax
Fits an AR model to \code{x} (or uses the existing fit) and computes
(and by default plots) the spectral density of the fitted model.
\end{Description}
\begin{Usage}
\begin{verbatim}
spec.ar(x, n.freq, order = NULL, plot = TRUE, na.action = na.fail,
        method = "yule-walker", ...)
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{x}] A univariate (not yet:or multivariate) time series or the
result of a fit by \code{\LinkA{ar}{ar}}.
\item[\code{n.freq}] The number of points at which to plot.
\item[\code{order}] The order of the AR model to be fitted.  If omitted,
the order is chosen by AIC.
\item[\code{plot}] Plot the periodogram?
\item[\code{na.action}] \code{NA} action function.
\item[\code{method}] \code{method} for \code{ar} fit.
\item[\code{...}] Graphical arguments passed to \code{\LinkA{plot.spec}{plot.spec}}.
\end{ldescription}
\end{Arguments}
\begin{Value}
An object of class \code{"spec"}.
The result is returned invisibly if \code{plot} is true.
\end{Value}
\begin{Section}{Warning}
Some authors, for example Thomson (1990), warn strongly
that AR spectra can be misleading.
\end{Section}
\begin{Note}\relax
The multivariate case is not yet implemented.
\end{Note}
\begin{References}\relax
Thompson, D.J. (1990) Time series analysis of Holocene climate data.
\emph{Phil. Trans. Roy. Soc. A} \bold{330}, 601--616.

Venables, W.N. and Ripley, B.D. (2002) \emph{Modern Applied
Statistics with S.} Fourth edition. Springer. (Especially
page 402.)
\end{References}
\begin{SeeAlso}\relax
\code{\LinkA{ar}{ar}}, \code{\LinkA{spectrum}{spectrum}}.
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
spec.ar(lh)

spec.ar(ldeaths)
spec.ar(ldeaths, method="burg")
\end{ExampleCode}
\end{Examples}

