Início / coding / Perguntas / 79566422
Accepted
Saïd Maanan
Asked: 2025-04-10 18:50:46 +0800 CST

Como criar um gráfico de tendência distinto com o AutoSTR?

  • 772

Tenho a seguinte carga elétrica diária para a Espanha:

> head(spain_daily)
# A tibble: 6 × 2
  Date       Daily_Load
  <date>          <dbl>
1 2022-01-01     20872.
2 2022-01-02     22729.
3 2022-01-03     27358.
4 2022-01-04     28262.
5 2022-01-05     28048.
6 2022-01-06     24182.
> tail(spain_daily)
# A tibble: 6 × 2
  Date       Daily_Load
  <date>          <dbl>
1 2024-12-26     25210.
2 2024-12-27     26363.
3 2024-12-28     25041.
4 2024-12-29     24182.
5 2024-12-30     27073.
6 2024-12-31     25924.

Como o consumo elétrico é cíclico, decompus a série usando os seguintes comandos:

library(forecast)
library(stR)

autoplot(mstl(msts(spain_daily$Daily_Load, seasonal.periods = c(7, 365.25))))

O que me deu o seguinte enredo:

decomposição_de_séries

Tentei recriar o gráfico acima usando a função stRdo pacote AutoSTRcom o seguinte comando:

plot(AutoSTR(msts(spain_daily$Daily_Load, seasonal.periods = c(7, 365.25))))

E me dá o seguinte gráfico:

segunda_decomposição

Como tornar este gráfico mais semelhante ao primeiro, ou seja, exibir a série de tendências em seu próprio gráfico e remover a série de previsão (azul)?

Obrigado.

1 respostas

  1. Best Answer

    O primeiro gráfico usa ggplot2, enquanto o segundo usa gráficos base. Portanto, você precisa extrair os componentes do objeto produzido por AutoSTRe usar ggplot2 para plotá-los. Felizmente, existem funções que facilitam isso, como mostrado no exemplo a seguir.

    ``` r
library(stR)
library(forecast)
library(ggplot2)

# synthetic data
spain_daily <- msts(
  sin(2*pi*seq(3*365)/365) + sin(2*pi*seq(3*365)/7) + rnorm(3*365, 0, 0.5),
  seasonal.periods = c(365, 7)
)
fit <- AutoSTR(spain_daily)
components(fit) |>
  autoplot(colour = FALSE) +
  facet_grid(series ~ ., scales = "free_y")

    Criado em 2025-04-11 com reprex v2.1.1

    • 1

relate perguntas