如何解决利润率与lincom / nlcom在Stata中有什么不同?
Stata中的margins
和lincom
之间到底有什么区别?如何调整nlcom
中的手动公式以匹配margins
的结果?谢谢
* load data
use http://www.stata-press.com/data/r13/nlswork
* set panel structure
xtset idcode year
* fixed effects regression
xtreg ln_wage c.wks_ue##c.wks_ue union age,fe coeflegend
margins,dydx(wks_ue)
------------------------------------------------------------------------------
| Delta-method
| dy/dx Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
wks_ue | -.00594 .0009744 -6.10 0.000 -.0078497 -.0040303
------------------------------------------------------------------------------
* check with lincom/nlcom
lincom _b[wks_ue] + 2*_b[c.wks_ue#c.wks_ue]
nlcom _b[wks_ue] + 2*_b[c.wks_ue#c.wks_ue]
------------------------------------------------------------------------------
ln_wage | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_nl_1 | -.0062406 .0010319 -6.05 0.000 -.0082632 -.0042181
------------------------------------------------------------------------------
解决方法
我错过了wks_ue
的平均值。从估计的样本中获取它,将其另存为本地,并将其包含在lincom
中。系数和标准误差将相同。 lincom
使用t分布,margins
使用z分布。
* load data
use http://www.stata-press.com/data/r13/nlswork
* set panel structure
xtset idcode year
* fixed effects regression
xtreg ln_wage c.wks_ue##c.wks_ue union age,fe coeflegend
margins,dydx(wks_ue)
qui sum wks_ue if e(sample)
local wks_ue_mean = r(mean)
lincom wks_ue + 2*c.wks_ue#c.wks_ue*`wks_ue_mean'
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。