R Удалить несколько сезонностей в одном временном ряду

Question

Я имею дело с почасовыми данными, и мои данные имеют сезонность 24 часа и 168 часов.

Когда я использую этот код. Я удалил 24-часовую сезонность.

deseasonalized_ec_hourly<-electricity_consumption_ts_hourly/ec_hourly_multip$seasonal

Но 168-часовая сезонность осталась. Как убрать 168 часов сезонности

Answer 1

Вот пример получасового потребления электроэнергии из Англии и Уолса (июнь - август 2000 г.). Сезонность (ежедневная и еженедельная) оценивается с помощью функции mstl() из пакета forecast. Функция seasadj() удаляет сезонные компоненты.

library(forecast)

autoplot(taylor)

adjust <- taylor %>% mstl() %>% seasadj()
autoplot(adjust)

^{Создано 07.05.2020 пакетом REPEX (v0.3.0)}