Subject: Re: arima0() (PR#314)
From: Prof Brian D Ripley (firstname.lastname@example.org)
Date: Tue 09 Nov 1999 - 18:11:17 EST
On Sun, 7 Nov 1999 email@example.com wrote:
> Full_Name: Ahmad Abu Hammour
> Version: rw0651
> OS: windows 95
> Submission from: (NULL) (188.8.131.52)
> Although I know that "ts package" is preliminary, I wanted to compare the
> results from R and SPSS. I ran ARIMA(2,1,2) in both softwares. I got NaN in
> standard errors of coefficients from R and real figures from SPSS. I changed
> "delta" in R to match that used by SPSS, I received results with no NaNs but
> different from those from SPSS specially in MA coefficients. They came out
> negative numbers while positive from SPSS.
Mr Hammour sent me the data privately. Running
arima0(y0, order=c(2,1,2), delta=-1)
arima0(x = y0, order = c(2, 1, 2), delta = -1)
ar1 ar2 ma1 ma2
-0.7686 0.1519 -0.1143 -0.8857
Approx standard errors:
ar1 ar2 ma1 ma2
0.1917 0.0912 0.1765 0.1762
sigma^2 estimated as 109859: log likelihood = -890.65, aic = 1789.3
[Using the default delta of 0.01 in this case gives inaccurate results:
but it is a very badly determined model.]
Number of residuals 123
Standard error 337.30778
Log likelihood -890.6952
Variables in the Model:
B SEB T-RATIO APPROX. PROB.
AR1 -.78319350 .16013797 -4.8907421 .00000318
AR2 .16191191 .09435030 1.7160719 .08875107
MA1 .07993592 .54769107 .1459508 .88420721
MA2 .91762463 .49337224 1.8599032 .06536851
>Warning # 16567. Command name: ARIMA
>Our tests have determined that the estimated model lies close to the
>boundary of the invertibility region. Although the moving average
>parameters are probably correctly estimated, their standard errors and
>covariances should be considered suspect.
So (up to the MA sign change, a matter of definition) the estimates
differ by 1 se or so. However, SPSS has given a warning, and the
log-likelihoods are similar but R's is larger. This suggests that
- There is no practical difference between the results, and
- R's optimizer has done a better job.
- R works on transformed scale and so has more chance of getting se's
close to the boundary (but the error distribution will be very far from
- This is not a good model for the data.
If there is a bug at all here, it would appear to be in SPSS.
Moral: take warnings seriously.
-- Brian D. Ripley, firstname.lastname@example.org 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 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: email@example.com _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
This archive was generated by hypermail 2b25 : Tue 04 Jan 2000 - 14:16:09 EST