bmlm: An R package for Bayesian MultiLevel Mediation models

bmlm is an R package providing convenient methods for Bayesian estimation of multilevel mediation models using Stan.

For more details, visit the package's website.


To install the latest stable version of bmlm from CRAN, run


Development version

The latest development version of bmlm requires devtools for installation. If you don't have the devtools package installed in R, first run this line:


Then proceed to install bmlm:

devtools::install_github("mvuorre/bmlm", args = "--preclean")


bmlm ships with an example data set from Intensive Longitudinal Methods: An Introduction to Diary and Experience Sampling Research (Bolger & Laurenceau, 2013, chapter 9). To estimate the multilevel mediation model presented in that chapter, run:

fit <- mlm(BLch9)

After a while, you will have a joint posterior distribution of plausible parameter values from the model applied to this data. Inspect the model:



Please contact the author of the package for questions and suggestions. I recommend creating a new issue on GitHub.


If you use this software, please cite it:

#> To cite package bmlm in publications, please use:
#>   Vuorre, M., (2016). bmlm: Bayesian Multilevel Mediation. R
#>   package version 1.2.1.
#> A BibTeX entry for LaTeX users is
#>   @Manual{,
#>     title = {bmlm: Bayesian Multilevel Mediation},
#>     author = {Matti Vuorre},
#>     year = {2016},
#>     url = {},
#>   }