egen statefip2 = group(statefip id)

/* calls program "table3_np_boot", which is defined in t34_programs.do.  That program re-does the "deviations from means" and then 
re-runs the regression, to get the betas and (clustered) standard errors. */
qui bootstrap b=r(mybeta) se=r(myse) , reps(`bsreps') nodots cluster(statefip) idcluster(pseudo_cluster_2) 
	seed(10101) saving(`pairs_bs_out' , double)
	: table3_np_boot , lhs(`lhs') rhs(`rhs') keyrhs(`keyrhs') clusterid(pseudo_cluster_2) `yearonly'

	if "`if'"~="" di _column(50) "If condition      =    `if'"

		    di "          (If you don't want this behavior, use the 'noeigenfix' option,"

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

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

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

local IK_dof_p5 = r(p5)

	if "`in'"~="" di _column(50)     "In condition      =    `in'"

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

		di "    If some variable in your options list contains the string 'robust', you will"

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

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

di "policy turns on starting in year >= `yrmean'"

				local included= `included' + 1

drop if incwage == 0

gen incwage1999 = incwage * CPI99

		    		di "            'Robust Inference with Multi-Way Clustering'."

	qui by `by' `varlist' `touse': gen int `count' = _N if _n == 1

syntax  , [mcreps(integer 3) bsreps(integer 4) numstates(integer 6) sourcedata(string) savedata(string)]

syntax  , [bsreps(integer 3) sourcedata(string) lhs(string) rhs(string) keyrhs(string) beta0(integer 0) yearonly]

di "intraclass correation = `rho'"

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

di in gr "Number of records is  "in ye "`N'"

		di "    but it can cause problems if we leave it in."

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

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