gen y_wild_rad_restricted = yhat_restricted + resid_wild_rad_restricted

egen y_wild_rad_restricted_avg_y = mean(y_wild_rad_restricted) , by(year)

	/* to do this, make year and state specific averages, then call "make_deviations" code */
cap drop y_wild_rad_restricted_avg_s

cap drop y_wild_rad_restricted_avg_y

local shortcommand = subinstr("`command'","`reglhs'","y_wild_rad_restricted_deviation",.)

gen y_wild_rademacher_restricted = yhat_restricted + resid_wild_rad_restricted

gen y_wild_rademacher_restricted = yhat_restricted + resid_wild_rad_restricted

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

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

gen y_wild_webb_restricted = yhat_restricted + resid_wild_webb_restricted

gen y_wild_webb_restricted = yhat_restricted + resid_wild_webb_restricted

egen y_wild_webb_restricted_avg_s = mean(y_wild_webb_restricted) , by(pseudo_cluster)     

make_deviations y_wild_rad_restricted y_wild_webb_restricted

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

egen y_wild_webb_restricted_avg_s = mean(y_wild_webb_restricted) , by(pseudo_cluster)     

egen y_wild_webb_restricted_avg_y = mean(y_wild_webb_restricted) , by(year)

egen yavg_ds_`var' = mean(ds_`var') , by(year)

cap drop ds_`var' yavg_ds_`var'

di in gr "Number of records is  "in ye "`N'"

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

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

	if "`eigenfix'"=="yes" {

/* generate residuals.  Also, generate "restricted, impose the null hypothesis" residuals */
predict yhat , xb

predict yhat_restricted , xb

gen y_wild_webb_restricted = yhat_restricted + resid_wild_webb_restricted

replace yhat_restricted = yhat_restricted + (`beta_hypothesis' * `key_rhs')

gen y_wild_webb_restricted = yhat_restricted + resid_wild_webb_restricted

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

gen y_wild_webb_restricted = yhat_restricted + resid_wild_webb_restricted

gen y_wild_rad_restricted = yhat_restricted + resid_wild_rad_restricted

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

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

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

local yrmean = r(mean)