<HTML dir=ltr><HEAD>
<META content="text/html; charset=unicode" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 11.00.9600.17280"></HEAD>
<BODY>
<DIV id=idOWAReplyText89502 dir=ltr>
<DIV dir=ltr><FONT color=#000000 size=2 face=Arial>Dear Subscribers,</FONT></DIV>
<DIV dir=ltr><FONT size=2 face=Arial></FONT> </DIV>
<DIV dir=ltr><FONT size=2 face=Arial>I've encountered the following error message. Edad_c is numeric. Estudic, sexo, deprinf and benvii are dichotomic. Regalf has three categories.</FONT></DIV>
<DIV dir=ltr><FONT size=2 face=Arial></FONT> </DIV>
<DIV dir=ltr><FONT size=2 face=Arial>Any help would be very valuable.</FONT></DIV>
<DIV dir=ltr><FONT size=2 face=Arial></FONT> </DIV>
<DIV dir=ltr><FONT size=2 face=Arial>Best regards,</FONT></DIV>
<DIV dir=ltr><FONT size=2 face=Arial></FONT> </DIV>
<DIV dir=ltr><FONT size=2 face=Arial>Angel Rodriguez-Laso</FONT></DIV>
<DIV dir=ltr><FONT size=2 face=Arial></FONT> </DIV>
<DIV dir=ltr><SPAN class=Apple-style-span style="WHITE-SPACE: pre-wrap; WORD-SPACING: 0px; BORDER-COLLAPSE: separate; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); FONT: 13px/15px 'Lucida Console'; ORPHANS: 2; WIDOWS: 2; LETTER-SPACING: normal; BACKGROUND-COLOR: rgb(225,226,229); TEXT-INDENT: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px"><PRE tabIndex=0 id=rstudio_console_output class=GEWYW5YBFEB style="BORDER-TOP-STYLE: none; FONT-SIZE: 10pt !important; FONT-FAMILY: 'Lucida Console'; WHITE-SPACE: pre-wrap !important; BORDER-BOTTOM-STYLE: none; WORD-BREAK: break-all; BORDER-RIGHT-STYLE: none; OUTLINE-STYLE: none; MARGIN: 0px; BORDER-LEFT-STYLE: none; OUTLINE-COLOR: ; LINE-HEIGHT: 1.2; -webkit-user-select: text"><SPAN class="GEWYW5YBJEB ace_keyword" style="WHITE-SPACE: pre; COLOR: blue; -webkit-user-select: text">> </SPAN><SPAN class="GEWYW5YBMDB ace_keyword" style="COLOR: blue">med.fit <- glm(estudic ~ edad_c + sexo + regalf + deprinf, family="binomial" ,data=child65)
</SPAN><SPAN class="GEWYW5YBJEB ace_keyword" style="WHITE-SPACE: pre; COLOR: blue; -webkit-user-select: text">> </SPAN><SPAN class="GEWYW5YBMDB ace_keyword" style="COLOR: blue">out.fit <- glm(benvii ~ edad_c + sexo + regalf + deprinf + estudic, family="binomial" ,data=child65)
</SPAN><SPAN class="GEWYW5YBJEB ace_keyword" style="WHITE-SPACE: pre; COLOR: blue; -webkit-user-select: text">> </SPAN><SPAN class="GEWYW5YBMDB ace_keyword" style="COLOR: blue">library("mediation")
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">Loading required package: MASS
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">Loading required package: Matrix
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">Loading required package: lpSolve
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">Loading required package: mvtnorm
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">mediation: Causal Mediation Analysis
Version: 4.4.2
</SPAN><SPAN class="GEWYW5YBJEB ace_keyword" style="WHITE-SPACE: pre; COLOR: blue; -webkit-user-select: text">> </SPAN><SPAN class="GEWYW5YBMDB ace_keyword" style="COLOR: blue">library("sandwich")
</SPAN><SPAN class="GEWYW5YBJEB ace_keyword" style="WHITE-SPACE: pre; COLOR: blue; -webkit-user-select: text">> </SPAN><SPAN class="GEWYW5YBMDB ace_keyword" style="COLOR: blue">set.seed(2014)
</SPAN><SPAN class="GEWYW5YBJEB ace_keyword" style="WHITE-SPACE: pre; COLOR: blue; -webkit-user-select: text">> </SPAN><SPAN class="GEWYW5YBMDB ace_keyword" style="COLOR: blue">med.out <- mediate(med.fit, out.fit, treat = "deprinf", mediator = "estudic", robustSE = TRUE, sims = 100)
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">Error in if (xhat == 0) out <- 1 else { :
missing value where TRUE/FALSE needed
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">In addition: </SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">Warning message:
</SPAN><SPAN class="GEWYW5YBAEB ace_constant" style="COLOR: rgb(197,6,11)">In mediate(med.fit, out.fit, treat = "deprinf", mediator = "estudic", :
treatment and control values do not match factor levels; using no and sí as control and treatment, respectively
</SPAN></SPAN><BR><HR tabIndex=-1></PRE></DIV></DIV></BODY></HTML>