library(tidyverse)
library(fpp3)
library(plotly)
library(readxl)
library(tsibble)
library(tibble)
library(zoo)
library(lubridate)Prueba
Datos
datos <- read_excel("V.xlsx", sheet = 2)datos$Date <- as.Date(datos$Date)datos <- as_tsibble(datos, index = Date)
datosEstimación de modelos
datos_fit <- datos |>
model(
media = MEAN(Sales),
naive = NAIVE(Sales),
snaive = SNAIVE(Sales),
drift = RW(Sales ~ drift())
)
datos_fitAnálisis de residuos
datos_fit |>
augment()datos_fit |>
augment() |>
features(.innov, ljung_box, lag = 10, dof = 0)Media
datos_fit |>
select(media) |>
gg_tsresiduals()Naive
datos_fit |>
select(naive) |>
gg_tsresiduals()Seasonal Naive
datos_fit |>
select(snaive) |>
gg_tsresiduals()Drift
datos_fit |>
select(drift) |>
gg_tsresiduals()Forecast
datos_fc <- datos_fit |>
forecast(h = "8 days")
datos_fcdatos_fc |>
autoplot(datos) +
facet_wrap(~.model)datos_fc |>
autoplot(datos) +
facet_wrap(~.model) +
theme(legend.position = "bottom")