<div dir="ltr"><div>Hi,<br><br></div>I've previously used this package to fit simple linear models and it worked very well. But recently I've been trying to do this with multilevel models (lmer) and I'm getting this message: <br><br>Error in mediation::mediate(model.m = mediatorModel, model.y = outcomeModel,  : <br>  mediator model is not yet implemented<br><div><br></div><div>Below are my models. Can anyone tell why I'm getting that error message?<br></div><div><br>mediatorModel <- lmer(m ~ x1 + x2 + (1 | id), data = df)<br><br>outcomeModel <- lmer(y ~ m + x1 + x2 + (1 | id), data = df)<br><br>mediation <- mediate(model.m = mediatorModel, model.y = outcomeModel, sims = 100, robustSE = F, treat = "x2", mediator = "m")<br><br></div><div>I really appreciate any help anyone can provide. Thanks!<br><br></div><div>Best,<br>Hause<br></div><div><br></div></div>