*
* Example AR1P.DAT from pp 381-382
*
open data ar1p.dat
data(format=free,org=columns) 1 100 time ar1p
*
boxjenk(ar=1,const,define=ar1eq) ar1p
*
uforecast(equation=ar1eq,stderrs=efse) forecast 101 124
set lower 101 124 = forecast+%invnormal(.025)*efse
set upper 101 124 = forecast+%invnormal(.975)*efse
*
print(picture="*.###") 101 124 forecast efse lower upper