在 hts 包中，为什么 gts 在创建组结构时会硬编码 2 组合？

Question

在hts包中，分组时间序列结构是使用gts创建的。该函数调用另一个名为 CreateGmat 的函数，其中一个 2 组合

cn <- combn(1L:total.len, 2)

被计算并用于指导构建由"g-matrix"表示的组结构。对于类别较少的较小结构，仅使用 2 组合来创建结构具有直观意义，但我想知道是否应该始终如此。

我还没有完成所有代码，但是我的直觉是 G 矩阵决定了 S 矩阵，它影响调节矩阵，调节矩阵影响预测使用最佳组合方法计算。

除非我对协调过程如何发生的理解有误，否则 G 矩阵的公式化方式至少对预测的计算方式有些重要。对于具有许多类别的大型结构，包括大于两个的组合是否会以任何显着积极的方式影响预测？仅使用 2 组合是否有统计或实验原因？

Answer 1

你对G矩阵的理解是完全正确的。如果有超过 3 个分组变量，我们需要创建一个具有 2 向和 3 向交互作用的 G 矩阵。去年我有一些代码可以为一个项目自动执行超过 3 种方式的组合，但最终没有把它放到 hts 包中。

nvars <- length(chr_key)
possible_len <- seq_len(nvars - 1)
list_comb <- lapply(possible_len, function(x) combn(seq_len(nvars), x,
  simplify = FALSE))
cs_comb <- c(0, cumsum(choose(nvars, possible_len)))
len_comb <- cs_comb[nvars]

在 hts 包中，为什么 gts 在创建组结构时会硬编码 2 组合？

In the hts package, why does the gts function hard code 2-combinations when creating a group structure?

r

time-series

forecasting