/* drop if earnings allocated */
drop if qincwage ~= 0

graph export quantiles_`type'.png , replace

/* the main loop */
qui forvalues mc = 1/`mcreps' {

qui forvalues bb = 1/`bsreps' {

/* the main loop */
qui forvalues mc = 1/`mcreps' {

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

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)

qui reg policy i.state i.year

qui reg policy i.year

qui use statefip using "`sourcedata'"

qui use "`sourcedata'" , replace

qui use `bsout'

qui use "`mysourcedata'"

	qui predict double `resid' if `touse'==1, residual

qui gen `recnum' = _n

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

	qui by `by': replace `uniq' = sum(`uniq')

qui keep if uniform() < 0.20

qui gen rej = .

qui by statefip: gen subid = _n

qui replace `reglhs' = `reglhs'- `beta_hypothesis' * `regrhs'

qui replace `reglhs' = `reglhs'  + (`beta_hypothesis' * `regrhs' )

qui gen rej = .

qui by statefip: gen subid = _n

qui gen datetime = "$S_DATE $S_TIME" 		

		qui unique `clusvar' if `touse'

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

/* calls program "table4_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(pseudo_cluster) idcluster(pseudo_cluster_2) 
	seed(10101) saving(`pairs_bs_out' , double)
	: table4_np_boot , lhs(`lhs') rhs(`rhs') keyrhs(`keyrhs') clusterid(pseudo_cluster_2)

qui save `statelist'

qui save `to_be_sampled'

qui save `rad'

qui save "`statelist'"

qui save `tobsample'

qui save `statelist'

qui save `to_be_sampled' , replace

qui save `bsout' , replace

qui save `main_data' , replace

qui save `savedata' , replace

qui save "`mysourcedata'" , replace

qui save t2_mcout_G_`ns' , replace

qui save "`mysourcedata'" , replace

qui save "`mcout_`ns''" , replace

qui save t4_mcout_all , replace

qui log using table4_out.csv , text replace name(t4out)

qui log using table3_out.csv , text replace name(t3out)

qui summ p_`mthd'

qui summ dof_CR2_IK , detail

qui summ numobs

qui summ bs_misreps

qui summ b

qui summ rank if abs(t_`var' - `main_t') < 0.0001  

qui summ numbsreps

qui summ bs_misreps_t

qui summ b

qui log off t2out

qui log off t2out

qui log on t3out

qui postclose bs_output

qui log off t3out

qui log on t4out

qui log close t4out

qui log on t1out

qui log close t1out

qui log close t3out