Common Title for two columns and the mean of dependent variable in a regression table

[Olivia Emma]

Hi I have to present the results of following regression commands in a table form

I am unable to add the Common Title for two columns and the mean of dependent variable.

xi: ivreg2 y Z x Z_x , r cluster(sample_cluster)
estimates store column1
estadd local Controls "yes"
estadd local Region "yes"
estadd local Survey "yes"

*** Column 2


xi: ivreg2 y Z x Z_x i.survey i.region i.urban n_siblings , r cluster(sample_cluster)
estimates store column2
estadd local Controls "yes"
estadd local Region "yes"
estadd local Survey "yes"
*** Column 3


xi: ivreg2 y Z x Z_x , r cluster(sample_cluster)
estimates store column3
estadd local Controls "yes"
estadd local Region "yes"
estadd local Survey "yes"

*** Column 4


xi: ivreg2 y Z x Z_x i.survey i.region i.urban n_siblings , r cluster(sample_cluster)
estimates store column4
estadd local Controls "yes"
estadd local Region "yes"
estadd local Survey "yes"


*** Column 5

xi: ivreg2 y Z x Z_x , r cluster(sample_cluster)
estimates store column5
estadd local Controls "yes"
estadd local Region "yes"
estadd local Survey "yes"

*** Column 6


xi: ivreg2 y Z x Z_x i.survey i.region i.urban n_siblings , r cluster(sample_cluster)
estimates store column6
estadd local Controls "yes"
estadd local Region "yes"
estadd local Survey "yes"

how I can add the mean of dependent variable and common Title for two column please

esttab column1 column2 column3 column4 column5 column6 using Table1.rtf , replace s(N Controls Region Survey, label("N" "Controls " "Region")) ar2 se noconstant star(* .10 ** .05 *** .01) keep(Z x ) varlabels(Z "Control 1" x "Control 2" )