Re: [R] Time series reliability questions

From: Achim Zeileis <>
Date: Thu, 24 Jul 2008 01:15:26 +0200 (CEST)

On Wed, 23 Jul 2008, John Theal wrote:

> EViews and Gretl give comparable (and I am inclined to presume, correct)
> results. R on the other hand, has the exogenous regressor with a negative
> coefficient. If I use other data I encounter the same problem - agreement
> between EViews and Gretl, disagreement with R (for identical data sets). Are
> there any known bugs with arima estimation in R?

You can check at

> If I use the Zelig package, I
> get the same results as the arima{stats} function call.

...because zelig() is just an interface to arima() for ARIMA models.

> If I remove the
> exogenous regressor from the estimations then I have agreement between R, Gretl
> and EViews, but with the exogenous regressor (basically an ARMAX model) the
> estimation results are substantially different.

Log-likelihood, innovation variance and AIC seem to be similar though. It might be a difference in parametrization between R and gretl/EViews. Z mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Wed 23 Jul 2008 - 23:30:35 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 23 Jul 2008 - 23:32:17 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive