use "rates80_08.dta", clear rename marrrate marrate gen pop2=round(pop/1000) gen t=year-1980 gen t2=t^2 sort state year by state: gen gonrate_1=gonrate[_n-1] by state: gen syprate_1=syprate[_n-1] by state: gen marrate_1=marrate[_n-1] gen marratet=marrate_1*t gen marratet2=marrate_1*t2 gen repeal=bldtest-1 replace repeal=1 if repeal==-1 sort state year by state: gen dropind=repeal[_n-1] drop if dropind==1 ***b/c won't have lags, and don't want states w/pre-80 repeal: drop if year==1980 *** probit repeal marrate_1 if state~=12 & state~=29 & state~=5, cluster(state) r est store A mfx compute probit repeal marrate_1 gonrate_1 syprate_1 if state~=12 & state~=29 & state~=5, cluster(state) r est store B mfx compute probit repeal marrate_1 gonrate_1 syprate_1 t t2 marratet marratet2 if state~=12 & state~=29 & state~=5, cluster(state) r est store C test marrate_1 gonrate_1 syprate_1 marratet marratet2 mfx compute est table A B C, b(%9.4f) se(%9.4f)