Stata中最大似然模型辅助参数的提取

在Stata中，当使用最大似然（

ml

）估计时，我可以使用

diparm

（在

ado

文件中）来显示辅助（辅助）参数。但是，当我使用

ereturn list

时，我没有找到提取这些系数的方法。辅助参数是否有类似于材料清单e（b）的内容

ereturn list

scalars:
               e(rank) =  13
                  e(N) =  6378
                 e(ic) =  17
                  e(k) =  13
               e(k_eq) =  5
               e(k_dv) =  2
          e(converged) =  1
                 e(rc) =  0
          e(k_autoCns) =  0
                 e(ll) =  -7465.044032705097
         e(k_eq_model) =  1
               e(df_m) =  4
               e(chi2) =  65.60818581991281
                  e(p) =  1.91572317062e-13
              e(k_aux) =  2

macros:
                e(cmd) : "mypro"
           e(chi2type) : "Wald"
                e(opt) : "moptimize"
            e(predict) : "ml_p"
                e(vce) : "oim"
               e(user) : "mypro_ll"
           e(crittype) : "log likelihood"
          e(ml_method) : "lf"
    e(singularHmethod) : "hybrid"
          e(technique) : "nr"
              e(which) : "max"
             e(depvar) : "y1 y2"
         e(properties) : "b V"

matrices:
                  e(b) :  1 x 13
                  e(V) :  13 x 13
               e(ilog) :  1 x 20
           e(gradient) :  1 x 13
               e(ml_h) :  1 x 13
           e(ml_scale) :  1 x 13

functions:
             e(sample)   

---

辅助参数通常存储在

e（b）

矩阵中它们自己的“方程”中。见：

M.L.Buis（2011）“Stata技巧97：达到$\rho$s和$\sigma$s”，Stata期刊，11（2），第315-317页

---

而且：

非常感谢。这就是我真正想要的。在您的论文中，您还提到了另一篇关于测试边界值约束的论文的链接。我有一个参数，它的值介于0和1之间。如果我施加0或1的约束，Stata会产生“数字溢出错误代码1400”问题。但是，如果我施加0.0002或0.9999的约束，它会起作用。我想知道这是测试0还是1的可接受方法。Stata中的约束是线性等式；不平等被排除在外。重新参数化到罗吉特规模是你所寻求的。关于这一点有一个Stata FAQ。不幸的是，没有。这些测试通常假设该参数在零假设下为正态抽样分布。如果你在（或非常接近）边界上，那么（大约）一半的假设采样分布是不可能的。因此，通常，空值下正确的采样分布是某种混合分布。你应该使用哪个确切的发行版不是我的专长，所以你需要问其他人。我想你提到的另一篇文章是：古铁雷斯、R、S.卡特和D.M.德鲁克。2001关于边值似然比检验。Stata技术公报，60:15-18。请注意，它是免费提供的，因此无需猜测其内容：