Hi and thank you so much for your response. I tried to incorporate your code into mine and it is working somewhat. However, I am not sure how to interpret. My logit model is as follows: proc surveylogistic data=FINALCITYDATA3_means_wins; cluster auditor_fkey; model Switch(event='1') = sanc_PY post_violation sanc_PYPost cpa_acm cpaPost growth absdacc invar gc modop icw tenure roa loss leverage chglev cash size chgsize m_a abnrml_fee yfe2006 yfe2007 yfe2008 yfe2009 yfe2010 yfe2011 yfe2012 yfe2013 yfe2014 yfe2015 sic_2 sic_7 sic_8 sic_9 sic_10 sic_12 sic_13 sic_14 sic_15 sic_16 sic_17 sic_20 sic_21 sic_22 sic_23 sic_24 sic_25 sic_26 sic_27 sic_28 sic_29 sic_30 sic_31 sic_32 sic_33 sic_34 sic_35 sic_36 sic_37 sic_38 sic_39 sic_40 sic_41 sic_42 sic_44 sic_45 sic_46 sic_47 sic_48 sic_49 sic_50 sic_51 sic_52 sic_53 sic_54 sic_55 sic_56 sic_57 sic_58 sic_59 sic_70 sic_72 sic_73 sic_75 sic_78 sic_79 sic_80 sic_82 sic_83 sic_87 sic_89/rsq DF=infinity; output out = SwitchwCPA_CITYPY p = prob xbeta = logit; run; in the paper that I am following the first thing they compare is the coefficient on the intercept to the coefficient on the intercept plus the sanc_PY variable . Obviously, the difference between those two is the sanc_PY variable. They are looking at if that difference is significant or not. Does the code you provide accomplish this? Sorry, but this is beyond my level of expertise. I am trying! Thank you again!!
... View more