I have a electricity consumption daily time series. Data are from the last 3 years.
I have adjusted a (0,1,1) model after difference (1) the log transformed series. (checked stationarity). I have identified some outliers but any of them resulted significant when I included it as regressors.
It seems that there is not autocorrelation in the residuals (tests OK) but the heteroscedasticity tests are not fulfilled. Tried to model the residuals with arch models but cannot fix the problem of heteroscedasticity (conditional and unconditional).
Any alternatives?
proc arima data=work.series;
identify var=lvalor(1) stationarity=(dickey) nlag=30;
estimate p=(1) q=(1) noint method=ml;
*outlier maxnum=5 alpha=0.01 id=fecha;
forecast id=fecha interval=day out=b;
run;
quit;
proc model data=b;
parms const;
residual = const ;
fit residual / white breusch=(1 time);
run;
quit;
proc autoreg data=b;
model residual= / garch=(Q=2, p=2) maxit=50;
output out=out R=pepe;
run;