Dear Statalisters,

I would like to get the following table for my summarized statistics. However, I don't know the right code to get the right column with p-value for difference.

Array

I'm using the following code so far:

Code:
by lobby: eststo: quietly estpost tabstat lev-ceocob , stat(mean median sd) col(stat)
eststo: quietly estpost ttest lev-ceocob, by(lobby)
esttab, cells ("mean(fmt(3)) p50(fmt(3)) sd(fmt(3))")
esttab using Tabel2.rtf, noobs cells("mu_2 mu_1 b(star)") label replace
The variable lobby indicates the type of firms with a binary variable, whether it is a Lobbying or Non-Lobbying Firm.

Is there anyone who could help me? Thanks in advance!