Question

我想使用xtreg来获取各个组的随机效果拦截及其预测值。但是，所有预测命令都只是使用常量值填充所有组。

使用1年的数据：

xtset group // set panel
xtreg outcome, re
predict u
predict xb

Answer 1

这部分是一个统计问题。这里的xt上下文对回归的标准特征没有影响：如果没有指定预测变量，则返回的预测必然是常数并且等于平均响应。我们看不到您的数据，但原则是普遍的。这是一个可重复的例子：

. webuse grunfeld

. xtset
       panel variable:  company (strongly balanced)
        time variable:  year, 1935 to 1954
                delta:  1 year

. xtreg mvalue

Random-effects GLS regression                   Number of obs     =        200
Group variable: company                         Number of groups  =         10

R-sq:                                           Obs per group:
     within  = 0.0000                                         min =         20
     between = 0.0000                                         avg =       20.0
     overall = 0.0000                                         max =         20

                                                Wald chi2(0)      =          .
corr(u_i, X)   = 0 (assumed)                    Prob > chi2       =          .

------------------------------------------------------------------------------
      mvalue |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
       _cons |   1081.681   422.1377     2.56   0.010     254.3064    1909.056
-------------+----------------------------------------------------------------
     sigma_u |  1332.6401
     sigma_e |  348.51426
         rho |  .93598465   (fraction of variance due to u_i)
------------------------------------------------------------------------------

. predict predict
(option xb assumed; fitted values)

. su predict

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     predict |        200    1081.681           0   1081.681   1081.681

. su mvalue

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
      mvalue |        200    1081.681     1314.47      58.12     6241.7

请注意，您使用的名称（案例xb和u中的示例）对predict生成的内容没有任何影响。在任何一种情况下，没有指定预测因子的默认值predict是预测的平均结果。

在这种情况下，您应该查看predict的选项，例如

predict u, u

Stata使用xtreg进行群集随机效应模型

1 个答案: