标准化系数该如何利用coefplot命令作图？

coefplot命令可以将回归系数及置信区间绘制在图中，更直观地展示回归结果。但同时为了方便比较系数之间的大小，文献中更多地展示了标准化之后的系数图，如下。如果回归是ols，则使用coefplot命令展示标准化的系数图非常方便，只需要在reg后加入beta选项即可。但如果是其他回归，则回归命令不能直接支持beta选项。
有个帖子提供了回归系数标准化的很好的方法，见https://bbs.pinggu.org/thread-10674527-1-1.html
但如何根据标准化之后的系数利用coefplot来作图呢？根据我的尝试，coefplot命令，est sto命令均不支持beta选项。
想请教大家，我该如何将标准化之后的系数保存并作图呢？

葱葱饼干 发表于 2023-4-11 14:51
特别感谢您的回答！还想请教一下，如果想把多个回归的标准化系数展示在同一个图中，用这种矩阵的方法该怎 ...

这样应该比较接近你提供的论文的形式了
webuse grunfeld,clear
reghdfe invest mvalue kstock,a(company) cl(company)
mat V1 = r(table)[2,1...]
esttab ,beta
mat B1 =e(beta)

reghdfe invest mvalue kstock,a(company time) cl(company)
mat V2 = r(table)[2,1...]
esttab ,beta
mat B2 =e(beta)
coefplot (matrix(B1), se(V1) drop(_cons) offset(0.05) label(Model A)) (matrix(B2), se(V2) drop(_cons) offset(-0.05) label(Model A)), ciopts(recast(rcap))

可以先将变量标准化之后再回归就可以避免这个问题

还有一种方法是按照你展示的那个帖子的做法，先回归，然后用esttab,beta得到标准化的系数，最后通过矩阵传递给coefplot

1. webuse grunfeld,clear
   
2. reghdfe invest mvalue kstock,a(company) cl(company)
   
3. mat V = r(table)[2,1...]
   
4. esttab ,beta
   
5. mat B =  e(beta)
   
6. coefplot matrix(B), se(V) drop(_cons)

复制代码

917968079 发表于 2023-4-11 10:42
还有一种方法是按照你展示的那个帖子的做法，先回归，然后用esttab,beta得到标准化的系数，最后通过矩阵传 ...

特别感谢您的回答！还想请教一下，如果想把多个回归的标准化系数展示在同一个图中，用这种矩阵的方法该怎么做呢？

917968079 发表于 2023-4-11 10:13
可以先将变量标准化之后再回归就可以避免这个问题

谢谢您的解答~但是我使用变量标准化来回归的结果，和先回归再求其标准化系数不一致呢？命令如下：
//变量先标准化处理，再进行回归：
norm y x1 x2 x3 ,method(mmx)
xtreg  mmx_y mmx_x1 mmx_x2 mmx_x3 ,fe vce (cluster xxx)
//直接进行回归，并求标准化系数：
xtreg  y x1 x2 x3 ,fe vce (cluster xxx)
esttab ,beta

这两个方法求得的系数并不一致。请问我是哪里出错了呢？

葱葱饼干 发表于 2023-4-11 14:51
特别感谢您的回答！还想请教一下，如果想把多个回归的标准化系数展示在同一个图中，用这种矩阵的方法该怎 ...

以此类推

1. webuse grunfeld,clear
   
2. reghdfe invest mvalue kstock,a(company) cl(company)
   
3. mat V1 = r(table)[2,1...]
   
4. esttab ,beta
   
5. mat B1 =e(beta)
   
6. 
   
7. reghdfe invest mvalue kstock,a(company time) cl(company)
   
8. mat V2 = r(table)[2,1...]
   
9. esttab ,beta
   
10. mat B2 =e(beta)
    
11. coefplot (matrix(B1), se(V1) drop(_cons)) (matrix(B2), se(V2) drop(_cons))

复制代码

葱葱饼干 发表于 2023-4-11 15:17
谢谢您的解答~但是我使用变量标准化来回归的结果，和先回归再求其标准化系数不一致呢？命令如下：
//变量 ...

norm method要选zee

1. norm y x1 x2 x3 ,method(zee)

复制代码

917968079 发表于 2023-4-10 23:40
这样应该比较接近你提供的论文的形式了
webuse grunfeld,clear
reghdfe invest mvalue kstock,a(compan ...

非常棒，再次感谢！