分层预测(hts包)R不均匀组和自定义预测

时间:2016-03-24 20:02:51

标签: r forecasting

我正在阅读我的hts包。我有兴趣尝试自定义预测部分,该部分需要使用指定的预测函数循环遍历每个系列。提供的示例显示了偶数组。我想知道如何运行具有不均匀组数的自定义预测,然后将它们重新组合成gts对象。

bts <- ts(5 + matrix(sort(rnorm(500)), ncol=5, nrow=100))
y <- hts(bts, nodes=list(2, c(3, 2)))

allts_y <- aggts(y)
allf <- matrix(allts_y, nrow=100, ncol=ncol(allts_y))

for(i in 1:ncol(allts_y))
  allf[,i] <- forecast(auto.arima(allts_y[,i]), h=10)$mean
allf <- ts(allf, start=1)

# below code is where I run into a hang-up
g <- matrix(c(rep(2, 5), rep(3, 5), rep(1:5, 2)), nrow = 2, byrow = T)
y.f <- combinef(allf, groups = g)

1 个答案:

答案 0 :(得分:1)

您的组结构gnodeshts()参数指定的层次结构不匹配。用

替换最后一行
y.f <- combinef(allf, nodes=y$nodes)