Fitting Linear Mixed-Effects Models Using lme4

Published on Oct 7, 2015in Journal of Statistical Software13.642
· DOI :10.18637/JSS.V067.I01
Douglas M. Bates38
Estimated H-index: 38
Martin Mächler13
Estimated H-index: 13
+ 1 AuthorsSteven C. Walker17
Estimated H-index: 17
Maximum likelihood or restricted maximum likelihood (REML) estimates of the parameters in linear mixed-effects models can be determined using the lmer function in the lme4 package for R. As for most model-fitting functions in R, the model is described in an lmer call by a formula, in this case including both fixed- and random-effects terms. The formula and data together determine a numerical representation of the model from which the profiled deviance or the profiled REML criterion can be evaluated as a function of some of the model parameters. The appropriate criterion is optimized, using one of the constrained optimization functions in R, to provide the parameter estimates. We describe the structure of the model, the steps in evaluating the profiled deviance or REML criterion, and the structure of classes or types that represents such a model. Sufficient detail is included to allow specialization of these structures by users who wish to write functions to fit specialized linear mixed models, such as models incorporating pedigrees or smoothing splines, that are not easily expressible in the formula language used by lmer.
📖 Papers frequently viewed together
14.5k Citations
5 Authors (Alain F. Zuur, ..., Graham M. Smith)
7,927 Citations
1 Author (Hadley Wickham)
16.8k Citations
#1C. G. KhatriH-Index: 1
involving an unknown function ?5 of a single variable t. Conditions under which the unknown functions in these two types of equations are polynomials of an assigned degree are given. The third, on the characterization of normal and gamma distributions, extends the earlier work of the authors (Rao, 1967 and Khatri and Rao, 1968*). We consider two sets of functions Li9 ..., Lq and Mlf ..., Mp of independent random variables Xl9 .... Xn with the condition
100 Citations
#1Douglas M. BatesH-Index: 38
#2Martin MaechlerH-Index: 10
Last. Steven WalkerH-Index: 2
view all 4 authors...
5,886 Citations
#1Kyle Klein (UCSB: University of California, Santa Barbara)H-Index: 7
#2Julian Neira (UCSB: University of California, Santa Barbara)H-Index: 5
The Nelder-Mead simplex method is an optimization routine that works well with irregular objective functions. For a function of $n parameters, it compares the objective function at the n+1 vertices of a simplex and updates the worst vertex through simplex search steps. However, a standard serial implementation can be prohibitively expensive for optimizations over a large number of parameters. We describe an implementation of the Nelder-Mead method in parallel using a distributed memory. F...
27 CitationsSource
Over the last twenty years there have been numerous developments in diagnostic pro- cedures for hierarchical linear models; however, these procedures are not widely imple- mented in statistical software packages, and those packages that do contain a complete framework for model assessment are not open source. The lack of availability of diagnostic procedures for hierarchical linear models has limited their adoption in statistical practice. The R package HLMdiag provides diagnostic tools targetin...
28 CitationsSource
124k Citations
#1Benjamin M. Bolker (McMaster University)H-Index: 62
#2Beth Gardner (Patuxent Wildlife Research Center)H-Index: 29
Last. Elise F. Zipkin (Patuxent Wildlife Research Center)H-Index: 27
view all 22 authors...
1. Ecologists often use nonlinear fitting techniques to estimate the parameters of complex ecological models, with attendant frustration. This paper compares three open-source model fitting tools and discusses general strategies for defining and fitting models. 2. R is convenient and (relatively) easy to learn, AD Model Builder is fast and robust but comes with a steep learning curve, while BUGS provides the greatest flexibility at the price of speed. 3. Our model-fitting suggestions range from ...
81 CitationsSource
#1Yeojin Chung (Kookmin University)H-Index: 4
#2Sophia Rabe-Hesketh (University of California, Berkeley)H-Index: 67
Last. Jingchen Liu (Columbia University)H-Index: 18
view all 5 authors...
Group-level variance estimates of zero often arise when fitting multilevel or hierarchical linear models, especially when the number of groups is small. For situations where zero variances are implausible a priori, we propose a maximum penalized likelihood approach to avoid such boundary estimates. This approach is equivalent to estimating variance parameters by their posterior mode, given a weakly informative prior distribution. By choosing the penalty from the log-gamma family with shape param...
199 CitationsSource
#1Jeff Bezanson (MIT: Massachusetts Institute of Technology)H-Index: 6
#2Stefan KarpinskiH-Index: 7
Last. Alan Edelman (MIT: Massachusetts Institute of Technology)H-Index: 43
view all 4 authors...
Dynamic languages have become popular for scientic computing. They are generally considered highly productive, but lacking in performance. This paper presents Julia, a new dynamic language for technical computing, designed for performance from the beginning by adapting and extending modern programming language techniques. A design based on generic functions and a rich type system simultaneously enables an expressive programming model and successful type inference, leading to good performance for...
394 Citations
#1Rense NieuwenhuisH-Index: 10
#2M. te GrotenhuisH-Index: 2
Last. Ben PelzerH-Index: 15
view all 3 authors...
influence.ME provides tools for de- tecting influential data in mixed effects mod- els. The application of these models has become common practice, but the development of diag- nostic tools has lagged behind. influence.ME calculates standardized measures of influential data for the point estimates of generalized mixed effects models, such as DFBETAS, Cook's dis- tance, as well as percentile change and a test for changing levels of significance. influence.ME calculates these measures of influence...
202 CitationsSource
#1John C. Nash (U of O: University of Ottawa)H-Index: 9
#2Ravi VaradhanH-Index: 44
R users can often solve optimization tasks easily using the tools in the optim function in the stats package provided by default on R installations. However, there are many other optimization and nonlinear modelling tools in R or in easily installed add-on packages. These present users with a bewildering array of choices. optimx is a wrapper to consolidate many of these choices for the optimization of functions that are mostly smooth with parameters at most bounds-constrained. We attempt to prov...
195 CitationsSource
Cited By30237
#1Charlotte Garcia (University of Nottingham)
#2Alexandra Țurcan (University of Nottingham)H-Index: 3
Last. Ruth Filik (University of Nottingham)H-Index: 23
view all 4 authors...
Abstract null null There is evidence for an age-related decline in the ability to understand non-literal language such as sarcasm. There is also evidence to suggest that devices such as emoticons/emojis may influence sarcasm comprehension in younger adults. However, research examining whether such devices may improve written sarcasm comprehension in older adults is scarce. The present study used an online rating task to investigate the influence of the winking face emoji on both the interpretati...
#1David W. Braithwaite (FSU: Florida State University)H-Index: 7
#2Jake McMullen (UTU: University of Turku)H-Index: 13
Last. Michelle Hurst (U of C: University of Chicago)H-Index: 5
view all 0 authors...
Understanding fractions and decimals requires not only understanding each notation separately, or within-notation knowledge, but also understanding relations between notations, or cross-notation knowledge. Multiple notations pose a challenge for learners but could also present an opportunity, in that cross-notation knowledge could help learners to achieve a better understanding of rational numbers than could easily be achieved from within-notation knowledge alone. This hypothesis was tested by r...
#1Bojan Luc Nys (Katholieke Universiteit Leuven)
#2Janie Brisson (UQAM: Université du Québec à Montréal)H-Index: 8
Last. Walter Schaeken (Katholieke Universiteit Leuven)H-Index: 27
view all 3 authors...
Abstract null null The false dilemma or dichotomy is a logical fallacy that occurs when interlocuters accept the premises in an incompatibility statement as being jointly exhaustive (i.e., leaving no third option), whereas that is in fact not the case. Brisson et al. [Memory & Cognition (2018), Vol. 46, pp. 657–670] investigated this fallacy in an adult sample and discovered a content effect that influenced participants’ performance. The current study aimed to elaborate on these findings by esta...
#1J. Pérez–González (University of Extremadura)
#2G. Rey Gozalo (University of Extremadura)H-Index: 4
Last. J.M. Barrigón Morillas (University of Extremadura)H-Index: 9
view all 5 authors...
#1Shuo Chen (AAFC: Agriculture and Agri-Food Canada)H-Index: 6
#2Barbara J. Cade-Menun (AAFC: Agriculture and Agri-Food Canada)H-Index: 41
Last. Qing Chen (CAU: China Agricultural University)H-Index: 3
view all 6 authors...
Abstract null null Efficient phosphorus (P) management is important for crop production and environmental sustainability of cropping systems. The effects of the agricultural management practices of crop rotation and fertilization on soil P forms and cycling were investigated in plots in Swift Current, SK, Canada, under three crop rotation phases [fallow (F), wheat after fallow (WF), and wheat after wheat (WW)], with four nitrogen (N) and/or P fertilizer treatments (+N+P, -N+P, +N-P, and -N-P), u...
#1Kaire Põder (Estonian Business School)H-Index: 6
#2Triin Lauri (University of Konstanz)H-Index: 1
This study investigates civic and citizenship education in a unique post-Communist context–in the bilingual education system of Estonia. Estonia continues to have a bilingual school system where there are Estonian and Russian language schools in parallel. While Estonian language school students are ranked very high in international comparisons, there is a significant difference between the achievement of Estonian and Russian language school students. We claim that this minority achievement gap i...
#1Marja I. Roslund (UH: University of Helsinki)H-Index: 10
#2Riikka Puhakka (UH: University of Helsinki)H-Index: 14
Last. Olli H. Laitinen (UTA: University of Tampere)H-Index: 32
view all 0 authors...
Abstract null null Background null In modern urban environments children have a high incidence of inflammatory disorders, including allergies, asthma, and type 1 diabetes. The underlying cause of these disorders, according to the biodiversity hypothesis, is an imbalance in immune regulation caused by a weak interaction with environmental microbes. In this 2-year study, we analyzed bacterial community shifts in the soil surface in day-care centers and commensal bacteria inhabiting the mouth, skin...
#1Paolo Zuccarini (CSIC: Spanish National Research Council)H-Index: 9
#2Dolores Asensio (CSIC: Spanish National Research Council)H-Index: 21
Last. Josep Peñuelas (CSIC: Spanish National Research Council)H-Index: 135
view all 0 authors...
Abstract null null N deposition is a key driver of global environmental change (GEC). The effects of N addition on soil enzymes which play a key role in biogeochemical cycles and functioning of terrestrial ecosystems remain unclear, especially in vulnerable ecosystems such as P- or/and water-limited Mediterranean ecosystems. Here, we studied seasonal effects of N deposition on the activities of five soil extracellular enzymes in a semi-arid, P-deficient Mediterranean shrubland, and the associati...
#1Christina Knudson (University of St. Thomas (Minnesota))H-Index: 5
#2Sydney Benson (UMN: University of Minnesota)H-Index: 1
Last. Galin L. Jones (UMN: University of Minnesota)H-Index: 26
view all 4 authors...
3 CitationsSource
#1K.H.T. SchrevenH-Index: 1
Last. Bart A. NoletH-Index: 44
view all 4 authors...