Quantile regression in linear mixed models: A stochastic approximation EM approach

Christian E. Galarza, Victor H. Lachos, Dipankar Bandyopadhyay

Research output: Contribution to journalArticlepeer-review

12 Scopus citations


This paper develops a likelihood-based approach to analyze quantile regression (QR) models for continuous longitudinal data via the asymmetric Laplace distribution (ALD). Compared to the conventional mean regression approach, QR can characterize the entire conditional distribution of the outcome variable and is more robust to the presence of outliers and misspecification of the error distribution. Exploiting the nice hierarchical representation of the ALD, our classical approach follows a Stochastic Approximation of the EM (SAEM) algorithm in deriving exact maximum likelihood estimates of the fixed-effects and variance components. We evaluate the finite sample performance of the algorithm and the asymptotic properties of the ML estimates through empirical experiments and applications to two real life datasets. Our empirical results clearly indicate that the SAEM estimates outperforms the estimates obtained via the combination of Gaussian quadrature and non-smooth optimization routines of the Geraci and Bottai (2014) approach in terms of standard errors and mean square error. The proposed SAEM algorithm is implemented in the R package qrLMM().

Original languageEnglish (US)
Pages (from-to)471-482
Number of pages12
JournalStatistics and its Interface
Issue number3
StatePublished - 2017


  • Asymmetric laplace distribution
  • Linear mixed-effects models
  • Quantile regression
  • SAEM algorithm


Dive into the research topics of 'Quantile regression in linear mixed models: A stochastic approximation EM approach'. Together they form a unique fingerprint.

Cite this