*
* ANOVA table from pp 33-34
*
open data tablef3-1[1].txt
calendar 1968
data(format=prn,org=columns) 1968:1 1982:1 year gnp invest cpi interest
*
set rgnp = gnp/cpi*(.1)
set rinv = invest/cpi*(.1)
set trend = t
*
* Compute the inflation rate. Given the data on tablef3-1, we can't get an
* inflation value for the first observation, so that gets copied from table 3.1.
*
set infl  = 100*(cpi/cpi{1}-1)
set infl 1968:1 1968:1 = 4.40
*
linreg rinv
# constant trend rgnp interest infl
*
* The REGANOVA procedure can be used after running a regression to get an
* analysis of variance table like Table 3.4.
*
@reganova