/* *** This do file creates the replication results for Implementing Multilateral Environmental Agreements: An Analysis of EU Directives (with Richard Perkins), Global Environmental Politics, 7 (3), 2007, pp. 13-31 Eric Neumayer (LSE) Please report any errors to e.neumayer@lse.ac.uk */ use "Article for GEP (EU compliance).dta", clear * Table 2 quietly xi: nbreg reasop gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears i.year, robust su reasop gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears if e(sample) * Table 3 xi: nbreg reasop gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears i.year, robust xi: nbreg reasop lreasop gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears i.year, robust * Non-Nested AIC and BIC results global XLISTMODEL1 gpiflorianklein manufact firsttwoyears i.year global XLISTMODEL2 lnpop lneumembershiplength firsttwoyears i.year global XLISTMODEL3 polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch firsttwoyears i.year global XLISTMODEL4 eupositivetrendfile meashare firsttwoyears i.year xi: nbreg reasop $XLISTMODEL1 estimates store model1 scalar ll1 = e(ll) scalar q1 = e(k) scalar N1 = e(N) scalar aic1 = -2*ll1 + 2*q1 scalar bic1 = -2*ll1 + ln(N1)*q1 scalar caic1 = -2*ll1 + (1 + ln(N1))*q1 xi: nbreg reasop $XLISTMODEL2 if e(sample) estimates store model2 scalar ll2 = e(ll) scalar q2 = e(k) scalar N2 = e(N) scalar aic2 = -2*ll2 + 2*q2 scalar bic2 = -2*ll2 + ln(N2)*q2 scalar caic2 = -2*ll2 + (1 + ln(N2))*q2 xi: nbreg reasop $XLISTMODEL3 if e(sample) estimates store model3 scalar ll3 = e(ll) scalar q3 = e(k) scalar N3 = e(N) scalar aic3 = -2*ll3 + 2*q3 scalar bic3 = -2*ll3 + ln(N2)*q3 scalar caic3 = -2*ll3 + (1 + ln(N2))*q3 xi: nbreg reasop $XLISTMODEL4 if e(sample) estimates store model4 scalar ll4 = e(ll) scalar q4 = e(k) scalar N4 = e(N) scalar aic4 = -2*ll4 + 2*q4 scalar bic4 = -2*ll4 + ln(N2)*q4 scalar caic4 = -2*ll4 + (1 + ln(N2))*q4 estimates table model1 model2 model3, stats(N k ll aic bic) di "Model 1: " _n "lnL: " ll1 " q: " q1 _n " N: " N1 di "-2lnL: " -2*ll1 _n "AIC: " aic1 _n " BIC: " bic1 _n "caic: " caic1 di "Model 2: " _n "lnL: " ll2 " q: " q2 _n " N: " N2 di "-2lnL: " -2*ll2 _n "AIC: " aic2 _n " BIC: " bic2 _n "caic: " caic2 di "Model 3: " _n "lnL: " ll3 " q: " q3 _n " N: " N3 di "-2lnL: " -2*ll3 _n "AIC: " aic3 _n " BIC: " bic3 _n "caic: " caic3 di "Model 4: " _n "lnL: " ll4 " q: " q4 _n " N: " N4 di "-2lnL: " -2*ll4 _n "AIC: " aic4 _n " BIC: " bic4 _n "caic: " caic4 * Nested AIC and BIC results global XLISTMODEL1 gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears i.year global XLISTMODEL2 polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears i.year global XLISTMODEL3 gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch eupositivetrendfile meashare firsttwoyears i.year global XLISTMODEL4 gpiflorianklein manufact eupositivetrendfile meashare lnpop lneumembershiplength firsttwoyears i.year global XLISTMODEL5 gpiflorianklein manufact polcon3 leg_british leg_french leg_german bureaucracyquality rgdpch lnpop lneumembershiplength firsttwoyears i.year xi: nbreg reasop $XLISTMODEL1 estimates store model1 scalar ll1 = e(ll) scalar q1 = e(k) scalar N1 = e(N) scalar aic1 = -2*ll1 + 2*q1 scalar bic1 = -2*ll1 + ln(N1)*q1 scalar caic1 = -2*ll1 + (1 + ln(N1))*q1 xi: nbreg reasop $XLISTMODEL2 if e(sample) estimates store model2 scalar ll2 = e(ll) scalar q2 = e(k) scalar N2 = e(N) scalar aic2 = -2*ll2 + 2*q2 scalar bic2 = -2*ll2 + ln(N2)*q2 scalar caic2 = -2*ll2 + (1 + ln(N2))*q2 xi: nbreg reasop $XLISTMODEL3 if e(sample) estimates store model3 scalar ll3 = e(ll) scalar q3 = e(k) scalar N3 = e(N) scalar aic3 = -2*ll3 + 2*q3 scalar bic3 = -2*ll3 + ln(N2)*q3 scalar caic3 = -2*ll3 + (1 + ln(N2))*q3 xi: nbreg reasop $XLISTMODEL4 if e(sample) estimates store model4 scalar ll4 = e(ll) scalar q4 = e(k) scalar N4 = e(N) scalar aic4 = -2*ll4 + 2*q4 scalar bic4 = -2*ll4 + ln(N2)*q4 scalar caic4 = -2*ll4 + (1 + ln(N2))*q4 xi: nbreg reasop $XLISTMODEL5 if e(sample) estimates store model5 scalar ll5 = e(ll) scalar q5 = e(k) scalar N5 = e(N) scalar aic5 = -2*ll5 + 2*q5 scalar bic5 = -2*ll5 + ln(N2)*q5 scalar caic5 = -2*ll5 + (1 + ln(N2))*q5 estimates table model1 model2 model3, stats(N k ll aic bic) di "Model 1: " _n "lnL: " ll1 " q: " q1 _n " N: " N1 di "-2lnL: " -2*ll1 _n "AIC: " aic1 _n " BIC: " bic1 _n "caic: " caic1 di "Model 2: " _n "lnL: " ll2 " q: " q2 _n " N: " N2 di "-2lnL: " -2*ll2 _n "AIC: " aic2 _n " BIC: " bic2 _n "caic: " caic2 di "Model 3: " _n "lnL: " ll3 " q: " q3 _n " N: " N3 di "-2lnL: " -2*ll3 _n "AIC: " aic3 _n " BIC: " bic3 _n "caic: " caic3 di "Model 4: " _n "lnL: " ll4 " q: " q4 _n " N: " N4 di "-2lnL: " -2*ll4 _n "AIC: " aic4 _n " BIC: " bic4 _n "caic: " caic4 di "Model 5: " _n "lnL: " ll5 " q: " q5 _n " N: " N5 di "-2lnL: " -2*ll5 _n "AIC: " aic5 _n " BIC: " bic5 _n "caic: " caic5