使用带有R hts包的Croston模型

Question

我认为这是与this post passing forecasting method to hts

类似的帖子

给出的示例是MAPE包。是否可以以类似的方式在预测包中使用croston函数？

我只是试过

{{1}}

但是这给了我一个错误，说要替换的项目数量不是替换长度的倍数。

Answer 1

查看croston返回的对象的结构。它不是一个简单的向量。以下代码是一个有效的示例。

library(hts)
nodes <- list(2, c(3, 2))
abc <- ts(5 + matrix(sort(rnorm(500)), ncol = 5, nrow = 100))
bts <- hts(abc, nodes)
all_ts <- aggts(bts)
allf <- matrix(NA, nrow = 3, ncol = ncol(all_ts))
for(i in 1:ncol(all_ts)){
  allf[,i] <- croston(all_ts[,i],h = 3)$mean
}
y.f <- combinef(allf, bts$nodes)