为何面板分位回归每次得到的结果不同？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

用的回归命令是qregpd
直接采用stata中help文件中给的例子
webuse nlswork
qregpd ln_wage tenure union, id(idcode) fix(year) optimize(mcmc) noisy draws(1000) burn(100) arate(.5)
mat list e(gamma)
结果每次重复执行这个回归指令都会得到不同的结果，请教大神这是怎么回事呢？应该如何处理呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：分位回归 结果不同 OPTIMIZE idcode Tenure

加一句  set seed   123456789
就可以了

1. webuse nlswork
   
2. set seed 123456
   
3. qregpd ln_wage tenure union, id(idcode) fix(year) optimize(mcmc) noisy  draws(1000) burn(100) arate(.5)
   
4. mat list e(gamma)

复制代码

蓝色 发表于 2018-1-1 13:38
加一句  set seed   123456789
就可以了

可以了，太感谢您了！

请问大家用qregpd做面板分位回归时，跑qregpd ln_wage tenure union, id(idcode) fix(year) optimize(mcmc) noisy draws(1000) burn(100) arate(.5)这个回归指令的结果有报告截距项吗？还是这种处理方法的结果本来就没有截距项？求了解的坛友解答一下，万分感谢！

蓝色 发表于 2018-1-1 13:38
加一句  set seed   123456789
就可以了

您好，请问设置任意seed，比如我set seed 1结果显著而seed 2结果不显著，那怎么说明seed 1的结果是可靠的呢？

chimny-1024 发表于 2018-1-1 12:10
用的回归命令是qregpd
直接采用stata中help文件中给的例子
webuse nlswork

同问～

Rita_P 发表于 2020-5-3 22:06
您好，请问设置任意seed，比如我set seed 1结果显著而seed 2结果不显著，那怎么说明seed 1的结果是可靠的 ...

不能因为自己想要结果是显著的，就认为显著结果是可靠的。
结果受随机数影响，只说明结果是不稳健的。
试试增加draws和burn的数量

heric221 发表于 2020-5-5 15:00
不能因为自己想要结果是显著的，就认为显著结果是可靠的。
结果受随机数影响，只说明结果是不稳健的。
...

增大这两个参数直到显著性不会改变吗？就算显著性不改变，每次系数也会改变，这个无所谓吗？谢谢！

楼主能解释一下mat list e(gamma) 什么意思吗，最近也在用这个命令跑实证，是不是只要变量通过显著性水平检验就可以了，还需要做什么检验吗？

凌晨两点喝红茶 发表于 2021-8-11 11:26
楼主能解释一下mat list e(gamma) 什么意思吗，最近也在用这个命令跑实证，是不是只要变量通过显著性水平检 ...

mat list e(gamma)应该是显示year时点固定效应