У меня есть агрегированный иерархический набор данных (агрегированный за час), как показано ниже
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
Я преобразовал данные в объект временного ряда (мтс) как
df_mat <- ts(data.matrix(df_agg), frequency = 525960)
Частота, установленная в час, с единицей агрегации во фрейме данных.
Поскольку число различных уровней level1 и level2 может варьироваться, я не установил nrow и ncol объект.
После этого я не очень уверен, как преобразовать этот объект в hts. Как сгенерировать аргументы узлов, имен и символов для передачи в метод hts?
hts(y, nodes, bnames = colnames(y), characters)
Есть предложения?
Спасибо