Stata 由于错误消息，无法计算sqrt函数

尝试在Stata 15.0上基于简单分位数回归进行单侧Wald检验

qreg volume if period == 2010 

---

我想测试我的

qreg的系数（常数）是否要执行单边测试，您需要先执行相应的双面Wald测试。然后可以使用结果计算单侧检验的检验统计量和p值

您的代码不工作，因为您没有执行Wald测试，所以没有定义r（df_r）和r（F）。当你试图取一个尚未计算的量的平方根时，Stata会给你一个神秘的错误。还有其他问题，比如在符号计算中使用1而不是10

我认为，在计算出问题之前，分解问题并显示所有部分是解决问题的好方法

下面是一个在cars数据集上显示这一点的示例（使用
dataex
的替代方法）：

。系统使用自动，清除
（1978年汽车数据）
. qreg价格重量长度i.国外，nolog
obs的回归中位数=74
偏差的原始总和71102.5（约4934）
最小偏差总和54411.29伪R2=0.2347
------------------------------------------------------------------------------
价格系数。标准错误。t P>| t |[95%配置间隔]
-------------+----------------------------------------------------------------
重量| 3.933588 1.328718 2.96 0.004 1.283543 6.583632
长度|-41.2519145.46469-0.910.367-131.928449.42456
|
外国的|
国外| 3377.771 885.4198 3.81 0.000 1611.857 5143.685
_cons | 344.6489 5182.394 0.07 0.947-9991.31 10680.61
------------------------------------------------------------------------------
. 本地符号=符号（_b[_cons]-10）
. 显示“`签名'”
1.
. 显示r（df_r）
.
. 显示器r（F）
.
. 测试b[_cons]=10
（1）_cons=10
F（1，70）=0.00
概率>F=0.9487
. 显示r（df_r）
70
. 显示器r（F）
.00416983
. 显示“Ho:_cons请向我们提供一个using Stata的
dataex
命令，以便我们可以尝试帮助您。没有样本数据，我们无法复制任何内容。我有超过一百万个观察结果，那么它是如何工作的？
帮助dataex
了解！
local to_test = sign(_b[_cons]-1 )
display "H0: volume <=10 p-value=" ttail(r(df_r),`to_test' * sqrt(r(F)))

H0: volume <=10 p-value=unknown function *sqrt()
r(133);

. sysuse auto, clear
(1978 Automobile Data)

. qreg price weight length i.foreign, nolog

Median regression                                   Number of obs =         74
  Raw sum of deviations  71102.5 (about 4934)
  Min sum of deviations 54411.29                    Pseudo R2     =     0.2347

------------------------------------------------------------------------------
       price |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      weight |   3.933588   1.328718     2.96   0.004     1.283543    6.583632
      length |  -41.25191   45.46469    -0.91   0.367    -131.9284    49.42456
             |
     foreign |
    Foreign  |   3377.771   885.4198     3.81   0.000     1611.857    5143.685
       _cons |   344.6489   5182.394     0.07   0.947     -9991.31    10680.61
------------------------------------------------------------------------------

. local sign_cons = sign(_b[_cons] - 10)

. display "`sign_cons'"
1

. display r(df_r)
.

. display r(F)
.

. test _b[_cons] = 10

 ( 1)  _cons = 10

       F(  1,    70) =    0.00
            Prob > F =    0.9487

. display r(df_r)
70

. display r(F)
.00416983

. display "Ho: _cons <= 10 p-value = " ttail(r(df_r),`sign_cons'*sqrt(r(F)))
Ho: _cons <= 10 p-value = .47434855