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

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

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

gen Webb = sum(one) / _N

graph twoway (line Rademacher t_stat) (line Webb t_stat) , xline(`main_t') ti("CDFs of Bootstrapped t-distributions") 
	note("Note: 6 Clusters.  999 bootstrap replications.  Vertical line at main t-statistic.")

sort t_stat Rademacher Webb

di "percentile T Webb6, ,"											_column(35) %6.4f `p_webb_res'

foreach mthd in rob_Tnk clu_Tnk clu_Tg1 clu_Tg2 CR2_Tg1 CR3_Tg1 CR2_Tik clu_Tcss bsse_Tg1 pairs_bs rad_res rad_res_high rad_res_low webb_res rad_unres {

foreach var in pairs_bs rad_res rad_unres webb_res {

foreach mthd in def_Tnk clu_Tnk clu_Tg1 clu_Tg2 bsse_Tg1 pairs_bs rad_res webb_res {

foreach var in pairs_bs rad_res webb_res {

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

di "se CR2 wIK,"			_column(20) %8.4f `se_fe_CR2' "," 		_column(35) %6.4f `p_CR2_Tik'

keep statefip wild_rademacher wild_webb

gen wild_rademacher = -1 + 2 * (my_uniform >= 0.5)

	/* create transformed residuals and new wild-outcome-variables */
gen resid_wild_rad_restricted = resid_restricted * wild_rademacher

gen wild_rademacher = -1 + 2 * (my_uniform >= 0.5)

	/* create transformed residuals and new wild-outcome-variables */
gen resid_wild_rad_restricted = resid_restricted * wild_rademacher

keep pseudo_cluster wild_rademacher wild_webb

gen wild_webb = 	(-1) * sqrt(1.5) * (my_uniform > (0) & my_uniform <= (1/6)) +  
						(-1) * sqrt(1) * (my_uniform > (1/6) & my_uniform <= (2/6))  + 
						(-1) * sqrt(0.5) * (my_uniform > (2/6) & my_uniform <= (3/6)) + 
						(+1) * sqrt(0.5) * (my_uniform > (3/6) & my_uniform <= (4/6)) + 
						(+1) * sqrt(1) * (my_uniform > (4/6) & my_uniform <= (5/6))  + 
						(+1) * sqrt(1.5) * (my_uniform > (5/6) & my_uniform <= (6/6))

gen resid_wild_webb_restricted = resid_restricted * wild_webb

keep pseudo_cluster wild_rademacher wild_webb

gen wild_webb = 	(-1) * sqrt(1.5) * (my_uniform > (0) & my_uniform <= (1/6)) +  
						(-1) * sqrt(1) * (my_uniform > (1/6) & my_uniform <= (2/6))  + 
						(-1) * sqrt(0.5) * (my_uniform > (2/6) & my_uniform <= (3/6)) + 
						(+1) * sqrt(0.5) * (my_uniform > (3/6) & my_uniform <= (4/6)) + 
						(+1) * sqrt(1) * (my_uniform > (4/6) & my_uniform <= (5/6))  + 
						(+1) * sqrt(1.5) * (my_uniform > (5/6) & my_uniform <= (6/6))

gen resid_wild_webb_restricted = resid_restricted * wild_webb

keep pseudo_cluster wild_rademacher wild_webb

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

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

	ereturn local title =             "Linear regression with CGM covariance estimation"

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

/* generate ln(earnings/hour) */
gen numweeks = 43.5 if WKSWORK2 == 4

replace numweeks = 50 if WKSWORK2 == 6

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

	ereturn local wtype =             "`wtype'"