mat `b' = e(b)

*	scalar N_clust =	e(N_clust)
	scalar `N' =           	e(N) 

*	scalar `df_r' =         e(df_r)
*	scalar `F' =            e(F)
	scalar `r2' =           e(r2)

	scalar `rss' =          e(rss)

	local title =           e(title)

	local predict =         e(predict)

	local estat_cmd =       e(estat_cmd)

	local n = e(N)

	ereturn post `b' `running_sum' , e(`touse') depname(`depname') 

local K = e(rank)

matrix `betavec' = e(b)

local N_ar1 = e(N)

local N_fe = e(N)

local G_fe = e(N_clust)

local bs_misreps = e(N_misreps)

local numobs = e(N)

local N_ols = e(N) 						

local N_re = e(N)

foreach vv in educ age sex year state incwage {

replace yrseduc = 3 if educ == 10

replace yrseduc = 2 if educ == 12

replace yrseduc = 4 if educ == 14

replace yrseduc = 5 if educ == 21

replace yrseduc = 8 if educ == 30

replace yrseduc = 8 if educ == 32

replace yrseduc = 10 if educ == 50

replace yrseduc = 11 if educ == 71

replace yrseduc = 12 if educ == 73

replace yrseduc = 14 if educ == 81

replace yrseduc = 14 if educ == 91

replace yrseduc = 15 if educ == 100

replace yrseduc = 16 if educ == 111

replace yrseduc = 16 if educ == 122

replace yrseduc = 19 if educ == 124

di "CSS effctv clustrs,"	_column(20) "-"  ","

di "Fixed Effects, `lhs', `keyrhs'"

	local eigenfix "no"

	ereturn local eigenfix  = 	  "`eigenfix'"

		    di "          -cgmreg- is replacing any/all negative eigenvalues with 0."

	mat `running_sum' = `eigenvectors'*diag(`eigenvalues')*`eigenvectors''

			local element = `element'[1,1]

				di "          Because you used the -noeigenfix- option, -cgmreg- must end."

	local estat_cmd =       e(estat_cmd)

		    di "          Raw, non-psd covariance estimate will be available "

		    di " -> NOTE: Raw estimated variance matrix was non-positive semi-definite."

	mat `rows' = rowsof(e(V))

local my_top_path = subinstr(c(pwd), "doug work\Colin\JHR cluster paper\work Summer 2014\empirical example code", "",.)

global rawdat = "`my_top_path'data\CPS\July 2013 extract for cluster JHR"