面板分位数回归 qregpd命令

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

想请教一下大佬们，为什么这个分位数回归的命令，每次跑出来的回归系数结果均不同，是什么原因造成的呢？应该如何解决？  具体代码如下：
clear
use "跨国数据全表.dta"

winsor2 y lnFS1_ lnfs2 lnKr import export fdi_net inovation Hi1 Found1 WGI
  xtset id year                                                                                                                                                                                                                                                                                                                                          
* 定义分位数
local quantiles 0.10 0.25 0.50 0.75 0.90
* 进行面板分位数回归
foreach q of local quantiles {
    qregpd y lnFS1_ lnfs2  lnKr import export fdi_net inovation Hi1 Found1 WGI, quantile(`q') id(id)  fix(year)  optimize(mcmc) noisy draws(1000) burn(100) arate(0.5)
  estimates store q`=int(`q'*100)'
}

esttab q10 q25 q50 q75 q90 using "面板分位数.rtf", replace   

同样的代码和数据跑两遍后的结果如下：
（是否是因为迭代次数导致的问题呢，还请各位大佬指点）

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：面板分位数回归 分位数回归 面板分位数 REG QRE 面板分位数回归； 分位数回归；

估计是里面缩尾的问题

可能是缩尾一次 和缩尾两次的差别，

wdlbcj 发表于 2024-10-21 13:59
估计是里面缩尾的问题

可能是缩尾一次 和缩尾两次的差别，

谢谢友友，昨天晚上鼓捣半天解决啦，是因为没有设置随机种子，蒙特卡洛模拟是随机抽样，如果种子不同就会导致结果不同。

暗香疏影_ 发表于 2024-10-21 15:11
谢谢友友，昨天晚上鼓捣半天解决啦，是因为没有设置随机种子，蒙特卡洛模拟是随机抽样，如果种子不同就会 ...

你好，请问解决了吗，代码应该需要修改哪里？

请问怎么修改呢？

温衡427 发表于 2025-1-20 20:58
你好，请问解决了吗，代码应该需要修改哪里？

local quantiles 0.10 0.25 0.50 0.75 0.90
set seed 100000

* 进行面板分位数回归，此处不设定随机种子跑出来的结果不稳定，每次跑一次结果都不一样
foreach q of local quantiles {
    qregpd lnGDP_per lnFS1_ lnfs2  lnKr import export fdi_net inovation Hi1 Found1 WGI, quantile(`q') id(id) fix(year) optimize(mcmc) noisy draws(1000) burn(100) arate(0.5)
  estimates store q`=int(`q'*100)'
}

local quantiles 0.10 0.25 0.50 0.75 0.90
set seed 100000

* 进行面板分位数回归，此处不设定随机种子跑出来的结果不稳定，每次跑一次结果都不一样
foreach q of local quantiles {
    qregpd lnGDP_per lnFS1_ lnfs2  lnKr import export fdi_net inovation Hi1 Found1 WGI, quantile(`q') id(id) fix(year) optimize(mcmc) noisy draws(1000) burn(100) arate(0.5)
  estimates store q`=int(`q'*100)'
}修改后的代码