predicted probabilities lmer
Well it depends on what exactly do you mean by 'FIXED part only'. If you mean the predicted probabilities for the "average" subject (i.e., a subject with random effect value zero), then you can simply obtain these probabilities from the estimated coefficients of the model. However, if you mean the marginal/population probabilities, i.e., the probabilities averaged over the subjects, then you in fact need, P(y = 1) = \int P(y = 1 | b) f(b) db, where f(b) is the normal density for the random effects. One option to get those is using sampling with a bit of extra programming. I hope it helps. Best, Dimitris
On 6/15/2011 2:47 PM, Ben Pelzer wrote:
Hi all, Is there a simple method to obtain the predicted probabilities based on the FIXED part only of a random intercept logistic regression model using lmer? Thanks for any help! Ben.
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
Dimitris Rizopoulos Assistant Professor Department of Biostatistics Erasmus University Medical Center Address: PO Box 2040, 3000 CA Rotterdam, the Netherlands Tel: +31/(0)10/7043478 Fax: +31/(0)10/7043014 Web: http://www.erasmusmc.nl/biostatistiek/