0
我想用xtreg获得随机效应拦截个别团体和他们的预测值。然而,所有的预测命令只是填充所有组的常数值。Stata的使用xtreg集群随机效应模型
利用1年数据:
xtset group // set panel
xtreg outcome, re
predict u
predict xb
我想用xtreg获得随机效应拦截个别团体和他们的预测值。然而,所有的预测命令只是填充所有组的常数值。Stata的使用xtreg集群随机效应模型
利用1年数据:
xtset group // set panel
xtreg outcome, re
predict u
predict xb
这部分是一个统计问题。这里的上下文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
感谢您的答复。的u是各组(各组间不同而不同)中的相同,但XB仍在与在所有的观察整体恒定填充。我想获得各组的个体随机效应 - 你说,这是不可能的统计?是否还有一种方法可以单独运行每年的模型?我收到的观察错误不足。 – username
如果你检查我的例子,那么'u'因面板而异。我们无法检查你的例子。我已经解释过,'xb'对于你所做的事情来说必然是不变的。你也在评论中提出一个新问题;你应该发布一个新问题,展示你的代码并给出一个可重现的例子。请查看http://stackoverflow.com/help/mcve –