*
* Example 6.3 from page 109
* First printing of text has errors
*
open data tablef5-1[1].txt
calendar(q) 1950
data(format=prn,org=columns) 1950:1 2000:4 year qtr realgdp realcons realinvs realgovt realdpi $
  cpi_u m1 tbilrate unemp pop infl realint
*
set logc = log(realcons)
set logy = log(realdpi)
*
linreg(vcv) logc
# constant logy logc{1}
*
* The linearized variance can be computed using SUMMARIZE.
*
* The variance in the text seems to be incorrect, as the cov(b,c) used
* in the calculations there isn't correct.
*
summarize(title="Long-run MPC") %beta(2)/(1-%beta(3))
*
* Direct test of long-run MPC=1
*
restrict(title="Direct Test of Long-Run MPC=1") 1
# 2 3
# 1.0 1.0 1.0