我有一个聚合的分层数据集(每小时聚合)如下
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)
有什么建议么 ?
谢谢