Prevendo agendamentos de chamadas
Outra série temporal com múltiplos períodos sazonais é calls, que contém 20 dias consecutivos de dados de volume de chamadas a cada 5 minutos de um grande banco norte-americano. Existem 169 períodos de 5 minutos em um dia útil e, portanto, a frequência sazonal semanal é 5 x 169 = 845. A sazonalidade semanal é relativamente fraca, então aqui você vai modelar apenas a sazonalidade diária. calls já está carregado no seu ambiente de trabalho.
Os resíduos, neste caso, ainda não passam nos testes de ruído branco, mas suas autocorrelações são minúsculas, embora significativas. Isso acontece porque a série é muito longa. É comum que resíduos não passem nos testes em séries tão extensas. O efeito das correlações restantes nas previsões será desprezível.
Este exercício faz parte do curso
Previsão em R
Instruções do exercício
- Faça o gráfico dos dados
callspara ver a forte sazonalidade diária e a fraca sazonalidade semanal. - Monte a matriz
xregusando ordem10para a sazonalidade diária e0para a sazonalidade semanal. Observe que, se você especificar o vetor incorretamente, sua sessão pode expirar! - Ajuste um modelo de regressão dinâmica chamado
fitusandoauto.arima()comseasonal = FALSEestationary = TRUE. - Verifique os resíduos do modelo ajustado.
- Crie as previsões para 10 dias úteis à frente como
fce, em seguida, faça o gráfico. A descrição do exercício deve ajudar você a determinar o valor adequado deh.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Plot the calls data
___
# Set up the xreg matrix
xreg <- fourier(___, K = ___)
# Fit a dynamic regression model
fit <- auto.arima(___, xreg = ___, ___, ___)
# Check the residuals
___
# Plot forecasts for 10 working days ahead
fc <- forecast(fit, xreg = fourier(calls, c(10, 0), h = ___))
___