*
* Example SERIESC.DAT from pp 388-389
*
open data seriesc.dat
data(format=free,org=columns) 1 48 c
*
graph(footer="SERIESC Microcomputer Demand")
# c
*
boxjenk(diffs=1,const,ar=1,ma=1,define=ceq) c
*
uforecast(equation=ceq,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
*
graph(footer="Figure 9-8. 95% prediction interval for SERIESC.DAT") 4
# c
# forecast
# lower / 3
# upper / 3