Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: [Fwd: Hierarchical model]


From   Evans Jadotte <evans.jadotte@uab.es>
To   statalist@hsphsun2.harvard.edu
Subject   st: [Fwd: Hierarchical model]
Date   Fri, 02 Oct 2009 18:26:37 +0200

Hello listers!

I find it hard to believe that no one has never used hierarchical model out there. The last 3 inquiries I made here were left with ZERO feedback. May be something is wrong with my email. any way, if somebody has any idea I would so much appreciate some feedback from the thread in this mail.

Thanks!

EJ
--- Begin Message ---
From   Evans Jadotte <evans.jadotte@uab.es>
To   statalist@hsphsun2.harvard.edu
Subject   Hierarchical model
Date   Fri, 02 Oct 2009 11:26:36 +0200
Dear statalisters,

I am estimating a linear three-level hierarchical model via "restricted maximum likelihood-REML" with the /xtmixed/ command, with households (level-1) nested into villages (level-2) nested into regions (level-3). Specification seems alright. However, the empirical Bayes random intercept at level-3, which according to theory should be inferior or equal to the corresponding calculated REML estimate, do not reflect the REML estimate that I calculated at that level. Specifically, I have an empirical Bayes (the BLUP r.e. computed by Stata) that is more than 300 times the REML estimate calculated at level-3. At level-2 it seems reasonable, empirical Bayes = 0.87 REML estimate.

Moreover, I expected the sum of the variances (that of raw residuals at level-1 + that of residuals REML level-2 + that of residuals REML level-3). When I compute them however, the result gives an estimated total variance that is lower than that of the raw residuals at level-1.

I am pretty confident about the model specification and the computation but evidently things do not match as they should. Can anyone familiar with these models give me any feedback? I would so much appreciate and thanks in advance.

Evans


--- End Message ---



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index