MCMC Methods for Multi-Response Generalized Linear Mixed Models: The MCMCglmm R Package
TL;DR: The R package MCMCglmm implements Markov chain Monte Carlo methods for generalized linear mixed models, which provide a flexible framework for modeling a range of data, although with non-Gaussian response variables the likelihood cannot be obtained in closed form.
read more
Abstract: Generalized 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.
read more
Chat with Paper
AI Agents for this Paper
Find similar papers on Google Scholar, PubMed and Arxiv
Write a critical review of this paper
Analyze citations of this paper to find unaddressed research gaps
Citations
A general and simple method for obtaining R2 from generalized linear mixed-effects models
TL;DR: In this article, the authors make a case for the importance of reporting variance explained (R2) as a relevant summarizing statistic of mixed-effects models, which is rare, even though R2 is routinely reported for linear models and also generalized linear models (GLM).
9.5K
Random effects structure for confirmatory hypothesis testing: Keep it maximal
TL;DR: It is argued that researchers using LMEMs for confirmatory hypothesis testing should minimally adhere to the standards that have been in place for many decades, and it is shown thatLMEMs generalize best when they include the maximal random effects structure justified by the design.
9.3K
glmmTMB Balances Speed and Flexibility Among Packages for Zero-inflated Generalized Linear Mixed Modeling
Mollie Elizabeth Brooks,Kasper Kristensen,Koen J. van Benthem,Arni Magnusson,Casper Willestofte Berg,Anders Nielsen,Hans J. Skaug,Martin Mächler,Benjamin M. Bolker +8 more
TL;DR: The glmmTMB package fits many types of GLMMs and extensions, including models with continuously distributed responses, but here the authors focus on count responses and its ability to estimate the Conway-Maxwell-Poisson distribution parameterized by the mean is unique.
brms: An R Package for Bayesian Multilevel Models Using Stan
TL;DR: The brms package implements Bayesian multilevel models in R using the probabilistic programming language Stan, allowing users to fit linear, robust linear, binomial, Poisson, survival, ordinal, zero-inflated, hurdle, and even non-linear models all in a multileVEL context.
7.4K
Advanced Bayesian Multilevel Modeling with the R Package brms
TL;DR: Brms provides an intuitive and powerful formula syntax, which extends the well known formula syntax of lme4, which is introduced in detail and demonstrated its usefulness with four examples, each showing other relevant aspects of the syntax.
References
Bayesian measures of model complexity and fit
TL;DR: In this paper, the authors consider the problem of comparing complex hierarchical models in which the number of parameters is not clearly defined and derive a measure pD for the effective number in a model as the difference between the posterior mean of the deviances and the deviance at the posterior means of the parameters of interest, which is related to other information criteria and has an approximate decision theoretic justification.
•Journal Article
Bayesian measures of model complexity and fit
TL;DR: The posterior mean deviance is suggested as a Bayesian measure of fit or adequacy, and the contributions of individual observations to the fit and complexity can give rise to a diagnostic plot of deviance residuals against leverages.
7.6K
R Development Core Team (2010): R: A language and environment for statistical computing
Peter Dalgaard
- 01 Jan 2010
TL;DR: In this article, the R Foundation for Statistical Computing (RFC) gave permission to make and distribute verbatim copies of this manual provided the copyright notice and this permission notice are preserved on all copies.
7.3K
JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling
Martyn Plummer
- 01 Jan 2003
TL;DR: JAGS is a program for Bayesian Graphical modelling which aims for compatibility with Classic BUGS and could eventually be developed as an R package.
Approximate inference in generalized linear mixed models
Norman E. Breslow,D. G. Clayton +1 more
TL;DR: In this paper, generalized linear mixed models (GLMM) are used to estimate the marginal quasi-likelihood for the mean parameters and the conditional variance for the variances, and the dispersion matrix is specified in terms of a rank deficient inverse covariance matrix.
4.6K