[R] Time series misalignment

Fernando Saldanha fsaldan1 at gmail.com
Wed Apr 13 00:47:21 CEST 2005


Can one also predetermine a set and then estimate all the models one
wants to compare using the zoo package? Or can that be done only with
the tseries package?

Thanks.

FS 

On 4/12/05, Achim Zeileis <Achim.Zeileis at wu-wien.ac.at> wrote:
> Fernando:
> 
> > This maybe a basic question, but I have spent several hours
> > researching and I could not get an answer, so please bear with me. The
> > problem is with time series in the package tseries.
> 
> BTW: the `tseries' package is not involved here.
> 
> > As the example
> > below shows, the time series can get misaligned, so that bad results
> > are obtained when doing regressions.
> 
> lm() per se has only very limited support for time series regression.
> Therefore, there are currently several tools under development for
> addressing this issue. In particular, Gabor Grothendieck and myself are
> working on different approaches to this problem.
> 
> <snip>
> 
> > To fix this problem I did the following:
> >
> > > tsf <- ts.intersect(y1, x1, z1)
> >
> > Now I can do:
> >
> > > lm1 <- lm(tsf[,3] ~ tsf[,2])
> 
> It is probably simpler to just do
>   lm1 <- lm(z1 ~ x1, data = tsf)
> 
> Another approach is implemented in the zoo package. This implements an
> formula dispatch and you can do
>   lm1 <- lm(I(z1 ~ x1))
> *without* computing tsf first.
> 
> Depending on what you want to do with the fitted models, one of the two
> approaches might be easier, currently. In particular, if you want to fit
> and compare several models, then I would compute the intersection first
> and fit all models of interest on this data set.
> 
> Furthermore, note that the dispatch implementation via I() in zoo is
> still under development and likely to change in future versions. (But
> this mainly means that improved implementations will become available
> soon, stay tuned :-)
> Z
> 
>




More information about the R-help mailing list