<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Dear R mediation Folk—I would like to start by saying that I have been enjoying the functionality and flexibility of the mediation package that Kosuke, Dustin, Teppei, Luke, et al. have developed. Bravo and keep up the important work of disseminating this tool! <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Nevertheless, I ran into a surprising error just this morning when running mediate on some lme4 data runs. It looks like an issue with the coding rather than my data because I get the same error when running the multilevel modeling example on the “jobs” data from the reference manual. This is what I get:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>> model.m <- glmer(job_dich ~ treat + econ_hard + (1 | educ), family = binomial(link = "probit"), data = jobs)<o:p></o:p></p><p class=MsoNormal>> model.y <- glmer(work1 ~ treat + job_dich + econ_hard + (1 + treat | occp), family = binomial(link = "probit"), data = jobs)<o:p></o:p></p><p class=MsoNormal>> multilevel <- mediate(model.m, model.y, treat = "treat", mediator = "job_dich", sims=50, group.out="educ")<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Error in mvrnorm(sims, mu = MModel.coef, Sigma = MModel.var.cov) : incompatible arguments<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am still somewhat new to R, so any insights on what needs to be tweaked would be much appreciated.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Best, J-P<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Jean-Philippe Laurenceau, Ph.D.<o:p></o:p></p><p class=MsoNormal>Department of Psychology<o:p></o:p></p><p class=MsoNormal>University of Delaware<o:p></o:p></p><p class=MsoNormal>110 Wolf Hall<o:p></o:p></p><p class=MsoNormal>Newark, DE 19716-2577<o:p></o:p></p><p class=MsoNormal>Voice: (302) 831-2309<o:p></o:p></p><p class=MsoNormal>Fax: (302) 831-3645<o:p></o:p></p><p class=MsoNormal><a href="mailto:jlaurenceau@psych.udel.edu">jlaurenceau@psych.udel.edu</a> <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>