local K_fe = e(rank)

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

/* Rows 4 and 5 */
CR23_IK_CSS , betavec(`betavec') lhs(`reglhs') rhs(`regrhs') key_rhs(`regrhs') 
	noconstant(1) main_data("`main_data'")

local key_rhs = "policy"

local key_rhs = "`keyrhs'"

/* get the dataset ready for estimation */
qui t3_setup , keylhs("lnwage_sy") sourcedata("${madedat}/CPS_panel") outdata(`main_data')

/* generate state and year averages of the LHS variable and the policy variable */
egen `keylhs'_avg_s = mean(`keylhs') , by(statefip)

predict `keylhs'_resid_styr , resid

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

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 ]

egen `keyrhs'_avg_s = mean(`keyrhs') , by(pseudo_cluster)