多级/层次时间序列

Question

我有一个数据框（df），其日期为过去两年

数据：

Country State County  2015-01-01   2015-02-01 2015-03-01
IN       A     a       400          500       300    
IN       A     b       231          535       214
IN       B     a       532          212       242
........So on  .......................

目的：

我想为每个组合建立时间序列模型，并获得未来3个月的预测

预期结果：

Country State County  2015-01-01 2015-02-01 2015-03-01  .... 2017-01-01  
IN       A     a       400          500       300            421 
IN       A     b       231          535       214            231
IN       B     a       532          212       242            313
   ......... so on for all combination......

问题：

如何在每个级别/组合（在国家/地区，国家和州级别）构建模型？

我厌倦了以下事情：

library(broom) 

time_series_models <- df %>% rowwise() %>% 
  mutate(forecasted_values = forecast(ets(ts(df[3:NCOL(df)], 
                                             start = c(2015, 5), frequency = 12), h= 4)))

问题：

错误：y应该是单变量时间序列