EmpezarEmpieza gratis

Regresión armónica para estacionalidad múltiple

Las regresiones armónicas también son útiles cuando las series temporales tienen varios patrones estacionales. Por ejemplo, taylor contiene la demanda eléctrica semihoraria en Inglaterra y Gales durante unos meses del año 2000. Los periodos estacionales son 48 (estacionalidad diaria) y 7 x 48 = 336 (estacionalidad semanal). No hay suficientes datos como para considerar la estacionalidad anual.

auto.arima() tardaría mucho en ajustar una serie larga como esta, así que en su lugar ajustarás un modelo de regresión estándar con términos de Fourier usando la función tslm(). Es muy similar a lm(), pero está diseñada para trabajar con series temporales. Con estacionalidad múltiple, tienes que especificar el orden \(K\) para cada uno de los periodos estacionales.

# El argumento formula es una descripción simbólica
# del modelo que se va a ajustar

> args(tslm)
function (formula, ...)

tslm() es una función recién introducida, así que deberías poder seguir el código preescrito en su mayor parte. Los datos taylor están cargados en tu espacio de trabajo.

Este ejercicio forma parte del curso

Forecasting en R

Ver curso

Instrucciones del ejercicio

  • Ajusta una regresión armónica llamada fit a taylor usando orden 10 para cada tipo de estacionalidad.
  • Pronostica 20 días laborables por delante como fc. Recuerda que los datos son semihorarios para establecer el valor correcto de h.
  • Crea un gráfico temporal de los pronósticos.
  • Revisa los residuales de tu modelo ajustado. Como verás, auto.arima() habría hecho un mejor trabajo.

ejercicio interactivo práctico

Prueba este ejercicio completando este código de ejemplo.

# Fit a harmonic regression using order 10 for each type of seasonality
fit <- tslm(taylor ~ fourier(___, K = c(10, 10)))

# Forecast 20 working days ahead
fc <- forecast(___, newdata = data.frame(fourier(___, K = ___, h = ___)))

# Plot the forecasts
___

# Check the residuals of fit
___
Editar y ejecutar código