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