https://github.com/cran/bayestestR
Raw File
Tip revision: d73e4a2afcfbd6402c11716877e8f7466f309ef4 authored by Dominique Makowski on 22 October 2020, 13:40:02 UTC
version 0.7.5
Tip revision: d73e4a2
mediation.html
<!DOCTYPE html>

<html xmlns="http://www.w3.org/1999/xhtml">

<head>

<meta charset="utf-8">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<meta name="generator" content="pandoc" />
<meta name="viewport" content="width=device-width, initial-scale=1">

<style type="text/css">
@font-face {
font-family: octicons-link;
src: url(data:font/woff;charset=utf-8;base64,d09GRgABAAAAAAZwABAAAAAACFQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABEU0lHAAAGaAAAAAgAAAAIAAAAAUdTVUIAAAZcAAAACgAAAAoAAQAAT1MvMgAAAyQAAABJAAAAYFYEU3RjbWFwAAADcAAAAEUAAACAAJThvmN2dCAAAATkAAAABAAAAAQAAAAAZnBnbQAAA7gAAACyAAABCUM+8IhnYXNwAAAGTAAAABAAAAAQABoAI2dseWYAAAFsAAABPAAAAZwcEq9taGVhZAAAAsgAAAA0AAAANgh4a91oaGVhAAADCAAAABoAAAAkCA8DRGhtdHgAAAL8AAAADAAAAAwGAACfbG9jYQAAAsAAAAAIAAAACABiATBtYXhwAAACqAAAABgAAAAgAA8ASm5hbWUAAAToAAABQgAAAlXu73sOcG9zdAAABiwAAAAeAAAAME3QpOBwcmVwAAAEbAAAAHYAAAB/aFGpk3jaTY6xa8JAGMW/O62BDi0tJLYQincXEypYIiGJjSgHniQ6umTsUEyLm5BV6NDBP8Tpts6F0v+k/0an2i+itHDw3v2+9+DBKTzsJNnWJNTgHEy4BgG3EMI9DCEDOGEXzDADU5hBKMIgNPZqoD3SilVaXZCER3/I7AtxEJLtzzuZfI+VVkprxTlXShWKb3TBecG11rwoNlmmn1P2WYcJczl32etSpKnziC7lQyWe1smVPy/Lt7Kc+0vWY/gAgIIEqAN9we0pwKXreiMasxvabDQMM4riO+qxM2ogwDGOZTXxwxDiycQIcoYFBLj5K3EIaSctAq2kTYiw+ymhce7vwM9jSqO8JyVd5RH9gyTt2+J/yUmYlIR0s04n6+7Vm1ozezUeLEaUjhaDSuXHwVRgvLJn1tQ7xiuVv/ocTRF42mNgZGBgYGbwZOBiAAFGJBIMAAizAFoAAABiAGIAznjaY2BkYGAA4in8zwXi+W2+MjCzMIDApSwvXzC97Z4Ig8N/BxYGZgcgl52BCSQKAA3jCV8CAABfAAAAAAQAAEB42mNgZGBg4f3vACQZQABIMjKgAmYAKEgBXgAAeNpjYGY6wTiBgZWBg2kmUxoDA4MPhGZMYzBi1AHygVLYQUCaawqDA4PChxhmh/8ODDEsvAwHgMKMIDnGL0x7gJQCAwMAJd4MFwAAAHjaY2BgYGaA4DAGRgYQkAHyGMF8NgYrIM3JIAGVYYDT+AEjAwuDFpBmA9KMDEwMCh9i/v8H8sH0/4dQc1iAmAkALaUKLgAAAHjaTY9LDsIgEIbtgqHUPpDi3gPoBVyRTmTddOmqTXThEXqrob2gQ1FjwpDvfwCBdmdXC5AVKFu3e5MfNFJ29KTQT48Ob9/lqYwOGZxeUelN2U2R6+cArgtCJpauW7UQBqnFkUsjAY/kOU1cP+DAgvxwn1chZDwUbd6CFimGXwzwF6tPbFIcjEl+vvmM/byA48e6tWrKArm4ZJlCbdsrxksL1AwWn/yBSJKpYbq8AXaaTb8AAHja28jAwOC00ZrBeQNDQOWO//sdBBgYGRiYWYAEELEwMTE4uzo5Zzo5b2BxdnFOcALxNjA6b2ByTswC8jYwg0VlNuoCTWAMqNzMzsoK1rEhNqByEyerg5PMJlYuVueETKcd/89uBpnpvIEVomeHLoMsAAe1Id4AAAAAAAB42oWQT07CQBTGv0JBhagk7HQzKxca2sJCE1hDt4QF+9JOS0nbaaYDCQfwCJ7Au3AHj+LO13FMmm6cl7785vven0kBjHCBhfpYuNa5Ph1c0e2Xu3jEvWG7UdPDLZ4N92nOm+EBXuAbHmIMSRMs+4aUEd4Nd3CHD8NdvOLTsA2GL8M9PODbcL+hD7C1xoaHeLJSEao0FEW14ckxC+TU8TxvsY6X0eLPmRhry2WVioLpkrbp84LLQPGI7c6sOiUzpWIWS5GzlSgUzzLBSikOPFTOXqly7rqx0Z1Q5BAIoZBSFihQYQOOBEdkCOgXTOHA07HAGjGWiIjaPZNW13/+lm6S9FT7rLHFJ6fQbkATOG1j2OFMucKJJsxIVfQORl+9Jyda6Sl1dUYhSCm1dyClfoeDve4qMYdLEbfqHf3O/AdDumsjAAB42mNgYoAAZQYjBmyAGYQZmdhL8zLdDEydARfoAqIAAAABAAMABwAKABMAB///AA8AAQAAAAAAAAAAAAAAAAABAAAAAA==) format('woff');
}
body {
-webkit-text-size-adjust: 100%;
text-size-adjust: 100%;
color: #333;
font-family: "Helvetica Neue", Helvetica, "Segoe UI", Arial, freesans, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", "Segoe UI Symbol";
font-size: 16px;
line-height: 1.6;
word-wrap: break-word;
}
a {
background-color: transparent;
}
a:active,
a:hover {
outline: 0;
}
strong {
font-weight: bold;
}
h1 {
font-size: 2em;
margin: 0.67em 0;
}
img {
border: 0;
}
hr {
box-sizing: content-box;
height: 0;
}
pre {
overflow: auto;
}
code,
kbd,
pre {
font-family: monospace, monospace;
font-size: 1em;
}
input {
color: inherit;
font: inherit;
margin: 0;
}
html input[disabled] {
cursor: default;
}
input {
line-height: normal;
}
input[type="checkbox"] {
box-sizing: border-box;
padding: 0;
}
table {
border-collapse: collapse;
border-spacing: 0;
}
td,
th {
padding: 0;
}
* {
box-sizing: border-box;
}
input {
font: 13px / 1.4 Helvetica, arial, nimbussansl, liberationsans, freesans, clean, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", "Segoe UI Symbol";
}
a {
color: #4078c0;
text-decoration: none;
}
a:hover,
a:active {
text-decoration: underline;
}
hr {
height: 0;
margin: 15px 0;
overflow: hidden;
background: transparent;
border: 0;
border-bottom: 1px solid #ddd;
}
hr:before {
display: table;
content: "";
}
hr:after {
display: table;
clear: both;
content: "";
}
h1,
h2,
h3,
h4,
h5,
h6 {
margin-top: 15px;
margin-bottom: 15px;
line-height: 1.1;
}
h1 {
font-size: 30px;
}
h2 {
font-size: 21px;
}
h3 {
font-size: 16px;
}
h4 {
font-size: 14px;
}
h5 {
font-size: 12px;
}
h6 {
font-size: 11px;
}
blockquote {
margin: 0;
}
ul,
ol {
padding: 0;
margin-top: 0;
margin-bottom: 0;
}
ol ol,
ul ol {
list-style-type: lower-roman;
}
ul ul ol,
ul ol ol,
ol ul ol,
ol ol ol {
list-style-type: lower-alpha;
}
dd {
margin-left: 0;
}
code {
font-family: Consolas, "Liberation Mono", Menlo, Courier, monospace;
font-size: 12px;
}
pre {
margin-top: 0;
margin-bottom: 0;
font: 12px Consolas, "Liberation Mono", Menlo, Courier, monospace;
}
.select::-ms-expand {
opacity: 0;
}
.octicon {
font: normal normal normal 16px/1 octicons-link;
display: inline-block;
text-decoration: none;
text-rendering: auto;
-webkit-font-smoothing: antialiased;
-moz-osx-font-smoothing: grayscale;
-webkit-user-select: none;
-moz-user-select: none;
-ms-user-select: none;
user-select: none;
}
.octicon-link:before {
content: '\f05c';
}
.markdown-body:before {
display: table;
content: "";
}
.markdown-body:after {
display: table;
clear: both;
content: "";
}
.markdown-body>*:first-child {
margin-top: 0 !important;
}
.markdown-body>*:last-child {
margin-bottom: 0 !important;
}
a:not([href]) {
color: inherit;
text-decoration: none;
}
.anchor {
display: inline-block;
padding-right: 2px;
margin-left: -18px;
}
.anchor:focus {
outline: none;
}
h1,
h2,
h3,
h4,
h5,
h6 {
margin-top: 1em;
margin-bottom: 16px;
font-weight: bold;
line-height: 1.4;
}
h1 .octicon-link,
h2 .octicon-link,
h3 .octicon-link,
h4 .octicon-link,
h5 .octicon-link,
h6 .octicon-link {
color: #000;
vertical-align: middle;
visibility: hidden;
}
h1:hover .anchor,
h2:hover .anchor,
h3:hover .anchor,
h4:hover .anchor,
h5:hover .anchor,
h6:hover .anchor {
text-decoration: none;
}
h1:hover .anchor .octicon-link,
h2:hover .anchor .octicon-link,
h3:hover .anchor .octicon-link,
h4:hover .anchor .octicon-link,
h5:hover .anchor .octicon-link,
h6:hover .anchor .octicon-link {
visibility: visible;
}
h1 {
padding-bottom: 0.3em;
font-size: 2.25em;
line-height: 1.2;
border-bottom: 1px solid #eee;
}
h1 .anchor {
line-height: 1;
}
h2 {
padding-bottom: 0.3em;
font-size: 1.75em;
line-height: 1.225;
border-bottom: 1px solid #eee;
}
h2 .anchor {
line-height: 1;
}
h3 {
font-size: 1.5em;
line-height: 1.43;
}
h3 .anchor {
line-height: 1.2;
}
h4 {
font-size: 1.25em;
}
h4 .anchor {
line-height: 1.2;
}
h5 {
font-size: 1em;
}
h5 .anchor {
line-height: 1.1;
}
h6 {
font-size: 1em;
color: #777;
}
h6 .anchor {
line-height: 1.1;
}
p,
blockquote,
ul,
ol,
dl,
table,
pre {
margin-top: 0;
margin-bottom: 16px;
}
hr {
height: 4px;
padding: 0;
margin: 16px 0;
background-color: #e7e7e7;
border: 0 none;
}
ul,
ol {
padding-left: 2em;
}
ul ul,
ul ol,
ol ol,
ol ul {
margin-top: 0;
margin-bottom: 0;
}
li>p {
margin-top: 16px;
}
dl {
padding: 0;
}
dl dt {
padding: 0;
margin-top: 16px;
font-size: 1em;
font-style: italic;
font-weight: bold;
}
dl dd {
padding: 0 16px;
margin-bottom: 16px;
}
blockquote {
padding: 0 15px;
color: #777;
border-left: 4px solid #ddd;
}
blockquote>:first-child {
margin-top: 0;
}
blockquote>:last-child {
margin-bottom: 0;
}
table {
display: block;
width: 100%;
overflow: auto;
word-break: normal;
word-break: keep-all;
}
table th {
font-weight: bold;
}
table th,
table td {
padding: 6px 13px;
border: 1px solid #ddd;
}
table tr {
background-color: #fff;
border-top: 1px solid #ccc;
}
table tr:nth-child(2n) {
background-color: #f8f8f8;
}
img {
max-width: 100%;
box-sizing: content-box;
background-color: #fff;
}
code {
padding: 0;
padding-top: 0.2em;
padding-bottom: 0.2em;
margin: 0;
font-size: 85%;
background-color: rgba(0,0,0,0.04);
border-radius: 3px;
}
code:before,
code:after {
letter-spacing: -0.2em;
content: "\00a0";
}
pre>code {
padding: 0;
margin: 0;
font-size: 100%;
word-break: normal;
white-space: pre;
background: transparent;
border: 0;
}
.highlight {
margin-bottom: 16px;
}
.highlight pre,
pre {
padding: 16px;
overflow: auto;
font-size: 85%;
line-height: 1.45;
background-color: #f7f7f7;
border-radius: 3px;
}
.highlight pre {
margin-bottom: 0;
word-break: normal;
}
pre {
word-wrap: normal;
}
pre code {
display: inline;
max-width: initial;
padding: 0;
margin: 0;
overflow: initial;
line-height: inherit;
word-wrap: normal;
background-color: transparent;
border: 0;
}
pre code:before,
pre code:after {
content: normal;
}
kbd {
display: inline-block;
padding: 3px 5px;
font-size: 11px;
line-height: 10px;
color: #555;
vertical-align: middle;
background-color: #fcfcfc;
border: solid 1px #ccc;
border-bottom-color: #bbb;
border-radius: 3px;
box-shadow: inset 0 -1px 0 #bbb;
}
.pl-c {
color: #969896;
}
.pl-c1,
.pl-s .pl-v {
color: #0086b3;
}
.pl-e,
.pl-en {
color: #795da3;
}
.pl-s .pl-s1,
.pl-smi {
color: #333;
}
.pl-ent {
color: #63a35c;
}
.pl-k {
color: #a71d5d;
}
.pl-pds,
.pl-s,
.pl-s .pl-pse .pl-s1,
.pl-sr,
.pl-sr .pl-cce,
.pl-sr .pl-sra,
.pl-sr .pl-sre {
color: #183691;
}
.pl-v {
color: #ed6a43;
}
.pl-id {
color: #b52a1d;
}
.pl-ii {
background-color: #b52a1d;
color: #f8f8f8;
}
.pl-sr .pl-cce {
color: #63a35c;
font-weight: bold;
}
.pl-ml {
color: #693a17;
}
.pl-mh,
.pl-mh .pl-en,
.pl-ms {
color: #1d3e81;
font-weight: bold;
}
.pl-mq {
color: #008080;
}
.pl-mi {
color: #333;
font-style: italic;
}
.pl-mb {
color: #333;
font-weight: bold;
}
.pl-md {
background-color: #ffecec;
color: #bd2c00;
}
.pl-mi1 {
background-color: #eaffea;
color: #55a532;
}
.pl-mdr {
color: #795da3;
font-weight: bold;
}
.pl-mo {
color: #1d3e81;
}
kbd {
display: inline-block;
padding: 3px 5px;
font: 11px Consolas, "Liberation Mono", Menlo, Courier, monospace;
line-height: 10px;
color: #555;
vertical-align: middle;
background-color: #fcfcfc;
border: solid 1px #ccc;
border-bottom-color: #bbb;
border-radius: 3px;
box-shadow: inset 0 -1px 0 #bbb;
}
.task-list-item {
list-style-type: none;
}
.task-list-item+.task-list-item {
margin-top: 3px;
}
.task-list-item input {
margin: 0 0.35em 0.25em -1.6em;
vertical-align: middle;
}
:checked+.radio-label {
z-index: 1;
position: relative;
border-color: #4078c0;
}
.sourceLine {
display: inline-block;
}
code .kw { color: #000000; }
code .dt { color: #ed6a43; }
code .dv { color: #009999; }
code .bn { color: #009999; }
code .fl { color: #009999; }
code .ch { color: #009999; }
code .st { color: #183691; }
code .co { color: #969896; }
code .ot { color: #0086b3; }
code .al { color: #a61717; }
code .fu { color: #63a35c; }
code .er { color: #a61717; background-color: #e3d2d2; }
code .wa { color: #000000; }
code .cn { color: #008080; }
code .sc { color: #008080; }
code .vs { color: #183691; }
code .ss { color: #183691; }
code .im { color: #000000; }
code .va {color: #008080; }
code .cf { color: #000000; }
code .op { color: #000000; }
code .bu { color: #000000; }
code .ex { color: #000000; }
code .pp { color: #999999; }
code .at { color: #008080; }
code .do { color: #969896; }
code .an { color: #008080; }
code .cv { color: #008080; }
code .in { color: #008080; }
</style>
<style>
body {
  box-sizing: border-box;
  min-width: 200px;
  max-width: 980px;
  margin: 0 auto;
  padding: 45px;
  padding-top: 0px;
}
</style>

</head>

<body>

<h1 id="summary-of-mediation-analysis-using-bayesian-regression-models">Summary of Mediation Analysis using Bayesian Regression Models</h1>
<p>This vignettes demonstrates the <code>mediation()</code>-function. Before we start, we fit some models, including a mediation-object from the <em>mediation</em>-package, which we use for comparison with <em>brms</em> and <em>rstanarm</em>.</p>
<div class="sourceCode" id="cb1"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb1-1"><a href="#cb1-1"></a><span class="kw">library</span>(bayestestR)</span>
<span id="cb1-2"><a href="#cb1-2"></a><span class="kw">library</span>(mediation)</span>
<span id="cb1-3"><a href="#cb1-3"></a><span class="kw">library</span>(brms)</span>
<span id="cb1-4"><a href="#cb1-4"></a><span class="kw">library</span>(rstanarm)</span>
<span id="cb1-5"><a href="#cb1-5"></a></span>
<span id="cb1-6"><a href="#cb1-6"></a><span class="co"># load sample data</span></span>
<span id="cb1-7"><a href="#cb1-7"></a><span class="kw">data</span>(jobs)</span>
<span id="cb1-8"><a href="#cb1-8"></a></span>
<span id="cb1-9"><a href="#cb1-9"></a><span class="kw">set.seed</span>(<span class="dv">123</span>)</span>
<span id="cb1-10"><a href="#cb1-10"></a><span class="co"># linear models, for mediation analysis</span></span>
<span id="cb1-11"><a href="#cb1-11"></a>b1 &lt;-<span class="st"> </span><span class="kw">lm</span>(job_seek <span class="op">~</span><span class="st"> </span>treat <span class="op">+</span><span class="st"> </span>econ_hard <span class="op">+</span><span class="st"> </span>sex <span class="op">+</span><span class="st"> </span>age, <span class="dt">data =</span> jobs)</span>
<span id="cb1-12"><a href="#cb1-12"></a>b2 &lt;-<span class="st"> </span><span class="kw">lm</span>(depress2 <span class="op">~</span><span class="st"> </span>treat <span class="op">+</span><span class="st"> </span>job_seek <span class="op">+</span><span class="st"> </span>econ_hard <span class="op">+</span><span class="st"> </span>sex <span class="op">+</span><span class="st"> </span>age, <span class="dt">data =</span> jobs)</span>
<span id="cb1-13"><a href="#cb1-13"></a></span>
<span id="cb1-14"><a href="#cb1-14"></a><span class="co"># mediation analysis, for comparison with brms</span></span>
<span id="cb1-15"><a href="#cb1-15"></a>m1 &lt;-<span class="st"> </span><span class="kw">mediate</span>(b1, b2, <span class="dt">sims =</span> <span class="dv">1000</span>, <span class="dt">treat =</span> <span class="st">&quot;treat&quot;</span>, <span class="dt">mediator =</span> <span class="st">&quot;job_seek&quot;</span>)</span></code></pre></div>
<div class="sourceCode" id="cb2"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb2-1"><a href="#cb2-1"></a><span class="co"># Fit Bayesian mediation model in brms</span></span>
<span id="cb2-2"><a href="#cb2-2"></a>f1 &lt;-<span class="st"> </span><span class="kw">bf</span>(job_seek <span class="op">~</span><span class="st"> </span>treat <span class="op">+</span><span class="st"> </span>econ_hard <span class="op">+</span><span class="st"> </span>sex <span class="op">+</span><span class="st"> </span>age)</span>
<span id="cb2-3"><a href="#cb2-3"></a>f2 &lt;-<span class="st"> </span><span class="kw">bf</span>(depress2 <span class="op">~</span><span class="st"> </span>treat <span class="op">+</span><span class="st"> </span>job_seek <span class="op">+</span><span class="st"> </span>econ_hard <span class="op">+</span><span class="st"> </span>sex <span class="op">+</span><span class="st"> </span>age)</span>
<span id="cb2-4"><a href="#cb2-4"></a>m2 &lt;-<span class="st"> </span><span class="kw">brm</span>(f1 <span class="op">+</span><span class="st"> </span>f2 <span class="op">+</span><span class="st"> </span><span class="kw">set_rescor</span>(<span class="ot">FALSE</span>), <span class="dt">data =</span> jobs, <span class="dt">cores =</span> <span class="dv">4</span>)</span></code></pre></div>
<div class="sourceCode" id="cb3"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb3-1"><a href="#cb3-1"></a><span class="co"># Fit Bayesian mediation model in rstanarm</span></span>
<span id="cb3-2"><a href="#cb3-2"></a>m3 &lt;-<span class="st"> </span><span class="kw">stan_mvmer</span>(</span>
<span id="cb3-3"><a href="#cb3-3"></a>  <span class="kw">list</span>(job_seek <span class="op">~</span><span class="st"> </span>treat <span class="op">+</span><span class="st"> </span>econ_hard <span class="op">+</span><span class="st"> </span>sex <span class="op">+</span><span class="st"> </span>age <span class="op">+</span><span class="st"> </span>(<span class="dv">1</span> <span class="op">|</span><span class="st"> </span>occp),</span>
<span id="cb3-4"><a href="#cb3-4"></a>       depress2 <span class="op">~</span><span class="st"> </span>treat <span class="op">+</span><span class="st"> </span>job_seek <span class="op">+</span><span class="st"> </span>econ_hard <span class="op">+</span><span class="st"> </span>sex <span class="op">+</span><span class="st"> </span>age <span class="op">+</span><span class="st"> </span>(<span class="dv">1</span> <span class="op">|</span><span class="st"> </span>occp)),</span>
<span id="cb3-5"><a href="#cb3-5"></a>  <span class="dt">data =</span> jobs,</span>
<span id="cb3-6"><a href="#cb3-6"></a>  <span class="dt">cores =</span> <span class="dv">4</span>,</span>
<span id="cb3-7"><a href="#cb3-7"></a>  <span class="dt">refresh =</span> <span class="dv">0</span></span>
<span id="cb3-8"><a href="#cb3-8"></a>)</span></code></pre></div>
<p><code>mediation()</code> is a summary function, especially for mediation analysis, i.e. for multivariate response models with casual mediation effects.</p>
<p>In the models <code>m2</code> and <code>m3</code>, <code>treat</code> is the treatment effect and <code>job_seek</code> is the mediator effect. For the <em>brms</em> model (<code>m2</code>), <code>f1</code> describes the mediator model and <code>f2</code> describes the outcome model. This is similar for the <em>rstanarm</em> model.</p>
<p><code>mediation()</code> returns a data frame with information on the <em>direct effect</em> (median value of posterior samples from treatment of the outcome model), <em>mediator effect</em> (median value of posterior samples from mediator of the outcome model), <em>indirect effect</em> (median value of the multiplication of the posterior samples from mediator of the outcome model and the posterior samples from treatment of the mediation model) and the <em>total effect</em> (median value of sums of posterior samples used for the direct and indirect effect). The <em>proportion mediated</em> is the indirect effect divided by the total effect.</p>
<p>The simplest call just needs the model-object.</p>
<div class="sourceCode" id="cb4"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb4-1"><a href="#cb4-1"></a><span class="co"># for brms</span></span>
<span id="cb4-2"><a href="#cb4-2"></a><span class="kw">mediation</span>(m2)</span>
<span id="cb4-3"><a href="#cb4-3"></a><span class="co">#&gt; # Causal Mediation Analysis for Stan Model</span></span>
<span id="cb4-4"><a href="#cb4-4"></a><span class="co">#&gt; </span></span>
<span id="cb4-5"><a href="#cb4-5"></a><span class="co">#&gt;   Treatment: treat</span></span>
<span id="cb4-6"><a href="#cb4-6"></a><span class="co">#&gt;   Mediator : job_seek</span></span>
<span id="cb4-7"><a href="#cb4-7"></a><span class="co">#&gt;   Response : depress2</span></span>
<span id="cb4-8"><a href="#cb4-8"></a><span class="co">#&gt; </span></span>
<span id="cb4-9"><a href="#cb4-9"></a><span class="co">#&gt; Effect                 | Estimate |          89% ETI</span></span>
<span id="cb4-10"><a href="#cb4-10"></a><span class="co">#&gt; ----------------------------------------------------</span></span>
<span id="cb4-11"><a href="#cb4-11"></a><span class="co">#&gt; Direct Effect (ADE)    |   -0.040 | [-0.110,  0.031]</span></span>
<span id="cb4-12"><a href="#cb4-12"></a><span class="co">#&gt; Indirect Effect (ACME) |   -0.015 | [-0.036,  0.004]</span></span>
<span id="cb4-13"><a href="#cb4-13"></a><span class="co">#&gt; Mediator Effect        |   -0.240 | [-0.285, -0.195]</span></span>
<span id="cb4-14"><a href="#cb4-14"></a><span class="co">#&gt; Total Effect           |   -0.055 | [-0.129,  0.018]</span></span>
<span id="cb4-15"><a href="#cb4-15"></a><span class="co">#&gt; </span></span>
<span id="cb4-16"><a href="#cb4-16"></a><span class="co">#&gt; Proportion mediated: 28.14% [-71.11%, 127.40%]</span></span>
<span id="cb4-17"><a href="#cb4-17"></a></span>
<span id="cb4-18"><a href="#cb4-18"></a><span class="co"># for rstanarm</span></span>
<span id="cb4-19"><a href="#cb4-19"></a><span class="kw">mediation</span>(m3)</span>
<span id="cb4-20"><a href="#cb4-20"></a><span class="co">#&gt; # Causal Mediation Analysis for Stan Model</span></span>
<span id="cb4-21"><a href="#cb4-21"></a><span class="co">#&gt; </span></span>
<span id="cb4-22"><a href="#cb4-22"></a><span class="co">#&gt;   Treatment: treat</span></span>
<span id="cb4-23"><a href="#cb4-23"></a><span class="co">#&gt;   Mediator : job_seek</span></span>
<span id="cb4-24"><a href="#cb4-24"></a><span class="co">#&gt;   Response : depress2</span></span>
<span id="cb4-25"><a href="#cb4-25"></a><span class="co">#&gt; </span></span>
<span id="cb4-26"><a href="#cb4-26"></a><span class="co">#&gt; Effect                 | Estimate |          89% ETI</span></span>
<span id="cb4-27"><a href="#cb4-27"></a><span class="co">#&gt; ----------------------------------------------------</span></span>
<span id="cb4-28"><a href="#cb4-28"></a><span class="co">#&gt; Direct Effect (ADE)    |   -0.040 | [-0.111,  0.031]</span></span>
<span id="cb4-29"><a href="#cb4-29"></a><span class="co">#&gt; Indirect Effect (ACME) |   -0.018 | [-0.037,  0.002]</span></span>
<span id="cb4-30"><a href="#cb4-30"></a><span class="co">#&gt; Mediator Effect        |   -0.241 | [-0.286, -0.197]</span></span>
<span id="cb4-31"><a href="#cb4-31"></a><span class="co">#&gt; Total Effect           |   -0.057 | [-0.130,  0.017]</span></span>
<span id="cb4-32"><a href="#cb4-32"></a><span class="co">#&gt; </span></span>
<span id="cb4-33"><a href="#cb4-33"></a><span class="co">#&gt; Proportion mediated: 30.59% [-75.65%, 136.82%]</span></span></code></pre></div>
<p>Typically, <code>mediation()</code> finds the treatment and mediator variables automatically. If this does not work, use the <code>treatment</code> and <code>mediator</code> arguments to specify the related variable names. For all values, the 89% credible intervals are calculated by default. Use <code>ci</code> to calculate a different interval.</p>
<p>Here is a comparison with the <em>mediation</em> package. Note that the <code>summary()</code>-output of the <em>mediation</em> package shows the indirect effect first, followed by the direct effect.</p>
<div class="sourceCode" id="cb5"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb5-1"><a href="#cb5-1"></a><span class="kw">summary</span>(m1)</span>
<span id="cb5-2"><a href="#cb5-2"></a><span class="co">#&gt; </span></span>
<span id="cb5-3"><a href="#cb5-3"></a><span class="co">#&gt; Causal Mediation Analysis </span></span>
<span id="cb5-4"><a href="#cb5-4"></a><span class="co">#&gt; </span></span>
<span id="cb5-5"><a href="#cb5-5"></a><span class="co">#&gt; Quasi-Bayesian Confidence Intervals</span></span>
<span id="cb5-6"><a href="#cb5-6"></a><span class="co">#&gt; </span></span>
<span id="cb5-7"><a href="#cb5-7"></a><span class="co">#&gt;                Estimate 95% CI Lower 95% CI Upper p-value</span></span>
<span id="cb5-8"><a href="#cb5-8"></a><span class="co">#&gt; ACME            -0.0157      -0.0387         0.01    0.19</span></span>
<span id="cb5-9"><a href="#cb5-9"></a><span class="co">#&gt; ADE             -0.0438      -0.1315         0.04    0.35</span></span>
<span id="cb5-10"><a href="#cb5-10"></a><span class="co">#&gt; Total Effect    -0.0595      -0.1530         0.02    0.21</span></span>
<span id="cb5-11"><a href="#cb5-11"></a><span class="co">#&gt; Prop. Mediated   0.2137      -2.0277         2.70    0.32</span></span>
<span id="cb5-12"><a href="#cb5-12"></a><span class="co">#&gt; </span></span>
<span id="cb5-13"><a href="#cb5-13"></a><span class="co">#&gt; Sample Size Used: 899 </span></span>
<span id="cb5-14"><a href="#cb5-14"></a><span class="co">#&gt; </span></span>
<span id="cb5-15"><a href="#cb5-15"></a><span class="co">#&gt; </span></span>
<span id="cb5-16"><a href="#cb5-16"></a><span class="co">#&gt; Simulations: 1000</span></span>
<span id="cb5-17"><a href="#cb5-17"></a></span>
<span id="cb5-18"><a href="#cb5-18"></a><span class="kw">mediation</span>(m2, <span class="dt">ci =</span> <span class="fl">.95</span>)</span>
<span id="cb5-19"><a href="#cb5-19"></a><span class="co">#&gt; # Causal Mediation Analysis for Stan Model</span></span>
<span id="cb5-20"><a href="#cb5-20"></a><span class="co">#&gt; </span></span>
<span id="cb5-21"><a href="#cb5-21"></a><span class="co">#&gt;   Treatment: treat</span></span>
<span id="cb5-22"><a href="#cb5-22"></a><span class="co">#&gt;   Mediator : job_seek</span></span>
<span id="cb5-23"><a href="#cb5-23"></a><span class="co">#&gt;   Response : depress2</span></span>
<span id="cb5-24"><a href="#cb5-24"></a><span class="co">#&gt; </span></span>
<span id="cb5-25"><a href="#cb5-25"></a><span class="co">#&gt; Effect                 | Estimate |          95% ETI</span></span>
<span id="cb5-26"><a href="#cb5-26"></a><span class="co">#&gt; ----------------------------------------------------</span></span>
<span id="cb5-27"><a href="#cb5-27"></a><span class="co">#&gt; Direct Effect (ADE)    |   -0.040 | [-0.124,  0.046]</span></span>
<span id="cb5-28"><a href="#cb5-28"></a><span class="co">#&gt; Indirect Effect (ACME) |   -0.015 | [-0.041,  0.008]</span></span>
<span id="cb5-29"><a href="#cb5-29"></a><span class="co">#&gt; Mediator Effect        |   -0.240 | [-0.294, -0.185]</span></span>
<span id="cb5-30"><a href="#cb5-30"></a><span class="co">#&gt; Total Effect           |   -0.055 | [-0.145,  0.034]</span></span>
<span id="cb5-31"><a href="#cb5-31"></a><span class="co">#&gt; </span></span>
<span id="cb5-32"><a href="#cb5-32"></a><span class="co">#&gt; Proportion mediated: 28.14% [-181.46%, 237.75%]</span></span>
<span id="cb5-33"><a href="#cb5-33"></a></span>
<span id="cb5-34"><a href="#cb5-34"></a><span class="kw">mediation</span>(m3, <span class="dt">ci =</span> <span class="fl">.95</span>)</span>
<span id="cb5-35"><a href="#cb5-35"></a><span class="co">#&gt; # Causal Mediation Analysis for Stan Model</span></span>
<span id="cb5-36"><a href="#cb5-36"></a><span class="co">#&gt; </span></span>
<span id="cb5-37"><a href="#cb5-37"></a><span class="co">#&gt;   Treatment: treat</span></span>
<span id="cb5-38"><a href="#cb5-38"></a><span class="co">#&gt;   Mediator : job_seek</span></span>
<span id="cb5-39"><a href="#cb5-39"></a><span class="co">#&gt;   Response : depress2</span></span>
<span id="cb5-40"><a href="#cb5-40"></a><span class="co">#&gt; </span></span>
<span id="cb5-41"><a href="#cb5-41"></a><span class="co">#&gt; Effect                 | Estimate |          95% ETI</span></span>
<span id="cb5-42"><a href="#cb5-42"></a><span class="co">#&gt; ----------------------------------------------------</span></span>
<span id="cb5-43"><a href="#cb5-43"></a><span class="co">#&gt; Direct Effect (ADE)    |   -0.040 | [-0.129,  0.048]</span></span>
<span id="cb5-44"><a href="#cb5-44"></a><span class="co">#&gt; Indirect Effect (ACME) |   -0.018 | [-0.042,  0.006]</span></span>
<span id="cb5-45"><a href="#cb5-45"></a><span class="co">#&gt; Mediator Effect        |   -0.241 | [-0.296, -0.187]</span></span>
<span id="cb5-46"><a href="#cb5-46"></a><span class="co">#&gt; Total Effect           |   -0.057 | [-0.151,  0.033]</span></span>
<span id="cb5-47"><a href="#cb5-47"></a><span class="co">#&gt; </span></span>
<span id="cb5-48"><a href="#cb5-48"></a><span class="co">#&gt; Proportion mediated: 30.59% [-221.09%, 282.26%]</span></span></code></pre></div>
<p>If you want to calculate mean instead of median values from the posterior samples, use the <code>centrality</code>-argument. Furthermore, there is a <code>print()</code>-method, which allows to print more digits.</p>
<div class="sourceCode" id="cb6"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb6-1"><a href="#cb6-1"></a>m &lt;-<span class="st"> </span><span class="kw">mediation</span>(m2, <span class="dt">centrality =</span> <span class="st">&quot;mean&quot;</span>, <span class="dt">ci =</span> <span class="fl">.95</span>)</span>
<span id="cb6-2"><a href="#cb6-2"></a><span class="kw">print</span>(m, <span class="dt">digits =</span> <span class="dv">4</span>)</span>
<span id="cb6-3"><a href="#cb6-3"></a><span class="co">#&gt; # Causal Mediation Analysis for Stan Model</span></span>
<span id="cb6-4"><a href="#cb6-4"></a><span class="co">#&gt; </span></span>
<span id="cb6-5"><a href="#cb6-5"></a><span class="co">#&gt;   Treatment: treat</span></span>
<span id="cb6-6"><a href="#cb6-6"></a><span class="co">#&gt;   Mediator : job_seek</span></span>
<span id="cb6-7"><a href="#cb6-7"></a><span class="co">#&gt;   Response : depress2</span></span>
<span id="cb6-8"><a href="#cb6-8"></a><span class="co">#&gt; </span></span>
<span id="cb6-9"><a href="#cb6-9"></a><span class="co">#&gt; Effect                 | Estimate |            95% ETI</span></span>
<span id="cb6-10"><a href="#cb6-10"></a><span class="co">#&gt; ------------------------------------------------------</span></span>
<span id="cb6-11"><a href="#cb6-11"></a><span class="co">#&gt; Direct Effect (ADE)    |  -0.0395 | [-0.1237,  0.0456]</span></span>
<span id="cb6-12"><a href="#cb6-12"></a><span class="co">#&gt; Indirect Effect (ACME) |  -0.0158 | [-0.0405,  0.0083]</span></span>
<span id="cb6-13"><a href="#cb6-13"></a><span class="co">#&gt; Mediator Effect        |  -0.2401 | [-0.2944, -0.1846]</span></span>
<span id="cb6-14"><a href="#cb6-14"></a><span class="co">#&gt; Total Effect           |  -0.0553 | [-0.1454,  0.0341]</span></span>
<span id="cb6-15"><a href="#cb6-15"></a><span class="co">#&gt; </span></span>
<span id="cb6-16"><a href="#cb6-16"></a><span class="co">#&gt; Proportion mediated: 28.60% [-181.01%, 238.20%]</span></span></code></pre></div>
<p>As you can see, the results are similar to what the <em>mediation</em> package produces for non-Bayesian models.</p>

</body>
</html>
back to top