HTS包:如何指定类似网络的预测层次?

age*_*nis 6 r graph time-series hierarchical forecasting

我正在尝试使用这个hts 包进行分层预测,以便一起预测水网的不同尺度,以便进行协调预测.水网络由流量计组成,流量计以正或负的方式测量流量,当我们对净流量进行聚集时,它给出了一个邻域的内部消耗.

这不是100%的等级敏感因为: - 流量计对于进入一个邻域的流量是积极的,但是在离开另一个邻域时是负的.

所以结构是这样的,简化了3个流量计和2个邻域:

在此输入图像描述

所以我试图复制F2仪表,一个正面和一个负面,但输出效果不理想.我怎么处理这个?

m=structure(c(28, 20, 20, 17, 19, 22, 26, 38, 45, 45, 42, 38, 37, 
            37, 33, 29, 28, 31, 37, 40, 36, 31, 26, 23, 22, 15, 13, 11, 11, 
            12, 16, 22, 25, 24, 21, 19, 17, 16, 14, 13, 13, 14, 16, 17, 16, 
            13, 11, 8, 7, 7, 8, 8, 8, 8, 8, 9, 11, 11, 11, 12, 12, 12, 11, 
            11, 11, 10, 10, 9, 9, 9, 8, 5), .Dim = c(24L, 3L), .Dimnames = list(
              NULL, NULL))
m=cbind(m, -m[, 3])
m[,2]<- -m[,2]
colnames(m) <- c("NaF1", "NbF3", "NbF2", "NaF2")
hm <- hts(m, characters=c(2, 2))
plot(hm); forecast(hm, fmethod="arima")$bts
#  
#        NaF1      NaF2     NbF2       NbF3
# 25 23.83012 -5.019900 4.799014  -8.505145
# 26 26.24205 -5.257814 4.989756 -10.778844
# 
Run Code Online (Sandbox Code Playgroud)

正如您所看到的,F2的预测在正负方向上是不同的,这是没有意义的.

为了澄清聚合的结构,这里是一个grafcet: 在此输入图像描述

小智 0

难道您不能将 Na 的 F.2 负值视为 Na 的子节点,将 Nb 的 F.2 正值视为 Nb 的子节点吗?

这会将其转变为更简单的层次结构。