2017-02-10 577 views
0

我想用xtreg获得随机效应拦截个别团体和他们的预测值。然而,所有的预测命令只是填充所有组的常数值。Stata的使用xtreg集群随机效应模型

利用1年数据:

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

回答

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 

注意,名称使用(在你的情况xbu例子),对于由predict产生任何影响。的predict在这两种情况下指定的默认,没有预测变量,是预测均值的结果。

你应该检查出的predict的选项在这种情况下,如

predict u, u 
+0

感谢您的答复。的u是各组(各组间不同而不同)中的相同,但XB仍在与在所有的观察整体恒定填充。我想获得各组的个体随机效应 - 你说,这是不可能的统计?是否还有一种方法可以单独运行每年的模型?我收到的观察错误不足。 – username

+0

如果你检查我的例子,那么'u'因面板而异。我们无法检查你的例子。我已经解释过,'xb'对于你所做的事情来说必然是不变的。你也在评论中提出一个新问题;你应该发布一个新问题,展示你的代码并给出一个可重现的例子。请查看http://stackoverflow.com/help/mcve –