gen p_rob_Tnk = 2 * ttail(dof_nmk,(abs(b - 0) / se_rob))

gen p_clu_Tg1 = 2 * ttail(dof_gm1,(abs(b - 0) / se_clu))

gen p_CR2_Tg1 = 2 * ttail(dof_gm1,(abs(b - 0) / se_CR2))

gen p_CR2_Tik = 2 * ttail(dof_CR2_IK,(abs(b - 0) / se_CR2))

gen p_bsse_Tg1 = 2 * ttail(dof_gm1,(abs(b - 0) / se_clu_bs))

gen p_clu_Tnk = 2 * ttail(dof_nmk,(abs(b - 0) / se_clu))

gen p_clu_Tg2 = 2 * ttail(dof_gm2,(abs(b - 0) / se_clu))

summ rank if abs(t_`var' - `main_t') < 0.0001  

summ rank if abs(t_`var' - `main_t') < 0.0001  

local p_clu_Tg1 = 2 * ttail(`dof_gm1',(abs(`b_ar1' - `beta0') / `se_ar1_clu'))

*set trace on ;
local p_def_Tnk = 2 * ttail(`dof_nmk',(abs(`b_fe' - `beta0') / `se_fe_def'))

local p_clu_Tg1 = 2 * ttail(`dof_gm1',(abs(`b_fe' - `beta0') / `se_fe_clu'))

local p_CR2_Tik = 2 * ttail(`DOF_fe_IK',(abs(`b_fe' - `beta0') / `se_fe_CR2'))

	di "Note: +/- means the corresponding matrix is added/subtracted"

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

		    di "            in which case -cgmreg- will throw an error)"

gen wage_per_hour = incwage1999 / annual_hours

	syntax anything [if] [in] [aweight fweight iweight pweight /], /*
		*/ Cluster(string) [NOEIGenfix *] 

xtreg `lhs' `rhs'  , pa corr(ar1) vce(conventional) noconstant

xtreg `lhs' `rhs'  , pa corr(ar1) vce(conventional) noconstant

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

merge 1:m statefip using "`sourcedata'" , assert(match using) keep(match) nogenerate

merge 1:m pseudo_cluster using `main_data' , assert(match) keep(match) nogenerate

merge 1:m statefip using "`sourcedata'" , assert(match using) keep(match) nogenerate

merge 1:m pseudo_cluster using `main_data' , assert(match) keep(match) nogenerate

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 "Avg miss prs BS reps, " _column(25) %6.0f  `mean_bsreps_miss'

di "Avg Rad Res p-high, " _column(25) %6.4f `mean_p_rad_res_high'

di "Avg Rad Res p-avg, " _column(25) %6.4f `mean_p_rad_res'

	syntax [if] [in] [aweight fweight iweight pweight /] , groupvar(string) xxinv(string) plusminus(string) resid(string) running_sum(string) touse(string)