/* main regression */
	qui regress `anything' if `touse' `weight', `options' mse1

areg lnwage age age2 sex yrseduc [pw=wtsupp] , a(styr)

reg lnwage policy age age2 yrseduc  , robust

/* the hypothesis is: beta on policy = 0.  So to estimtate imposing this restriction, we will just drop "policy" from the regression */
reg lnwage age age2 yrseduc

reg y_wild_rademacher_restricted policy age age2 yrseduc  , cluster(statefip)

reg lnwage policy age age2 yrseduc  , robust

qui bootstrap _b _se , reps(`bsreps') nodots cluster(pseudo_cluster) idcluster(pseudo_cluster_2) saving(`pairs_bs_out' , double)
	: reg lnwage policy age age2 yrseduc  , cluster(pseudo_cluster_2)

/* note that vce(robust) will do cluster robust on the i-level (statefip) cluster */
xtreg `lhs' `rhs'   , pa corr(ar1) vce(robust) noconstant

/* generate deviations-from state and year FE variables for LHS and RHS.
	use those transformed variables in later regressions */
qui reg `keylhs' i.state i.year

/* lnwage, lnwage_sy */
qui reg `keylhs' i.year

/* Row 1, default standard erorrs */
reg `reglhs' `regrhs' , vce(ols) noconstant

/* Row 3, Cluster on state */
reg `reglhs' `regrhs' , vce(cluster statefip) noconstant

cgmreg `reglhs' `regrhs' , cluster(statefip year) noconstant

/* get standard errors for rows 2-4 */
reg `reglhs' `regrhs' , noconstant cluster(pseudo_cluster)

reg lnwage policy age age2 yrseduc  , robust

reg lnwage policy age age2 yrseduc  , vce(bootstrap , cluster(statefip) reps(399) seed(10101) )

xtreg lnwage policy age age2 yrseduc , re

xtreg lnwage policy age age2 yrseduc , pa corr(exch) robust