*
* PHARMDEM.DAT example from pp 387-388
*
open data pharmdem.dat
cal(m) 1980
data(format=free,org=columns) 1980:1 1992:12 pharmdem
*
boxjenk(sdiffs=1,const,define=pharmeq) pharmdem
*
uforecast(equation=pharmeq,stderrs=efse) forecast 157 192
set lower 157 192 = forecast+%invnormal(.025)*efse
set upper 157 192 = forecast+%invnormal(.975)*efse
print(picture="*.####") 157 192 forecast efse lower upper
*
graph(footer="Figure 9-6. Forecast of animal pharmaceutical demand") 2
# pharmdem 1 156
# forecast
*
* The /3's on lower and upper force the upper and lower bounds to use the same
* color or pattern.
*
graph(footer="Figure 9-7. Prediction interval for PHARMDEM.DAT") 4
# pharmdem 132 156
# forecast
# lower / 3
# upper / 3