Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, is already up and running.

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

Re: st: Deriving Bayes estimates from xtmelogit

From   Jamie Fagg <>
Subject   Re: st: Deriving Bayes estimates from xtmelogit
Date   Fri, 10 Sep 2010 09:15:13 +0100

Dear Bobby and Stas,

Thanks for your replies, they were both very helpful.

Best wishes,


On 9 September 2010 17:11, Roberto G. Gutierrez, StataCorp
<> wrote:
> Jamie Fagg <> asks:
>> I am trying to derive Bayes estimates from a 3-level logistic regression
>> model in Stata version 10.1.
>> The model has the following structure, where psychometric items are nested
>> in individuals, nested in geographic areas.
>> Level 1: psychometric items
>> Level 2: individuals
>> Level 3: geographic areas
>> I read in Rabe-Hesketh and Skrondal (Multilevel and Longitudinal modelling
>> using Stata, 2008, p.162) that "Empirical Bayes predictions of the random
>> intercepts ... can be obtained" from xtmixed using -predict- with the
>> reffects option
>> I couldn't find any analogous advice in the section about multilevel
>> logistic regression models so I ran the model and ran the predict eb,
>> reffects (see code below).
>> Is this the correct way to derive Bayes predictions from xtmelogit?
> Yes this would be correct, with the one caveat that what you obtain are
> empirical Bayes _modal_ predictions, rather than empirical Bayes mean
> predictions.  In a logistic regression setting, the posterior distribution of
> the random effects is no longer symmetric, and thus the posterior modes and
> the posterior means, while very similar for most data, are not strictly equal.
> Posterior mean predictions after multilevel logistic regression are not
> available in current official Stata, but you could obtain these by
> alternatively fitting your model using -gllamm-; see -ssc describe gllamm- for
> details.
> --Bobby
> *
> *   For searches and help try:
> *
> *
> *

Room 2.09, City Centre
Dept. of Geography, Queen Mary, University of London
Mile End Rd
E1 4NS
Tel: 020 7882 2748

*   For searches and help try:

© Copyright 1996–2016 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index