*
* Example MA1P.DAT from page 385
*
open data ma1p.dat
data(format=free,org=columns) 1 100 time ma1p
*
boxjenk(const,ma=1,define=ma1eq) ma1p
*
uforecast(equation=ma1eq,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