Skip to content

tseries(arma) vs. stats(arima)

2 messages · Richard Saba, Brian Ripley

#
Hello,
The "arma"  function in the "tseries"  package allows estimation of models
with specific "ar" and  "ma" lags  with its "lag" argument. 
For example:  y[t] = a[0] + a[1]y[t-3] +b[1]e[t-2] + e[t]  can be estimated
with the following specification :   arma(y, lag=list(ar=3,ma=2)).

Is this possible with the "arima" function in the "stats" or in other time
series packages like fArima, forecast, or FinTS? They all take a "lag"
argument. I  would like to have the ability to estimate models like the one
above while utilizing the "xreg" argument available in the other arima
functions .
Thanks,
Richard Saba
sabaric at auburn.edu
#
arima() certainly does.  You missed the 'fixed' argument that allows ARMA 
coefficients to be set to any value, including 0.

Note that arma() does not do ML estimation.  For that you have to worry 
about invertibility, and it can be much harder to do the optimization with 
constrained parameters.
On Fri, 21 Mar 2008, Richard Saba wrote: