Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • coefplot title related

    Dear Stata users

    I am doing some plots using
    Code:
    coefplot
    . The code is below:

    Code:
    preserve
        summ     t if event_time == -1
        local     true_neg = r(mean)
        reghdfe        log_ref_dur delta i.t#c.delta o`true_neg'.t#o.delta ///
                    LWbas_male age LWbas_married LWbas_immigrant LWbas_child LWbas_pre LWbas_pre9 ///
                    LWbas_s_high LWbas_l_high LWbas_s_univ LWbas_l_univ LW5_emp_mon_last5y ///
                    LWhist_empbef1 LW5_av_dur LW5_last_dur LWhist_income_bef1 LW5_welfare5, ///
                    cluster(referal_CWid#referal_year) absorb(referal_year referal_CWid llm_kod)
        estimates store log_ref_dur
    
        coefplot (log_ref_dur, msymbol(Sh) msize(medium) label(log_ref_dur)),  level(90) ///
        keep(10.t#c.delta 11.t#c.delta 12.t#c.delta ///
        13.t#c.delta 14.t#c.delta 15.t#c.delta 16.t#c.delta 17.t#c.delta) ///
        rename(10.t#c.delta= "t=-4" 11.t#c.delta= "t=-3" 12.t#c.delta= "t=-2" ///
        13.t#c.delta ="t=-1" 14.t#c.delta= "t=0" 15.t#c.delta= "t=1" 16.t#c.delta ="t=2" 17.t#c.delta= "t=3") ///
        omitted baselevels xline(4, lcolor(gs5) lp(dash) lwidth(vlarge)) ytitle("Log duration from 1st referral", size(vlarge)) ///
        yline(0, lcolor(black) lp(dash) lwidth(vlarge))  name(d, replace) xtitle("Event time", size(vlarge)) leg(off) ///
        legend(size(vlarge)) ylabel(-2(1)3, labsize(vlarge)) xlabel(, labsize(vlarge)) recast(connected) ciopts(recast(rcap)) lcolor(red) lp(dash) vertical
    restore

    Code:
    preserve
        summ     t if event_time == -1
        local     true_neg = r(mean)
        reghdfe        ref_dur90 delta i.t#c.delta o`true_neg'.t#o.delta ///
                    LWbas_male age LWbas_married LWbas_immigrant LWbas_child LWbas_pre LWbas_pre9 ///
                    LWbas_s_high LWbas_l_high LWbas_s_univ LWbas_l_univ LW5_emp_mon_last5y ///
                    LWhist_empbef1 LW5_av_dur LW5_last_dur LWhist_income_bef1 LW5_welfare5, ///
                    cluster(referal_CWid#referal_year) absorb(referal_year referal_CWid llm_kod)
        estimates store ref_dur90
    
        coefplot (ref_dur90, msymbol(Sh) msize(medium) label(ref_dur90)),  level(90) ///
        keep(10.t#c.delta 11.t#c.delta 12.t#c.delta ///
        13.t#c.delta 14.t#c.delta 15.t#c.delta 16.t#c.delta 17.t#c.delta) ///
        rename(10.t#c.delta= "t=-4" 11.t#c.delta= "t=-3" 12.t#c.delta= "t=-2" ///
        13.t#c.delta ="t=-1" 14.t#c.delta= "t=0" 15.t#c.delta= "t=1" 16.t#c.delta ="t=2" 17.t#c.delta= "t=3") ///
        omitted baselevels xline(4, lcolor(gs5) lp(dash) lwidth(vlarge)) ytitle("Pr(duration<90 days)", size(large)) ///
        yline(0, lcolor(black) lp(dash) lwidth(vlarge))  name(e, replace) xtitle("Event time", size(vlarge)) leg(off) ///
        legend(size(vlarge)) ylabel(-.6(.2).4, labsize(vlarge)) xlabel(, labsize(vlarge)) recast(connected) ciopts(recast(rcap)) lcolor(red) lp(dash) vertical
    restore
    The corresponding plots are below:
    Click image for larger version

Name:	twfe_eventstudy_logrefdur_Norisk_time_variant.png
Views:	1
Size:	39.7 KB
ID:	1776195

    Click image for larger version

Name:	twfe_eventstudy_refdur90_Norisk_time_variant.png
Views:	1
Size:	41.0 KB
ID:	1776196




    I don't have any clue regarding why there is a large gap between ytitle and yaxis in the second graph. How can I reduce this gap to have plot like figure 1? Any help would be much appreciated.

    Thanks in advance,

    Zariab
    Uppsala University

  • #2
    coefplot is from SSC, as you are asked to explain (FAQ Advice #12). Please provide a reproducible example.

    Comment

    Working...
    X