*
* Example 6.5.5 from page 210
*
open data deaths.dat
calendar(m) 1973
data(format=free,org=columns) 1973:1 1978:12 deaths
*
boxjenk(diff=1,sdiff=1,demean,ma=1,sma=1,maxl,define=eq658) deaths
uforecast(equation=eq658,stderrs=s658) f658 1979:1 1979:6
*
boxjenk(diff=1,sdiff=1,demean,ma=||1,6,12,13||,maxl,define=eq659) deaths
uforecast(equation=eq659,stderrs=s659) f659 1979:1 1979:6
*
print(picture="*.") / f658 s658 f659 s659