Coefplot - multiple models with different dependent variables

[amandacpac]

All, I am running multiple regressions with the same independent variable and 29 different dependent variable, and then running these same models multiple times for different time points. I would like to get a plot of the regression coefficients where the 29 dependent variables are along the Y axis, and multiple plots, side-by-side for the different time points. I found a post where one of the comments was getting at something similar, but I can't figure out how to get it the way I need it.
So, using this code:

sysuse auto, clear  
reg trunk foreign gear_ratio if rep78 == 5  // Outcome = Trunk and sub-group = 4
estimates store trunk_s5
reg trunk foreign gear_ratio if rep78 == 4    // Outcome = Trunk and sub-group = 5
estimates store trunk_s4  

reg mpg foreign gear_ratio if rep78 ==5  // Outcome = MPG and sub-group = 4
estimates store mpg_s5
reg mpg foreign gear_ratio if rep78 ==4  // Outcome = MPG and sub-group = 5
estimates store mpg_s4  

reg turn foreign gear_ratio if rep78 ==5 // Outcome = Turn and sub-group = 4
estimates store turn_s5
reg turn foreign gear_ratio if rep78 ==4  // Outcome = Turn and sub-group = 5
estimates store turn_s4  

coefplot (trunk_s4, keep(foreign) \ mpg_s4 \ turn_s4, keep(gear_ratio) \) ///      
|| (trunk_s5, keep(foreign) \ mpg_s5 \ turn_s5, keep(gear_ratio) \) ///      
|| , drop(_cons) aseq swapnames eqrename(*_s4 = "" *_s5 = "") ///          
order(trunk turn mpg) bylabels("Subgroup 4" "Supgroup 5") ///            
eqlabels(, asheading)

I get the following plot:


What I would REALLY like to have happen is for the two trunk variables (trunk_s4 and trunk_s5) to be on the same line (i.e., side-by-side), AND to be able to label that as "Trunk". And the same goes for MPG, Turn, etc. Is there a way to do that?? I have tried everything I can think of, and have been at this for several days now.

[Hemanshu Kumar]

label var trunk "Trunk"
label var mpg "MPG"
label var turn "Turn"
#delimit ;
    coefplot (trunk_s4 \ mpg_s4 \ turn_s4, keep(foreign gear_ratio) )  
            || (trunk_s5 \ mpg_s5 \ turn_s5, keep(foreign  gear_ratio) ) ,
            eqrename(*_s4 = "" *_s5 = "")
            swapnames aseq eqlabels("Origin" "Gear Ratio")
            bylabels("Subgroup 4" "Subgroup 5")
            scheme(s2color)
                ;
#delimit cr

which produces:

[amandacpac]

So, that plot looks like exactly what I need, but when I run your code, I don't get that. Here is the code I am using:

sysuse auto, clear  
reg trunk foreign gear_ratio if rep78 == 5  // Outcome = Trunk and sub-group = 4
estimates store trunk_s5
reg trunk foreign gear_ratio if rep78 == 4    // Outcome = Trunk and sub-group = 5
estimates store trunk_s4  

reg mpg foreign gear_ratio if rep78 ==5  // Outcome = MPG and sub-group = 4
estimates store mpg_s5
reg mpg foreign gear_ratio if rep78 ==4  // Outcome = MPG and sub-group = 5
estimates store mpg_s4  

reg turn foreign gear_ratio if rep78 ==5 // Outcome = Turn and sub-group = 4
estimates store turn_s5
reg turn foreign gear_ratio if rep78 ==4  // Outcome = Turn and sub-group = 5
estimates store turn_s4  

label var trunk "Trunk"
label var mpg "MPG"
label var turn "Turn"

#delimit ;
    coefplot (trunk_s4 \ mpg_s4 \ turn_s4, keep(foreign gear_ratio) )  
            || (trunk_s5 \ mpg_s5 \ turn_s5, keep(foreign  gear_ratio) ) ,
            eqrename(*_s4 = "" *_s5 = "")
            swapnames aseq eqlabels("Origin" "Gear Ratio")
            bylabels("Subgroup 4" "Subgroup 5")
            scheme(s2color)
                ;
#delimit cr

And here is the plot I get:

[Hemanshu Kumar]

That's peculiar, I just copy-pasted your code into my Stata, and got what I showed you, not what you get. Perhaps this is a version issue?

If you do

which coefplot

do you get version 1.8.5?

If not, you might want to

ssc install coefplot, replace

or even do (to update all your ado files at one shot)

ado update, update

Also, for what it's worth, I'm on Stata 17/MP on a Mac.

[amandacpac]

That fixed it. Though, I had literally just updated Stata like 2 days ago. But, for some reason, I still had version 1.7.4. Thank you so much for your help!

[Hemanshu Kumar]

Great!

Updating Stata by itself does not update community-contributed commands you may have installed, including coefplot. For that, one needs to run the ado update command. Hence the problem.

[amandacpac]

I have a similar bit of code that I am using for a different dataset, but I am trying to get sample sizes listed on the right y-axis (each model has different sample sizes due to missing values and different time periods). I have code to pull the Ns, but can't figure out how to get it displayed. I am using the ymlab option with various local macros. But, for ease with this example, we can use simple text. So, for example, How can I modify the following code to display the strings: first, second, third, etc on the right y-axis, but differently for each plot. I hope that makes sense. I have attached a manually modified image to kind of show what I'm trying to do with the code below.

sysuse auto, clear  
reg trunk foreign gear_ratio if rep78 == 5  // Outcome = Trunk and sub-group = 4
estimates store trunk_s5
reg trunk foreign gear_ratio if rep78 == 4    // Outcome = Trunk and sub-group = 5
estimates store trunk_s4  

reg mpg foreign gear_ratio if rep78 ==5  // Outcome = MPG and sub-group = 4
estimates store mpg_s5
reg mpg foreign gear_ratio if rep78 ==4  // Outcome = MPG and sub-group = 5
estimates store mpg_s4  

reg turn foreign gear_ratio if rep78 ==5 // Outcome = Turn and sub-group = 4
estimates store turn_s5
reg turn foreign gear_ratio if rep78 ==4  // Outcome = Turn and sub-group = 5
estimates store turn_s4  

label var trunk "Trunk"
label var mpg "MPG"
label var turn "Turn"

#delimit ;
    coefplot (trunk_s4 \ mpg_s4 \ turn_s4, keep(foreign gear_ratio)), ymlab(1 "First" 2 "Second" 3 "Third", axis(2)) ||  
            (trunk_s5 \ mpg_s5 \ turn_s5, keep(foreign  gear_ratio)), ymlab(1 "Fourth" 2 "Fifth" 3 "Sixth", axis(2)) || ,
            eqrename(*_s4 = "" *_s5 = "")
            swapnames aseq eqlabels("Origin" "Gear Ratio")
            bylabels("Subgroup 4" "Subgroup 5") nooffset
    ;
#delimit cr