Fiz uma previsão usando o auto.arima
onde minha base de dados são valores mensais de jan/2018 a set/2019.
Minha base de treino é de jan/2018 a jun/2019:
VL_TR_treino_5S = window(VL_TR_TS_5S, start=c(2018,1), end=c(2019,6))
VL_TR_teste_5S = window(VL_TR_TS_5S, start=c(2019,6))
E para aplicar o auto.arima usei outras covariáveis em xreg
, como exemplo coloquei A, B e C:
VL_TR_modelo_5S = auto.arima(VL_TR_treino_5S, xreg = cbind(A,B,C), trace = T, stepwise = T, approximation = T, seasonal = T)
E depois usei o forecast usando um período de 24 meses:
VL_TR_Prev5S = forecast(VL_TR_modelo_5S, xreg = cbind(A,B,C), h = 24)
Porém quando vou visualizar os dados de VL_TR_Prev5S
, ao invés de me mostrar 24 valores previstos (que seria até dez/2020), me mostra apenas 13 valores que seria de fev/2019 até fev/2020.
print(VL_TR_Prev5S)
Point Forecast Lo 20 Hi 20
Feb 2019 7649351 7634063 7664639
Mar 2019 8260246 8244958 8275534
Apr 2019 8950091 8934803 8965380
May 2019 8657965 8642677 8673253
Jun 2019 8534740 8519451 8550028
Jul 2019 8349148 8333859 8364436
Aug 2019 7596208 7580920 7611496
Sep 2019 8515507 8500218 8530795
Oct 2019 8103160 8087871 8118448
Nov 2019 8143330 8128042 8158619
Dec 2019 7393488 7378199 7408776
Jan 2020 7007616 6992328 7022905
Feb 2020 6819635 6804346 6834923
Quando eu executo o script do auto.arima
, apesar do algoritmo rodar normalmente, o r
me deixa o seguinte aviso:
Warning message:
The chosen seasonal unit root test encountered an error when testing for the first difference.
From stl(): series is not periodic or has less than two periods
0 seasonal differences will be used. Consider using a different unit root test.
Não sei se esse aviso influência na questão, mas optei em mencionar o aviso por desencargo. Pesquisando em alguns fóruns parece que usando covariáveis em xreg
pode limitar os períodos, mas não sei o porque e também não sei como posso evitar isso.
Mas enfim, como posso usar o auto.arima
para prever 24 períodos ou mais?