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

savedresults compare ivhomo e(), include(macros: depvar scalar: N j matrix: b V) ///
             tol(1e-7) verbose

savedresults compare ivrobust e(), include(macros: depvar scalar: N j matrix: b V) /// 
             tol(1e-7) verbose

savedresults compare ivreg2cue e(), include(macros: depvar scalar: N j matrix: b V) ///
             tol(1e-4) verbose

savedresults compare ivreg2cue e(), include(macros: depvar scalar: N j matrix: b V) ///
             tol(1e-4) verbose

matrix list e(b)

matrix b = e(b)

matrix V = e(V)

// retrieve elasticities
   matrix Meyex = e(Xmfx_eyex)                         

local regressors: colnames e(b)

matrix nu =  e(Xmfx_eyex)

matrix nu =  e(Xmfx_eyex)

	matrix beta`i' = e(b)

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

		mat xbar = e(Xmfx_X)

// compute appropriate t-statistic  for 95% confidence interval
scalar tmfx = invttail(e(df_r), 0.975)

  cap noi spinst_itsp_wrk egenmore http://fmwww.bc.edu/repec/bocode/e

esttab one two, mtitles("OLS" "IV") stat(nu nuse, labels("Elas.(pop)" "SE"))

	esttab _all using lwage`y'.tex, replace ti("Wage equations for `y'") ///
	       nomtitles nodepvars not se noobs ar2 booktabs ///
	       scalar(eta_educ etase exper se) ///
	       addnotes("eta_educ: semi-elasticity of lwage with respect to educ" ///
	       "etase: standard error of the semi-elasticity" ///
	       "exper: effect of exper on lwage at mean exper") ///
	       substitute("_cons" "Constant" "eta_educ" "$\eta_{educ}$" ///
	       "etase" "$\eta_{se}$") 

mygmm2s lw, endog(iq) inexog(s expr tenure rns smsa) ///
            exexog(med kww age mrt)

mygmmcue lw, endog(iq) inexog(s expr tenure rns smsa) ///
            exexog(med kww age mrt)

label var end Endpoint

file write `hh' "\hline" _n "\end{tabular}" "\medskip" _newline

keep if eoq

	esttab _all using lwage`y'.tex, replace ti("Wage equations for `y'") ///
	       nomtitles nodepvars not se noobs ar2 booktabs ///
	       scalar(eta_educ etase exper se) ///
	       addnotes("eta_educ: semi-elasticity of lwage with respect to educ" ///
	       "etase: standard error of the semi-elasticity" ///
	       "exper: effect of exper on lwage at mean exper") ///
	       substitute("_cons" "Constant" "eta_educ" "$\eta_{educ}$" ///
	       "etase" "$\eta_{se}$") 

	esttab _all using lwage`y'.tex, replace ti("Wage equations for `y'") ///
	       nomtitles nodepvars not se noobs ar2 booktabs ///
	       scalar(eta_educ etase exper se) ///
	       addnotes("eta_educ: semi-elasticity of lwage with respect to educ" ///
	       "etase: standard error of the semi-elasticity" ///
	       "exper: effect of exper on lwage at mean exper") ///
	       substitute("_cons" "Constant" "eta_educ" "$\eta_{educ}$" ///
	       "etase" "$\eta_{se}$") 

estadd scalar nu = nu[1, 1]

estadd scalar nu = nu[1, 1]

		eststo, addscalars(eta_educ eta[1, 1] etase eta[2, 1] exper r(estimate) se r(se))

  cap noi spinst_itsp_wrk estout http://fmwww.bc.edu/repec/bocode/e

eststo: quietly regress so2 temp manuf pop

eststo: quietly regress so2 temp wind precip days

esttab using esttab_example.tex, ar2 label se nostar nodepvars brackets ///
nomtitles title("Models of sulphur dioxide concentration") booktabs ///
alignment(D{.}{.}{-1}) replace

   matrix eta = Meyex[1, "stratio"]             // for the stratio column

		eststo, add(eta_educ eta[1, 1] etase eta[2, 1])

		eststo, add(eta_educ eta[1, 1] etase eta[2, 1])

		eststo, addscalars(eta_educ eta[1, 1] etase eta[2, 1] exper r(estimate) se r(se))

		eststo, add(eta_educ eta[1, 1] etase eta[2, 1])

		eststo, addscalars(eta_educ eta[1, 1] etase eta[2, 1] exper r(estimate) se r(se))

		eststo, add(eta_educ eta[1, 1] etase eta[2, 1])

		eststo, addscalars(eta_educ eta[1, 1] etase eta[2, 1] exper r(estimate) se r(se))

save ex4.1, replace

nlcom -exp([a]_cons)

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"))

Bolivia exports 80 115 

mfx, eyex

mfx compute, eyex

tempvar y x eyex seyex1 seyex2