2017-02-10 2 views
0

Ich möchte mit XTREG die zufälligen Effekte Abschnitte für einzelne Gruppen und ihre vorhergesagten Werte erhalten. Alle Vorhersagebefehle füllen jedoch nur alle Gruppen mit dem konstanten Wert.Stata mit XTREG für Cluster zufällige Effekte Modelle

Mit 1 Jahr Daten:

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

Antwort

1

Dies ist zum Teil eine Frage Statistiken. Der Kontext macht hier keinen Unterschied zu einem Standardmerkmal der Regression: Wenn keine Prädiktoren angegeben sind, dann ist die zurückgegebene Vorhersage notwendigerweise konstant und gleich der mittleren Antwort. Wir können Ihre Daten nicht sehen, aber das Prinzip ist universell. Hier ist ein reproduzierbares Beispiel:

. 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 

Beachten Sie, dass die Namen, die Sie verwenden (Beispiele in Ihrem Fall xb und u) haben keine Auswirkungen für das, was von predict produziert wird. Der Standardwert von predict ist in beiden Fällen ohne prognostizierte Prädiktoren das prognostizierte mittlere Ergebnis.

Sie sollten die Optionen von predict in diesem Fall wie

predict u, u 
+0

Danke für die Antwort überprüfen. Das u ist in jeder Gruppe gleich (variiert über die Gruppen hinweg), aber xb wird immer noch mit der Gesamtkonstanten für alle Beobachtungen gefüllt. Ich möchte die einzelnen zufälligen Effekte für jede Gruppe bekommen - sagst du, dass dies statistisch nicht möglich ist? Gibt es auch eine Möglichkeit, Modelle für jedes Jahr einzeln zu betreiben? Ich bekomme einen unzureichenden Beobachtungsfehler. – username

+0

Wenn Sie mein Beispiel überprüfen, unterscheidet sich "u" je nach Panel. Wir können Ihr Beispiel nicht überprüfen. Ich habe bereits erklärt, dass 'xb' für das, was du getan hast, notwendigerweise konstant ist. Sie stellen auch eine neue Frage in einem Kommentar; Sie sollten eine neue Frage stellen, Ihren Code zeigen und ein reproduzierbares Beispiel geben. Bitte lesen Sie http://stackoverflow.com/help/mcve –

Verwandte Themen