MCMC Methods for Multi-Response Generalized Linear Mixed Models: The MCMCglmm R Package
AbstractGeneralized linear mixed models provide a flexible framework for modeling a range of data, although with non-Gaussian response variables the likelihood cannot be obtained in closed form. Markov chain Monte Carlo methods solve this problem by sampling from a series of simpler conditional distributions that can be evaluated. The R package MCMCglmm implements such an algorithm for a range of model fitting problems. More than one response variable can be analyzed simultaneously, and these variables are allowed to follow Gaussian, Poisson, multi(bi)nominal, exponential, zero-inflated and censored distributions. A range of variance structures are permitted for the random effects, including interactions with categorical or continuous variables (i.e., random regression), and more complicated variance structures that arise through shared ancestry, either through a pedigree or through a phylogeny. Missing values are permitted in the response variable(s) and data can be known up to some level of measurement error as in meta-analysis. All simu- lation is done in C/ C++ using the CSparse library for sparse linear systems.
Download InfoIf you experience problems downloading a file, check if you have the proper application to view it first. In case of further problems read the IDEAS help page. Note that these files are not on the IDEAS site. Please be patient as the files may be large.
Bibliographic InfoArticle provided by American Statistical Association in its journal Journal of Statistical Software.
Volume (Year): 33 ()
Issue (Month): i02 ()
Contact details of provider:
Web page: http://www.jstatsoft.org/
Please report citation or reference errors to , or , if you are the registered author of the cited work, log in to your RePEc Author Service profile, click on "citations" and make appropriate adjustments.:
- P. Damlen & J. Wakefield & S. Walker, 1999. "Gibbs sampling for Bayesian non-conjugate and hierarchical models by using auxiliary variables," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 61(2), pages 331-344.
- David J. Spiegelhalter & Nicola G. Best & Bradley P. Carlin & Angelika van der Linde, 2002. "Bayesian measures of model complexity and fit," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 64(4), pages 583-639.
- William J. Browne & Fiona Steele & Mousa Golalizadeh & Martin J. Green, 2009. "The use of simple reparameterizations to improve the efficiency of Markov chain Monte Carlo estimation for multilevel models with applications to discrete time survival models," Journal of the Royal Statistical Society Series A, Royal Statistical Society, vol. 172(3), pages 579-598.
- Irene Albarrán & J. Miguel Marín & Pablo J. Alonso, 2011. "Why using a general model in Solvency II is not a good idea : an explanation from a Bayesian point of view," Statistics and Econometrics Working Papers ws113729, Universidad Carlos III, Departamento de Estadística y Econometría.
- Jesse Shore & Ethan Bernstein & David Lazer, 2014. "Facts and Figuring: An Experimental Investigation of Network Structure and Performance in Information and Solution Spaces," Harvard Business School Working Papers 14-075, Harvard Business School, revised Jun 2014.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Christopher F. Baum).
If references are entirely missing, you can add them using this form.