*
* SERIESB.DAT example from pp 383-384
*
open data seriesb.dat
data(format=free,org=columns) 1 48 seriesb
*
boxjenk(diffs=1,ar=1,define=arima110) seriesb
*
uforecast(equation=arima110,stderrs=efse) forecast 49 84
set lower 49 84 = forecast+%invnormal(.025)*efse
set upper 49 84 = forecast+%invnormal(.975)*efse
print(picture="*.####") 49 84 forecast efse lower upper