Я провожу сезонный / трендовый анализ временных рядов в R. Я использую TSLM и forecast в пакете fable. Я написал регрессию для своей модели, чтобы фиксировать сезонность и тенденцию за месяц, но теперь я пытаюсь спрогнозировать 3 года, а затем построить фактические, подогнанные и прогнозируемые результаты. Когда я включаю функцию forecast, я получаю эту ошибку Error: Can't slice a scalar. Я не могу понять, почему и как решить эту проблему. Приветствуются любые ссылки / предложения.
df <- matrix(exp(rnorm(47 * 64)), ncol = 47)
colnames(df) <- c(letters, LETTERS)[1:47]
df <- as.data.frame(df)
df <- df %>%
mutate(
Month = yearmonth(seq(as.Date("2016-01-01"), by = "1 month", length = NROW(df)))
) %>%
pivot_longer(-Month, names_to = "Series", values_to = "value") %>%
as_tsibble(index = Month, key = Series)
reg<-df %>% model(
TSLM(log(value) ~ trend() + season())) %>%forecast(h = "3 years")