/* get the CR2, CR3 standard errors adjustments.  Also get I-K degrees of freedom, and CSS effective number of clusters */
CR23_IK_CSS , betavec(`betavec') lhs(`lhs') rhs(`rhs') key_rhs(`key_rhs') main_data("`main_data'")

/* 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'

/* the t3_ols code will partial out state and year fixed effects */
t3_ols , bsreps(`bsreps') sourcedata(`main_data') lhs(lnwage_sy) rhs(policy) keyrhs(policy)

/* the t3_ar1 code does not partial out any fixed effects, so need to directly pass these in as regressors */
t3_ar1 , bsreps(`bsreps') sourcedata(`main_data') lhs(lnwage_sy) rhs(policy i.statefip i.year) keyrhs(policy)

local lhs = "lnwage"

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

syntax  , lhs(string) rhs(string) keyrhs(string) clusterid(string) [ yearonly ]

/* make deviations from state and year Fixed Effects */
egen `lhs'_avg_s = mean(`lhs') , by(pseudo_cluster)     

gen resid_restricted = `lhs'_deviation

/* transform data to be "deviations from state and year dummies".  Then run cluster-robust regression */
cap drop `lhs'_deviation `rhs'_deviation

local shortcommand = subinstr("`command'","`lhs'_deviation","y_wild_rademacher_restricted",.)

local reglhs = "`lhs'_deviation"

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 " -> Removing string 'robust' from your options line: it's unnecessary as an option,"

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

	tempname N_clust N df_m F r2 mss rss r2_a ll ll_0 

	if "`if'"=="" local if "if 1"

	ereturn local eigenfix  = 	  "`eigenfix'"

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

	ereturn local cmd =               "cgmreg"

	ereturn local predict =           "regres_p"

	ereturn local estat_cmd =         "regress_estat"

	ereturn local wtype =             "`wtype'"

	ereturn local clusvar  =          "`cluster'"

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