在R中使用dshw()的非嵌套双季节性

时间:2013-05-20 19:57:07

标签: r forecasting decomposition stl-decomposition

我正在尝试使用dshw()来处理双重季节性 - 在我的情况下,每日数据包含一周(7天)和一年(365天)季节性。但是,当我运行代码时出现以下错误:

data<-msts(1:1000, seasonal.periods=c(7,365), ts.frequency=365, start=2012)
decompose<-dshw(data, period1=7, period2=365)
 -- Error in dshw(data, period1 = 7, period2 = 365) : Seasonal periods are not nested

您认为解决此问题的最佳做法是什么?我应该只对我的数据使用stl两次(7天和365天的频率)吗?或者以某种方式修改数据?

谢谢!

2 个答案:

答案 0 :(得分:2)

请尝试使用tbats()模型。它专门设计用于避免此问题。 DSHW是TBATS模型的特例。

decompose <- tbats(data)

答案 1 :(得分:0)

如果您根据period1定义了period2,则不会收到错误消息。

代替:

分解<-dshw(数据,period1 = 7,period2 = 365)

使用:

分解<-dshw(数据,period1 = 7,period2 = 7 * 52)