[R] GLS models - bootstrapping

From: Lillian Sandeman <l.sandeman_at_abdn.ac.uk>
Date: Mon 02 Oct 2006 - 12:59:09 GMT


Hello,

I am have fitted GLS models to time series data. Now I wish to bootstrap this data to produce confidence intervals for the model.

However, because this is time series data, normal bootstrapping is not applicable. Secondly, 'tsboot' appears to only be useful for ar models - and does not seem to be applicable to GLS models.

I have written code in R to randomly sample blocks of the data (as in Davison & Hinkley's book - bootstrap methods and their application) and use this resampling to re-run the model, but this does not seem to be the correct approach since Confidence Intervals produced do not show the underlying pattern (cycles) in the data [even when block length is increased, it only picks up a little of this variation].

Any help as to how to proceed with this would be greatly appreciated, as I cannot find anything applicable on the R pages. Alternatively, if there is another method to proceed with this (other than bootstrapping), I would also be happy to try it.

Thankyou,

Lillian.



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Mon Oct 02 23:13:35 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Mon 02 Oct 2006 - 13:30:07 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.