0

我有一个聚合的分层数据集(每小时聚合)如下

  TIME                level1 level2     n
  <dttm>              <fct>  <fct>  <int>
1 2019-11-29 00:00:00 A      a          1
2 2019-11-29 00:00:00 B      a          1
3 2019-11-29 00:00:00 B      b          2
4 2019-11-29 00:00:00 B      c          1
5 2019-11-29 01:00:00 A      a          2
6 2019-11-29 02:00:00 A      a          1
7 2019-11-29 02:00:00 B      a          1
8 2019-11-29 02:00:00 B      b          2
9 2019-11-29 02:00:00 B      c          1

我想将它转换为 hts 对象,并按照 hts 手册中的说明和交叉验证的这篇文章: https ://stats.stackexchange.com/questions/133894/get-groups-in-time-series-with-categorical -data-in-r-for-use-in-gts

我已将数据转换为时间序列(mts)对象

df_mat <- ts(data.matrix(df_agg), frequency = 525960)

频率设置为与数据框中的聚合单位一致的小时。

由于不同的 level1 和 level2 的数量可以变化,我没有设置 nrow 和 ncol 对象。

在此之后,我不太确定如何将此对象转换为 hts。如何生成要传递给 hts 方法的参数节点、列名和字符?

hts(y, nodes, bnames = colnames(y), characters)

有什么建议么 ?

谢谢

4

0 回答 0