*
*  Example 5.6
*  AKAIKE.PRG
*  Akaike Information Criterion, Schwarz Criterion for Distributed Lags
*
cal 1947 1 12
all 1999:3
open data rates.rat
data(format=rats)
*
cmom
# constant shortrate{0 to 24} longrate
report(action=define,hlabels=||'Lags','Akaike','Schwarz'||)
do maxlag=0,24
   linreg(cmom,noprint) longrate
   # constant shortrate{0 to maxlag}
   compute akaike =log(%rss/%nobs)+%nreg*2.0/%nobs
   compute schwarz=log(%rss/%nobs)+%nreg*log(%nobs)/%nobs
   report(row=new,atcol=1) maxlag akaike schwarz
end do
report(action=format,picture='*.###')
report(action=show)