[R] Problem with zoo and BootPR packages
Achim Zeileis
Achim.Zeileis at wu-wien.ac.at
Thu Nov 19 18:06:54 CET 2009
On Thu, 19 Nov 2009, Ricardo Gonçalves Silva wrote:
> Hi,
>
> I'm trying to plot the forecasts I generated using the Plot.Fore function of the BootPR package.
> But I got an error from zoo:
>
> My data:
>
> Time Series:
> Start = 1
> End = 18
> Frequency = 1
> [1] 38731 38628 39117 92809 71984 31226 58613 72360 107956 92066
> [11] 95208 99098 95848 120383 110717 105680 98469 101916
>
> Script:
>
> y1<-ts(y1);
> forey1<-BootBC(y1,p=2,h=3,nboot=5000,type="const+trend",correct="ssf")
> Plot.Fore(y1,forey1$forecast,start=1966,end=1984,frequency=1)
>
> The Error:
>
> Error in zooreg(x, start, end, frequency) :
> "data" : attempt to define illegal zoo object
>
> Any Help?
This is a bug in Plot.Fore() which does not use the zoo functions
correctly. It should be reported to the package maintainer.
To avoid it, you can do
Plot.Fore(as.vector(y1), ...)
instead of
Plot.Fore(y1, ...)
Note to the maintainer (Jae Kim, Cc now): Plot.Fore() calls zooreg(x)
where x is a "ts" object. This isn't the appropriate use of zooreg() which
expects a numeric vector/matrix (or a factor). In the case above
as.zoo(y1) would already be enough (and preserve all time information). Or
you can manually call zooreg(coredata(y1), start = ..., end = ..., ....).
See the zoo vignettes/examples for more details.
Best,
Z
> Thanks
>
> Rick
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org 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.
>
>
More information about the R-help
mailing list