[R] ARIMA0 with xreg,
Brian Scholl
brianscholl1973 at yahoo.com
Thu Dec 20 21:08:59 CET 2001
Fair enough. I won't bother to define things since I
think what you are saying is that if X is my Txr
matrix of time series data, make xreg=X[,2:r] for the
regression on X[,1], which means I misunderstood the
function specifications. Sorry, and thanks.
--- Prof Brian Ripley <ripley at stats.ox.ac.uk> wrote:
> You need to define some terms. But from `VAR' and
> `matrix of data' I
> guess that you have a multiple time series. If so,
> it just isn't valid
> to use arima0, which is univariate. You also did
> not say what `x' is,
> but the `x' in `xreg' means eXogenous.
>
> On Wed, 19 Dec 2001, Brian Scholl wrote:
>
> > Hi all,
> >
> > Using ar(), I fit a VAR to my time series that has
> a
> > reasonably 'nice' error spectrum and aic
> determines
> > p=7. But the output for ar isn't quite as
> convenient
> > as arima0, namely in that it takes me an extra
> step to
> > get the s.e.'s of parameters and it doesn't
> produce an
> > estimate of the log-likelihood for comparison to
> other
> > models. So I thought I'd use arimia0, with xreg=x,
> my
> > matrix of data and order=(7,0,0).
> >
> > When I do this, the results are funny. The
> residuals
> > look - well, just like the original series on a
> > readjusted scale. The spectrum for the residuals
> > looks slightly different from the spectrum of x
> and
> > bears no resemblance to the nice spectrum I got
> before
> > (this one looks like a long memory process, while
> the
> > earlier was closer to white noise). Certainly I
> don't
> > expect the spectra and residuals to look exactly
> alike
> > because of the different solution methods, but
> these
> > bear little resemblance to each other.
> >
> > I assume, of course that the error is on my part,
> > perhaps something I've overlooked, but I'm not
> able to
> > find it.
> >
> > Thanks,
> >
> > Brian
> >
> >
> >
> > __________________________________________________
> >
> >
> >
> >
> >
>
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> > r-help 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:
> r-help-request at stat.math.ethz.ch
> >
>
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> >
>
> --
> Brian D. Ripley,
> ripley at stats.ox.ac.uk
> 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-help 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:
> r-help-request at stat.math.ethz.ch
>
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
__________________________________________________
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list