想请教一下大佬们，为什么这个分位数回归的命令，每次跑出来的回归系数结果均不同，是什么原因造成的呢？应该如何解决？  具体代码如下：
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   

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