[R] partial analisys of a time series
PIKAL Petr
petr.pikal at precheza.cz
Tue Dec 4 11:07:48 CET 2012
Hi
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of Antonio Silva
> Sent: Tuesday, December 04, 2012 10:26 AM
> To: R-help at r-project.org
> Subject: [R] partial analisys of a time series
>
> Dear list members
>
> I want to analyze separately the months of a time series. In other
> words, I want to plot and fit models for each month separately.
>
> Taking the example of
> http://a-little-book-of-r-for-time-
> series.readthedocs.org/en/latest/src/timeseries.html
>
> births <- scan("http://robjhyndman.com/tsdldata/data/nybirths.dat")
> birthstimeseries <- ts(births, frequency=12, start=c(1946,1))
> birthstimeseries
> plot.ts(birthstimeseries)
> birthstimeseriesHW <- HoltWinters(birthstimeseries)
> plot(birthstimeseriesHW)
>
> How to proceed the plotting and HoltWinters smoothing considereing only
> Januarys, Februarys, etc. separately.
Split your data by months to a list, use lapply.
using zoo package
blist <-split(birthstimeseries, months(as.Date(birthstimeseries)))
l.blist <- lapply(blist, HoltWinters)
Regards
Petr
>
> Thanks in advance.
>
> Antonio Olinto
>
> [[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