*
* Example 5.5.1 from page 171
*
open data lake.dat
cal 1875
data(format=free,org=columns) 1875:1 1972:1 lake
report(action=define,hlabels=||"p","sigma**2","FPE"||)
do p=0,10
   boxjenk(ar=p,maxl,demean,noprint) lake
   report(atcol=1) p %rss/%nobs %rss/%nobs*(%nobs+p)/(%nobs-p)
end do p
report(action=format,atcol=2,tocol=3,picture="*.####")
report(action=format,tag=minimum,atcol=3,align=decimal,special=one)
report(action=show)