The following is an internal function in arima()
arCheck <- function(ar) {
p <- max(which(c(1, -ar) != 0)) - 1
if (!p)
return(TRUE)
all(Mod(polyroot(c(1, -ar[1:p]))) > 1)
}
On Tue, 10 Apr 2007, Leeds, Mark (IED) wrote:
> I've looked around but I can't find the method in R for testing whether
> the resulting estimated coefficients
> of an AR model imply that the model is invertible.
>
> To quote from eric zivot's blue book :
Which really doesn't help us. (I suspect you mean Zivot and Wang.)
This is standard times series material.
>
> " the AR(p) is invertible provided the rots of the characteristic
> equation
>
> Phi(z) = 1 - phi_1*z - phi_2*z^2 = phi_3*z^3 - ..... Phi_p*z^p = 0 lie
> outside
> the complex circle".
>
> I can't find a function nor do I know how to do the above myself. I
> think there is an equivalent method in which
> I can check whether the eigenvalues of some dual equation ( I forget
> what it is ) are less than one but I don't
> remember exactly what that equation is and , even if I did, I still
> wouldn't know how to do it.
>
> Maybe checking whether the absolute of the sum of the coefficients is
> less than one is okay ?
> I remember doing that in another life but I'm not sure if that's an
> approximation or an actual test.
>
>
> Thanks for any help.
> --------------------------------------------------------
>
> This is not an offer (or solicitation of an offer) to buy/se...{{dropped}}
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
