Question

在运行Tobit回归时，是否有人知道让Stata显示p值超过三位数的简单方法？

通常，Stata报告p值为.001或.065，但我希望看到更多数字，例如.0011123或.065320。

要明确，我不希望（必然）改变回归表中数据的生成方式。

我只希望能够让Stata为我感兴趣的那些p值显示更多数字。

Answer 1

使用

跟进tobit命令

est tab, p(%12.10g)

（例如）。这应该适用于Stata的旧版本。不那么容易就是编写自己的输出程序。

Answer 2

Stata 11.1 introduced一个set pformat命令，用于指定系数表中p值的输出格式。（我不知道STATA我害怕，因为我认为在20世纪80年代的某个时候已经停止了。）

顺便说一下，你最好在Statalist而不是在这里询问完全与Stata相关的问题。

Answer 3

很多时候，如果你通过内部名称知道你的p值，你就可以获得最大的精确度。我经常在我将认真使用的几乎所有命令后输入return list或ereturn list，然后获取可能看起来像e(p)或r(p)或e(p_chi2)的结果或无论包含p值的标量可能是什么。

Answer 4

在进行tobit回归后，您可以使用test命令从零假设x1 = 0获取p值：

sysuse auto
tobit weight trunk length headroom, ll(1500)
test trunk

r(p)

中返回的结果

return list

Answer 5

使用tobit命令帮助文件中的第一个示例：

. sysuse auto, clear
. generate wgt = weight / 1000

. tobit mpg wgt, ll(17)

Tobit regression                                Number of obs     =         74
                                                LR chi2(1)        =      72.85
                                                Prob > chi2       =     0.0000
Log likelihood = -164.25438                     Pseudo R2         =     0.1815

------------------------------------------------------------------------------
         mpg |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         wgt |   -6.87305   .7002559    -9.82   0.000    -8.268658   -5.477442
       _cons |   41.49856    2.05838    20.16   0.000     37.39621     45.6009
-------------+----------------------------------------------------------------
      /sigma |   3.845701   .3663309                      3.115605    4.575797
------------------------------------------------------------------------------
            18  left-censored observations at mpg <= 17
            56     uncensored observations
             0 right-censored observations

您可以轻松地从r()中返回的结果中获取任何p值：

. matrix list r(table)

r(table)[9,3]
             model:      model:      sigma:
               wgt       _cons       _cons
     b  -6.8730504   41.498557   3.8457011
    se   .70025591   2.0583803   .36633085
     t  -9.8150552   20.160782          .b
pvalue   5.610e-15   1.471e-31          .b
    ll  -8.2686584   37.396211   3.1156048
    ul  -5.4774424   45.600903   4.5757975
    df          73          73          73
  crit   1.9929971   1.9929971   1.9929971
 eform           0           0           0

然后相应地格式化：

. matrix results = r(table)

. display %18.17f results[4,1]
0.00000000000000561

从Stata的命令提示符输入help format以获取更多信息。

如何获得p值的更多数字？

5 个答案: