*
* Example 6.2.2 from page 193
*
open data lake.dat
cal 1875
data(format=free,org=columns) 1875:1 1972:1 lake
*
boxjenk(ar=1,ma=1,demean,maxl) lake
@bdindtests %resids
diff(center) lake / clake
@bjautofit(pmax=5,qmax=5) clake