predict double pricehat, xb

		scalar mu2 = 2 * xbar[1, 2]

label def ind 1 Kappa 2 Lambda 3 Nu 4 Xi 5 Tau 6 Upsilon 7 Phi 8 Chi 9 Psi

xi: regress bp i.agegrp i.sex

xi: regress bp i.agegrp*bp0

textbarplot compname beta_k, recast(dot) xline(0) ///
scheme(s2mono) nodraw name(fig9_2b,replace)

		mat eta = e(Xmfx_dyex) \ e(Xmfx_se_dyex)

// retrieve elasticities
   matrix Meyex = e(Xmfx_eyex)                         

matrix nu =  e(Xmfx_eyex)

		mat eta = e(Xmfx_dyex) \ e(Xmfx_se_dyex)

// retrieve standard errors of the elasticities
   matrix Seyex = e(Xmfx_se_eyex)                     

matrix nuse = e(Xmfx_se_eyex)

// save predictions at these points in tempvar y
   quietly replace `y' = e(Xmfx_y) in `i'   

scatter expp1992 expp1998 sal, m(oh x) scheme(s2mono) ///
  title("School District Expenditures vs. Average Salaries")  ///
  xtitle("Average salaries") ytitle("Expenditures") ///
  legend(position(5) ring(0) label(1 "1992") label(2 "1998"))