di _column(50) "Number of obs     =    `n'"

di "Num obs,"				_column(20) %8.0f `N_ar1'

di "Num obs,"				_column(20) %8.0f `N_fe'

di "Num obs,"		%8.0f `N_ols'

di "Mean num. observations," _column(25) %6.0f `nobs_mean'

end commented for version 2.3.1 fix of missing cluster values issue

graph twoway (line Rademacher t_stat if t_stat >= `lower' & t_stat <= `upper') (line Webb t_stat if t_stat >= `lower' & t_stat <= `upper') , xline(`main_t') ti("CDFs of Bootstrapped t-distributions") 
	note("Note: 6 Clusters.  999 bootstrap replications.  Vertical line at main t-statistic.  Zoomed in near main t-staistic.")

di in gr "Number of unique values of `varlist' is  " in ye _result(18)

	di in blu "variable `generate' contains number of unique values of `varlist' by `by'"

di "Table 4 loop.  MC Reps = `mcreps', BS Reps = `bsreps', Number of states sampled = `numstates'"

qui log off t2out

set more off

qui log off t3out

qui log off t4out

set more off

set more off

qui log off t3out

di "OLS"

	ereturn local model =             "ols"

estimates store ols_white

qui log on t3out

qui log on t3out

qui log on t1out

gen Rademacher = sum(one) / _N

gen one = 1

collapse  one , by(statefip)

gen byte one = 1

program define cgmreg, eclass byable(onecall) sortpreserve

local in "opt"

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

		di " -> Removing string 'robust' from your options line: it's unnecessary as an option,"

		local options = regexr("`options'", "robust", "")

/* get the dataset ready for estimation */
qui t3_setup , keylhs("lnwage_sy") sourcedata("${madedat}/CPS_panel") outdata(`main_data')

syntax  ,  keylhs(string)  sourcedata(string) outdata(string)