https://github.com/cran/bayestestR
Revision 40f7c88ddf855896018cf20ec8a7ac5fbd0ea2fb authored by Dominique Makowski on 27 January 2020, 05:30 UTC, committed by cran-robot on 27 January 2020, 05:30 UTC
1 parent d8462ad
Raw File
Tip revision: 40f7c88ddf855896018cf20ec8a7ac5fbd0ea2fb authored by Dominique Makowski on 27 January 2020, 05:30 UTC
version 0.5.1
Tip revision: 40f7c88
credible_interval.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="credible-intervals-ci">Credible Intervals (CI)</h1>
<ul>
<li><a href="#what-is-a-credible-interval">What is a <em>Credible</em> Interval?</a></li>
<li><a href="#why-is-the-default-89">Why is the default 89%?</a></li>
<li><a href="#different-types-of-cis">Different types of CIs</a></li>
<li><a href="#the-support-interval">The Support Interval</a></li>
<li><a href="#references">References</a></li>
</ul>
<p>This vignette can be referred to by citing the package:</p>
<ul>
<li>Makowski, D., Ben-Shachar, M. S., &amp; Lüdecke, D. (2019). <em>bayestestR: Describing Effects and their Uncertainty, Existence and Significance within the Bayesian Framework</em>. Journal of Open Source Software, 4(40), 1541. <a href="https://doi.org/10.21105/joss.01541">https://doi.org/10.21105/joss.01541</a></li>
</ul>
<hr />
<h1 id="what-is-a-credible-interval">What is a <em>Credible</em> Interval?</h1>
<p>Credible intervals are an important concept in Bayesian statistics. Its core purpose is to describe and summarise <strong>the uncertainty</strong> related to your parameters. In this regards, it could appear as quite similar to the frequentist <strong>Confidence Intervals</strong>. However, while their goal is similar, <strong>their statistical definition annd meaning is very different</strong>. Indeed, while the latter is obtained through a complex algorithm full of rarely-tested assumptions and approximations, the credible intervals are fairly straightforward to compute.</p>
<p>As the Bayesian inference returns a distribution of possible effect values (the posterior), the credible interval is just the range containing a particular percentage of probable values. For instance, the 95% credible interval is simply the central portion of the posterior distribution that contains 95% of the values.</p>
<p>Note that this drastically improve the interpretability of the Bayesian interval compared to the frequentist one. Indeed, the Bayesian framework allows us to say <em>“given the observed data, the effect has 95% probability of falling within this range”</em>, while the frequentist less straightforward alternative (the 95% <strong><em>Confidence</em> Interval</strong>) would be “<em>there is a 95% probability that when computing a confidence interval from data of this sort, the effect falls within this range</em>”.</p>
<h1 id="why-is-the-default-89">Why is the default 89%?</h1>
<p>Naturally, when it came about choosing the CI level to report by default, <strong>people started using 95%</strong>, the arbitrary convention used in the <strong>frequentist</strong> world. However, some authors suggested that 95% might not be the most apppropriate for Bayesian posterior distributions, potentially lacking stability if not enough posterior samples are drawn (Kruschke, 2014).</p>
<p>The proposition was to use 90% instead of 95%. However, recently, McElreath (2014, 2018) suggested that if we were to use arbitrary tresholds in the first place, why not use 89% as this value has the additional argument of being a prime number.</p>
<p>Thus, by default, the CIs are computed with 89% intervals (<code>ci = 0.89</code>), deemed to be more stable than, for instance, 95% intervals (Kruschke, 2014). An effective sample size (ESS; see <a href="https://easystats.github.io/bayestestR/reference/diagnostic_posterior.html">here</a>) of at least 10.000 is recommended if 95% intervals should be computed (Kruschke, 2014, p. 183ff). Moreover, 89 is the highest <strong>prime number</strong> that does not exceed the already unstable 95% threshold. What does it have to do with anything? <em>Nothing</em>, but it reminds us of the total arbitrarity of any of these conventions (McElreath, 2018).</p>
<h1 id="different-types-of-cis">Different types of CIs</h1>
<p>The reader might notice that <code>bayestestR</code> provides <strong>two methods</strong> to compute credible intervals, the <strong>Highest Density Interval (HDI)</strong> (<code>hdi()</code>) and the <strong>Equal-tailed Interval (ETI)</strong> (<code>eti()</code>). These methods can also be changed via the <code>method</code> argument of the <code>ci()</code> function. What is the difference? Let’s see:</p>
<div class="sourceCode" id="cb1"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb1-1" title="1"><span class="kw">library</span>(bayestestR)</a>
<a class="sourceLine" id="cb1-2" title="2"><span class="kw">library</span>(dplyr)</a>
<a class="sourceLine" id="cb1-3" title="3"><span class="kw">library</span>(ggplot2)</a>
<a class="sourceLine" id="cb1-4" title="4"></a>
<a class="sourceLine" id="cb1-5" title="5"><span class="co"># Generate a normal distribution</span></a>
<a class="sourceLine" id="cb1-6" title="6">posterior &lt;-<span class="st"> </span><span class="kw">distribution_normal</span>(<span class="dv">1000</span>)</a>
<a class="sourceLine" id="cb1-7" title="7"></a>
<a class="sourceLine" id="cb1-8" title="8"><span class="co"># Compute HDI and ETI</span></a>
<a class="sourceLine" id="cb1-9" title="9">ci_hdi &lt;-<span class="st"> </span><span class="kw">ci</span>(posterior, <span class="dt">method =</span> <span class="st">&quot;HDI&quot;</span>)</a>
<a class="sourceLine" id="cb1-10" title="10">ci_eti &lt;-<span class="st"> </span><span class="kw">ci</span>(posterior, <span class="dt">method =</span> <span class="st">&quot;ETI&quot;</span>)</a>
<a class="sourceLine" id="cb1-11" title="11"></a>
<a class="sourceLine" id="cb1-12" title="12"><span class="co"># Plot the distribution and add the limits of the two CIs</span></a>
<a class="sourceLine" id="cb1-13" title="13">posterior <span class="op">%&gt;%</span><span class="st"> </span></a>
<a class="sourceLine" id="cb1-14" title="14"><span class="st">  </span><span class="kw">estimate_density</span>(<span class="dt">extend=</span><span class="ot">TRUE</span>) <span class="op">%&gt;%</span><span class="st"> </span></a>
<a class="sourceLine" id="cb1-15" title="15"><span class="st">  </span><span class="kw">ggplot</span>(<span class="kw">aes</span>(<span class="dt">x=</span>x, <span class="dt">y=</span>y)) <span class="op">+</span></a>
<a class="sourceLine" id="cb1-16" title="16"><span class="st">  </span><span class="kw">geom_area</span>(<span class="dt">fill=</span><span class="st">&quot;orange&quot;</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb1-17" title="17"><span class="st">  </span><span class="kw">theme_classic</span>() <span class="op">+</span></a>
<a class="sourceLine" id="cb1-18" title="18"><span class="st">  </span><span class="co"># HDI in blue</span></a>
<a class="sourceLine" id="cb1-19" title="19"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_hdi<span class="op">$</span>CI_low, <span class="dt">color=</span><span class="st">&quot;royalblue&quot;</span>, <span class="dt">size=</span><span class="dv">3</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb1-20" title="20"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_hdi<span class="op">$</span>CI_high, <span class="dt">color=</span><span class="st">&quot;royalblue&quot;</span>, <span class="dt">size=</span><span class="dv">3</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb1-21" title="21"><span class="st">  </span><span class="co"># Quantile in red</span></a>
<a class="sourceLine" id="cb1-22" title="22"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_eti<span class="op">$</span>CI_low, <span class="dt">color=</span><span class="st">&quot;red&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb1-23" title="23"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_eti<span class="op">$</span>CI_high, <span class="dt">color=</span><span class="st">&quot;red&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>)</a></code></pre></div>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA8cAAAJACAMAAACXE+S8AAAAqFBMVEUAAAAAADoAAGYAOpAAZrYzMzM6ADo6AGY6kNtBaeFNTU1NTW5NTY5NbqtNjshmAABmtv9uTU1uTW5uTY5ubo5ubqtuq+SOTU2OTW6OTY6Obk2OyP+QOgCQtpCQ2/+rbk2rbm6ryKur5OSr5P+2ZgC2/7a2///Ijk3I///bkDrb///kq27k////AAD/pQD/tmb/yI7/25D/5Kv//7b//8j//9v//+T///+newIPAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAatklEQVR4nO3dDXMcR72FcRHQvUTghBDiJJiLAecSB5uNryKs7//Nrlev+3p6Zk7/u3u6n6dS1JqqrdKc6Z9Xu1KVL26JaO1d1P4CiMgOx0TrD8dE6w/HROsPx0TrD8dE62++Y+Qf9vnnn//735/+p/bX0XMPE7PxmXDsh+P4cKxTKm++vXrx/uHxx9evJjxjzHAcH451QuWW7rsvHv7w7grH58JxfDjWCZU337+9vf767d3j6z/9GcfnwnF8ONYJldffvL+9+e7N9uHHv/3j/vvqy8tLHB+G4/hwrBMqP7x4cvzuJe+Pz4fj+HCsm/R6/OkRjs+H4/hwrJv0/vjd1baXyWcMGo7jw7FOfl798vnzal6Pz4fj+HCsS//8+P4lGcfnw3F8ONbx+1x+OI4Pxzoc++E4PhzrcOyH4/hwrMOxH47jw7EOx344jg/HOhz74Tg+HOtw7Ifj+HCsw7EfjuPDsQ7HfjiOD8c6HPvhOD4c63Dsh+P4cKzDsR+O48OxDsd+OI4Pxzoc+63L8Wav2l/N1HCsw7HfihxvTlT7a5oUjnU49luN41OKV0IZxzoc+63E8VnFa4CMYx2O/VbhWCleAWUc63DstwLHScWtQ8axDsd+7TuewrhtyDjW4divecfTGDctGcc6HPs17niy4pYh41iHY7+2Hc9h3C5kHOtw7Ne043mMm5WMYx2O/Vp2PJtxo5BxrMOxX7uOFyhuFDKOdTj2a9bxMsZNQsaxDsd+3TluEDKOdTj2a9XxYsYNQsaxDsd+bTo2FDcIGcc6HPs16dhj3BxkHOtw7NeiY5dxa5BxrMOxX5+O24KMYx2O/Rp0nIFxW5BxrMOxX3uOszDG8YrCsV9zjvMwbgoyjnU49mvNcS7GLUHGsQ7Hfv06bgcyjnU49mvMcUbG7UDGsQ7Hfm05zsoYxysJx35NOc7LuBnIONbh2K8lx7kZtwIZxzoc+/XtuA3IONbh2K8hxwGM24CMYx2O/XAcH451OPZrx3EIYxyvIBz7NeM4hnETkHGsw7FfK46jGLcAGcc6HPv177g+ZBzrcOzXiONAxvUh41iHY78RHNeGjGMdjv3acBzLGMdth2O/JhwHM8Zx2+HYrwXH0YxrQ8axDsd+DTiOZ1wZMo51OPYbxHFVyDjW4dgPx/HhWIdjv/qOizCuChnHOhz7VXdciDGO2w3HfsM4rggZxzoc+9V2XIwxjpsNx36VHZdjXBEyjnU49hvIcTXIONbh2A/H8eFYh2O/uo6LMq4GGcc6HPtVdVyYMY7bDMd+QzmuBBnHOhz71XRcnDGOmwzHfhUdl2dcCTKOdTj2G8xxFcg41uHYD8fx4ViHY796jqswrgIZxzoc++E4PhzrcOxXzXElxjhuLxz71XJci3ENyDjW4dhvPMflIeNYh2O/So4rMsZxa+HYb0DHxSHjWIdjvxEdl4aMYx2O/eo4rssYx22FY78qjiszLg0Zxzoc++E4PhzrcOxXw3FtxThuKxz7jem4LGQc63DsV8FxbcPbCl4ujlPh2G9Qx0Uh41iHY79RHZeEjGMdjv3KO64N+KFyF4zjRDj2G9ZxQcg41uHYr7jj2nyfKnbFOE6EY7/Sjmvrfa7UFeM4FY79xnVcDjKOdTj2K+y4tt3dCl0yjlPh2G9gx8Ug41iHY7+yjmvL3a/MNeM4FY79RnZcCjKOdTj2G9pxIcg41uHYr6jj2myPKnLVOE6EYz8cx4djHY79Sjqurfa4EleN41Q49ivouDbaUxW4bBynwrHf4I6LQMaxDsd+5RzXFnu6+OvGcSoc+43uuARkHOtw7Ifj+HCsw7Hf8I4LQMaxDsd+xRzX5nq28CvHcSIc+5VyXFvr+aKvHMepcOyH43jIONbh2K+Q49pWVcGXjuNUOPbDcTxkHOtw7FfGcW2puthrx3EqHPvheBMOGcc6HPvheIPjyimVN99evXh/9+jD1dXv3054xpjheFvs1eNYJ1R+fP3q9t0X20fXX799eKSfMWhFHNdmmiz06nGcSKi8+f7tneD7nh7h+LASjmsrTRd59ThOJVRef/P+9ua7Nw9/un89vry8xPFhOL4r8vJxnEio/PDi2fH1V797k37GoBVwXNvolAIvH8epJr8ePz3C8WE4vi/w+nGcaPL749sfXyWfMWjxjmsLnVbc9eM4lfy8+uXDu+Ld77BxfBSOH4obAMeJ0j8/vvuh09UV74/PhuOH4gbAcSJ+n8sPx4/FLYBjHY79wh3X5jm5sAVwnAjHfjh+KmwCHOtw7BftuDbOGUVNgONEOPYLdlzb5pyCJsBxKhz74fi5oA1wnAjHfrGOa8ucV8wGOE6FYz8c7xQzAo4T4dgv1HFtl3MLGQHHqXDsh+PdQlbAcSIc++F4t5AVcJwIx3443itkBhzrcOwX6bg2ygVFzIDjRDj2w/FeETPgOBGO/QId1za5pIAZcJwKx3443i9gBxwnwrFfnOPaIpeVfwccp8KxH44Pyj8EjhPh2C/McW2PS8s+BI5T4dgPx4dlXwLHiXDsh+PDsi+B40Q49sPxUdmnwLEOx35RjmtjNMo9BY4T4dgPx0flngLHiXDsF+S4tkWnzFPgOBWO/WIc16bolXcLHKfCsR+OT5R3DBwnwrFfiOPaDt2yjoHjVDj2w/Gpsq6B40Q49sPxqbKugeNEOPbD8cmyzoFjHY79IhzXRpihnHPgOBGO/XB8spxz4DgRjv0CHNc2mKOMc+A4FY79cHy6jHvgOBGO/fI7ri0wT/n2wHEqHPvh+Ez5BsFxIhz7ZXdc21+usg2C41Q49sPxubItguNEOPbD8bmyLYLjRDj2w/HZsk2CYx2O/XI7ro0vY7kmwXEiHPvh+Hy5NsGxDsd+mR3Xppe1TJvgOBGO/XB8vkyb4DgRjv3yOq4tL295NsFxKhz74ViUZxQcJ8KxX1bHtd3lLssoOE6FYz8cq7KsguNEOPbL6bi2uvzlWAXHqXDsh2NZjllwnAjHfjiW5ZgFx4lw7IdjXY5dcKzDsV9Gx7XJhZRhFxwnwrEfjnUZdsFxIhz75XNcW1xM/i44ToVjPxwn8ofBcSIc+2VzXNtbVFk2xrEKx344TpVjYxyrcOyXy3FtbXHl2BjHKhz74ThZho1xrMKxH46TZdgYxyoc++E4nb8xjlU49svkuDa10PyNcazCsR+O0/kb41iFY788jmtLi83fGMcqHPvheEL2xjhW4dgvi+PazqKzN8axCsd+OJ6SuzGOVTj2y+G4trL43I1xrMKxH44nZW6MYxWO/XA8KXNjHKtw7IfjaXkb41iFY78MjmsTK5K3MY5VOPbD8bS8jXGswrGf77i2sDJ5G+NYhWM/HE/M2hjHKhz72Y5r+yqVtTGOVTj2w/HUnI1xrMKxn+u4tq5yORvjWIVjPxxPztgYxyoc++F4csbGOFbh2A/H01u+MY5VOPYzHdemVbTlG+NYhWM/HE9v+cY4VuHYz3NcW1bZlm+MYxWO/XA8o8Ub41iFYz/LcW1XpVu8MY5VOPbD8ZyWboxjFY79HMe1VZVv6cY4VuHYD8ezWrgxjlU49sPxrBZujGMVjv1wPK9lG+NYhWM/w3FtUlVatjGOVTj2w/HMFm2MYxWO/ZY7rg2qUos2xrEKx344ntmijXGswrHfYse1PdVq0cY4VuHYD8dzW7IxjlU49lvquLamei3ZGMcqHPvheHYLNsaxCsd+Cx3XtlSzBRvjWIVjPxzPb/7GOFbh2A/H85u/MY5VOPbD8YJmb4xjFY79ljmuDalyszfGsQrHfjhe0OyNcazCsd8ix7Ud1W72xjhW4dgPx0uauzGOVTj2W+K4tqL6zd0Yxyoc++F4UTM3xrEKx34LHNc21EIzN8axCsd+OF7WvI1xrMKxH46XNW9jHKtw7Ifjhc3aGMcqpfLm26sX7+8eXX91dfVqwjPGbL7j2oAaadbGOFYJlR9fv7p998X20c13b26v//gm+YxBw/HCZm2MY5VQefP929vrr99+evRhq/nHV8lnDNpsx7X9tNKsjXGsEiqvv3l/90p83/2jy8tLHB+G46XN2RjHKqHyw4sdxx9fv0w/Y9DmOq6tp53mbIxj1cTX45tvHxnj+CgcL27GxjhWTXp/fHv91aun/xvHh810XNtOS83YGMcq+Xn1y4fPq3cZ4/goHC9v+sY4VqV/fvzpJfnd1TY+rz4Tjpc3fWMcq/h9Lj8cG03eGMeqXZW/fPmrv857Bm2b57g2nMaavDGOVfsqf7i4+PW/Zj2DcGw1eWMcqw5V/vLlxcVvZz2DZjmu7aa1Jm+MY9Wxyq3kz/455xmjh2OnqRvjWHWo8qeLi998+v5afHeN48PmOK6tpr2mboxj1Z7K//zl4uIP2wc/ixdkHB+GY6uJG+NYtf95tfp++tQzaNsMx7XNtNjEjXGs4ufHfjj2mrYxjlU49sOx17SNcazCsR+OzSZtjGMVjv2mO64NptEmbYxjFY79cOw2ZWMcq3DsN9lxbS7NNmVjHKtw7Idjtykb41iFY7+pjmtrabcpG+NYhWM/HNtN2BjHKhz7TXRc20rLTdgYxyoc++HYL70xjlU49pvmuLaUtktvjGMVjv1wnKHkxjhW4dgPxxlKboxjFY79Jjmu7aT5UhvjWIVjPxznKLUxjlU49pviuLaS9kttjGMVjv1wnKXExjhW4dhvguPaRtZQYmMcq3Dsh+M86Y1xrMKxX9pxbSHrSG+MYxWO/XCcKbkxjlU49sNxpuTGOFbh2A/HuVIb41iFY7+k49o8VpPaGMcqHPvhOFdqYxyrcOyXclxbx3pSG+NYhWM/HGdLbIxjFY79Eo5r21hTYmMcq3Dsh+N8nd8Yxyoc+2nHtWWsq/Mb41iFYz8cZ+zsxjhW4dgPxzk7tzGOVTj2k45rs1hd5zbGsQrHfjjO2bmNcazCsZ9yXFvF+jq3MY5VOPbDcdbObIxjFY79hOPaJtbYmY1xrMKxH47zdnpjHKtw7HfecW0R6+z0xjhW4dgPx5k7uTGOVTj2w3HuTm2MYxWO/c46rs1htZ3aGMcqHPvhOHenNsaxCsd+5xzX1rDeTm2MYxWO/XCcvRMb41iFY78zjmtbWHMnNsaxCsd+OM7f8cY4VuHY77Tj2hLW3fHGOFbh2A/HAR1tjGMVjv1wHNHhxjhW4djvpOPaDFbf4cY4VuHYD8chHWyMYxWO/U45ro2ggw42xrEKx344DulgYxyrcOx3wnFtAz10sDGOVTj2O3Zcm0Af7W+MYxWO/XAc1N7GOFbh2O/Ice3z30t7G+NYhWM/HEe1uzGOVTj2w3FUuxvjWIVjPxyHtbMxjlU49jt0XPvwd9TOxjhW4djvwHHts99TOxvjWIVjPxzH9bwxjlU49tt3XPvk99XzxjhW4dgPx4E9bYxjFY799hzXPve99bQxjlU49sNxZI8b41iFYz8ch/awMY5VOPbbdVz70HfYw8Y4VuHYD8ehPWyMYxWO/XYc1z7zPfawMY5VOPZ7dlz7yPfZ/cY4VuHYD8fB3eI4FY79nhzXPu+9dovjVDj2w3F0OE6FYz8cR4fjVDj2e3Rc+7R3HI4T4dgPx+HhOBGO/XAcHo4T4dgPx/HhWIdjPxzHh2Mdjv1wHB+OdTj2w3F8ONbh2A/H8eFYh2M/HMeHYx2O/XAcH451OPbDcXw41uHYD8fx4ViHYz8cx4djHY79cBwfjnU49sNxfDjW4dgPx/HhWIdjPxzHh2Mdjv1wHB+OdTi22+A4PhzrcOy2wXGBcKzDsRuOS/ToeJO+IUOGY7MNjkuEYx2Ove7PGI6je3IM5JPh2Ov+jOE4umfHQD4Vjr3uzxiOo8OxDsdWD2cMx9HtOAbyiXDs9HjGcBwdjnU4dno8YziObtcxkI/DsdHTGcNxdDjWKZU33169eP/w+PrrtxOeMVbPZwzH0e05BvJRQuXH169u331x//jD1e9xfNjzGcNxdPuOgXyYUHnz/dvHl+Eff/d3Xo+Pej5jOI4Oxzqh8vqb97c33715+MO948vLSxw/tHPGcBzdgWMgHyRUfnhx7Fg/Y6x2zhiOo8OxbubrsX7GUO2eMRxHh2PdtPfHOD5q74zhOLpDx0DeT35e/fLp82ocH7R/xnAc3ZFjIO+V/vnxvWAc77d/xnAc3bFjIO/G73Mta/+M4Tg6HOtwvKiDM4bj6E44BvJOOF7S4RnDcXSnHAP5ORwv6fCM4Tg6HOtwvKCjM4bj6HCsw/H8js8YjqM76RjIT+F4fsdnDMfRnXYM5MdwPL/jM4bj6HCsw/HsTpwxHEd3xjGQH8Lx3E6dMRxHd84xkO/D8cxOnjEcR4djHY5ndvKM4Ti6s46BfBeO53X6jOE4OhzrcDyrM2cMx9GddwzkbTie1ZkzhuPohGMg3+J4XufOGI6jU46BjON5nTtjOI4Oxzocz+jsGcNxdDjW4Xh6588YjqOTjoGM4+mJM4bj6LRjION4cuKM4Ti6hOPhIeN4cuKM4Ti6lOPRIeN4auqM4Tg6HOtwPDF5xnAcHY51OJ6WPmM4ji7peHDIOJ5U4ozhOLq047Eh43hSiTOG4+gmOB4aMo6nlDpjOI5uiuORIeN4SqkzhuPocKzD8YSSZwzH0eFYh+N06TOG4+gmOR4YMo6TTThjOI5umuNxIeM42YQzhuPoJjoeFjKOU005YziObqrjUSHjONGkM4bj6CY7HhQyjnXTzhiOo8OxDseyiWcMx9FNdzwmZByrpp4xHEc3w/GQkHGsmnrGcBzdHMcjQsaxaPIZw3F0sxwPCBnH55t+xnAcHY51OD7f9DOG4+jmOR4PMo7PNuOM4Tg6HOtwfK45ZwzH0c10PBxkHJ9rzhnDcXRzHY8GGcenm3fGcBzdbMeDQcbxyWaeMRxHN9/xWJBxfKq5ZwzH0S1wPBRkHJ9q7hnDcXRLHI8EGccnmn3GcBzdIscDQcbxcfPPGI6jw7EOx0ctOGM4jm6Z43Eg4/iwJWcMx9EtdDwMZBwftuSM4Ti6pY5HgYzjgxadMRxHt9jxIJBxvN+yM4bj6JY7HgMyjvdaeMZwHB2OdTjeafEZw3F0huMhJON4p8VnDMfRWY4HgIzj55afMRxH5znuHzKOnzLOGI6jw7EOx485ZwzH0ZmOu4eM44esM4bj6FzHvUvG8X3eGcNxdL7jviHj+C7zjOE4OhzrcHxrKsZxiTI47loyjm9xvIKyOO4YMo5txjguUB7H/ULGcY4zhuPoMjnuFvLwjrOcMRxHl8txr5BxnOOM4Ti6bI47hTy4Y/9YbHBconyO+4Q8tuMMp2KD4xJldNwl5KEd5zgUGxyXKKfjHiWP7DjPmcBxifI67g/yuI4znYgNjkuU2XF3kId1nOtAbHBcotyOe4M8quNs52GD4xJld9yZ5DEdZzwNGxyXKMBxV5CHdJzzMGxwXKIIxz1BHtFx1rOwwXGJQhx3JHk8x5lPwgbHJQpy3A3k4RznPggbHJcoynEvkEdznP0cbHBcojDHnUgey3HAKdjguESBjruAPJLjiDOwDcfxRTruQfJAjmOOwAbHJYp1vH7JwziOOgAbHJco2vHaIQ/iOOz2b8NxfOGOVy55DMeBd3+D4xIVcLxqySM4Dr33GxyXqIjjFUPu33Hsnd+G4/jKON6slnLvjsPv+wbHJSrmeKWQO3ccf9s3OC5ROcfrlNyz4xL3fBuO4yvpeI2SO3Zc5pZvcFyiso7XJ7lXx6Xu9zYcx1fa8dok9+m43N3ehuP4yjverIpyh46L3uptOI6viuMVSe7OceEbvQ3H8VVyvFkL5b4cl7/L23AcXz3H65Dck+Ma93gbjuOr6XgNkrtxXOcG34Xj+Oo63lb7gOv6cFzv7m7DcXz1HW+aprx+x1Xv7F04jq8Jx5t2Ka/cce3beheO42vF8aZRyit2XPt+PoXj+BpyvK322T9qrY5r38jdcBxfY4631Saw1wod175/R+E4vgYdb6tt4amVOa59306G4/gadXxXTRCPrcZx7ZslwnF8LTu+q46Kp9bguPYtSoXj+Jp3fF9xG4+17bj2bZkWjuNbieP7Cgp5rFHHte/ErHAc36ocP1QCymNK5c23Vy/eHzwKdlx7+kXhOL41On4sEsxjQuXH169u332x/0g/Y2G1d3bDcXxrdvxcdjpPCZU337+9vf767d4j/Yyp1Z4zcziOrw/Hh/mUHhMqr795f3vz3Zu9R5eXl9V/D6S5Hh3X/jp67tFx7a+j1YTKDy8e9T4/0s8YNBzHh2PdzNdj/YxBw3F8ONZVeX/cWTiOD8c6+Xn1y6fPq19Gfl699nAcH4516Z8fb1+Ii/38eJXhOD4c6xr9fa5VheP4cKzDsR+O48OxDsd+OI4Pxzoc++E4PhzrcOyH4/hwrMOxH47jw7EOx344jg/HOhz74Tg+HOtw7Ifj+HCsw7EfjuPDsQ7HfjiOD8c6HPvhOD4c63Dsh+P4cKzDsR+O48OxDsd+OI4Pxzoc++E4PhzrcOyH4/hwrMOxH47jw7FugePHLi+6jstbc6Nc3nLHT10uf+oa4vLW3GiXh+NzcXlrbrTLw/G5uLw1N9rl8akV0frDMdH6wzHR+sMx0frDMdH6cxx/fP0q29fRXNdfXV11e3m7/15Xh3V977Ydy3Mcv+t4rO0/93z9xze1v4yYtsfg6d/P7K+u791dx/IMx9d/+nO/jj9sj/mPnV7f3r9n3V9d37ttJ+Qtd/zxb//o+fvq2/u/17vs+pv3/V7cfT1f3il5yx2/e9n1++P7f729zz686N1xv/fu9rS8RY5/vLr64tPf6b063l7e9rOgbo9C96/HHd+7u7uX8XOud1fb+p3r+qs+/5La1vn7467v3Rl5/NzpdF0fhe13nR1/Xt31vbsr78+denZ8/5der9fX98+P+7532/I6JqI2wjHR+sMx0frDMdH6wzHR+sMx0frDMdH6wzHR+sMx0frDMdH6w/GA/fSrv97+8uVva38ZlC8cj9gPv9n+R/2E4xH7v//6n//+Z+0vgjKG4yH76eIPtb8EyhmOh+yHC76t7iocj9jPn/3vl7wg9xSOB+yXT4h//ow3yB2F4wH74df/uv3PX/jOuqNwTLT+cEy0/nBMtP5wTLT+cEy0/nBMtP5wTLT+cEy0/nBMtP5wTLT+/h8uq8cwBeh6SwAAAABJRU5ErkJggg==" /><!-- --></p>
<blockquote>
<p><strong>These are exactly the same…</strong></p>
</blockquote>
<p>But is it also the case for other types of distributions?</p>
<div class="sourceCode" id="cb2"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb2-1" title="1"><span class="kw">library</span>(bayestestR)</a>
<a class="sourceLine" id="cb2-2" title="2"><span class="kw">library</span>(dplyr)</a>
<a class="sourceLine" id="cb2-3" title="3"><span class="kw">library</span>(ggplot2)</a>
<a class="sourceLine" id="cb2-4" title="4"></a>
<a class="sourceLine" id="cb2-5" title="5"><span class="co"># Generate a beta distribution</span></a>
<a class="sourceLine" id="cb2-6" title="6">posterior &lt;-<span class="st"> </span><span class="kw">distribution_beta</span>(<span class="dv">1000</span>, <span class="dv">6</span>, <span class="dv">2</span>)</a>
<a class="sourceLine" id="cb2-7" title="7"></a>
<a class="sourceLine" id="cb2-8" title="8"><span class="co"># Compute HDI and Quantile CI</span></a>
<a class="sourceLine" id="cb2-9" title="9">ci_hdi &lt;-<span class="st"> </span><span class="kw">ci</span>(posterior, <span class="dt">method =</span> <span class="st">&quot;HDI&quot;</span>)</a>
<a class="sourceLine" id="cb2-10" title="10">ci_eti &lt;-<span class="st"> </span><span class="kw">ci</span>(posterior, <span class="dt">method =</span> <span class="st">&quot;ETI&quot;</span>)</a>
<a class="sourceLine" id="cb2-11" title="11"></a>
<a class="sourceLine" id="cb2-12" title="12"><span class="co"># Plot the distribution and add the limits of the two CIs</span></a>
<a class="sourceLine" id="cb2-13" title="13">posterior <span class="op">%&gt;%</span><span class="st"> </span></a>
<a class="sourceLine" id="cb2-14" title="14"><span class="st">  </span><span class="kw">estimate_density</span>(<span class="dt">extend=</span><span class="ot">TRUE</span>) <span class="op">%&gt;%</span><span class="st"> </span></a>
<a class="sourceLine" id="cb2-15" title="15"><span class="st">  </span><span class="kw">ggplot</span>(<span class="kw">aes</span>(<span class="dt">x=</span>x, <span class="dt">y=</span>y)) <span class="op">+</span></a>
<a class="sourceLine" id="cb2-16" title="16"><span class="st">  </span><span class="kw">geom_area</span>(<span class="dt">fill=</span><span class="st">&quot;orange&quot;</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb2-17" title="17"><span class="st">  </span><span class="kw">theme_classic</span>() <span class="op">+</span></a>
<a class="sourceLine" id="cb2-18" title="18"><span class="st">  </span><span class="co"># HDI in blue</span></a>
<a class="sourceLine" id="cb2-19" title="19"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_hdi<span class="op">$</span>CI_low, <span class="dt">color=</span><span class="st">&quot;royalblue&quot;</span>, <span class="dt">size=</span><span class="dv">3</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb2-20" title="20"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_hdi<span class="op">$</span>CI_high, <span class="dt">color=</span><span class="st">&quot;royalblue&quot;</span>, <span class="dt">size=</span><span class="dv">3</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb2-21" title="21"><span class="st">  </span><span class="co"># Quantile in red</span></a>
<a class="sourceLine" id="cb2-22" title="22"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_eti<span class="op">$</span>CI_low, <span class="dt">color=</span><span class="st">&quot;red&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb2-23" title="23"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>ci_eti<span class="op">$</span>CI_high, <span class="dt">color=</span><span class="st">&quot;red&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>)</a></code></pre></div>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA8cAAAJACAMAAACXE+S8AAAAqFBMVEUAAAAAADoAAGYAOpAAZrYzMzM6ADo6AGY6kNtBaeFNTU1NTW5NTY5NbqtNjshmAABmtv9uTU1uTW5uTY5ubqtuq8huq+SOTU2OTW6OTY6ObquOyP+QOgCQtpCQ2/+rbk2rbm6rjk2ryKur5P+2ZgC2/7a2///Ijk3I///bkDrb///kq27k////AAD/pQD/tmb/yI7/25D/5Kv//7b//8j//9v//+T///9v35IHAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAWqklEQVR4nO3d63bcxpWGYcYOZzLMwVHisT1RJkrGjg+JRWtoRbz/O4u62WSfAHxVQNXeu6reZ3ll8Qe8ggK+12yJP3jzCKB1N943AGAzOgbaR8dA++gYaB8dA+2jY6B9qR3TexVv3+b+G/91UONuMGH+FUV6FXTsio7Do2NIdBweHUOi4/DoGBIdh0fHkOg4PDqGRMfh0TEkOg6PjiHRcXh0DImOw6NjSHQcHh1DouPw6BgSHYdHx5DoODw6hkTH4dExJDoOj44h0XF4dAyJjsOjY0h0HB4dQ6Lj8OgYEh2HR8eQ6Dg8OoZEx+HRMSQ6Do+OIdFxrvsL1f8P6RgSHee4bNikZTqGRMep5hqunjIdQ6LjNKLimiXTMSQ6TpAQccWS6RgSHUupFdcqmY4h0bGQU3GdkukYEh0vy824Qsh0DImOl+RXXCFkOoZEx/NWVVy+ZDqGRMezVmdcOGQ6hkTHMzZUXDhkOoZEx5O2VVw2ZDqGRMdTNmdcsmQ6hkTHE0pkXC5kOoZEx1fKVEzH265DFjq+VCrjYiXTMSQ6vlAw40Ih0zEkOj5TtOJCIdMxJDo+VTpjOl59HbLQ8aniHZcImY4h0fFR+YqLhEzHkOj4RZWMC4RMx5Do+FmljLeHTMeQ6PigWsabQ6ZjSHT8pGLGdLzmOmSh452aFW8OmY4h0fFO5Y63hUzHkOj4sX7G20KmY0h0bJHxppDpGBIdW1RMx9nXIcvwHdtkvCVkOoZEx3RcBB27Gr1jq4w3hEzHkAbv2C7j9SHTMaShO7aseH3IdAyJjsOHTMeQRu7YOmM6zrgOWQbu2DxjOs64DlnG7dg+45Uh0zGkYTv2yHhdyHQMadSOfTKmYzqug47Dh0zHkMbs2KviVSHTMaQhO3bMmI6r3sWwRuzYM2M6rnoXw6Lj8CHTMaQBO/bNOD9kOoY0XsfeGdMxyhuuY++K77NDpmNIo3Xs3fBe3i3TMaTBOvYu+EnePdMxJDr2kHXPdAxprI69+32WddN0DGmojr3zfZF113QMaaSOves9kXPbdAxpoI692z2Vc990DImOfWTcNx1DGqdj73IvpN84HUMapmPvbi+l3zkdQxqlY+9sryXfOh1DGqRj72gnJN87HUMao2PvZqck3zwdQ6JjN6k3T8eQhujYu9hpqXdPx5BG6Ng72DmJt0/HkAbo2DvXWYn3T8eQ+u/Yu9YFaQegY0h07CjtAHQMiY49JR2AjiF137F3qouSTkDHkHrv2LtUIeUIdAyp8469O1VSzkDHkPru2DtTKeUQdAyp6469K02QcAo6hkTHvhJOQceQeu7Yu9Ek+hh0DKnjjr0LTaPPQceQ+u3YO9BU8iB0DKnbjr3zTCZPQseQ6NidOgkdQ6Jjd+okdAyp146948ygjkLHkDrt2LvNHOosdAypz46908yiDkPHkLrs2LvMTOI0dAypx469u8wljkPHkOg4gOXj0DGkDjv2rjLf8nnoGBIdR7B4HjqG1F/H3k2usXggOobUXcfeSa6zdCI6htRbx95BrrR0JDqGRMchLB2JjiF11rF3j2stnYmOIfXVsXeO6y0cio4h0XEMC4eiY0hddewd4wYLp6JjSD117N3iJvPHomNIHXXsXeI28+eiY0h0HMXsuegYUj8de3e41ezB6BhSNx17Z7jd3MnoGBIdhzF3MjqG1EvH3hEWMHc0OobUScfeDRYxczY6hkTHccycjY4h9dGxd4FlzByOjiF10bF3gIXMnI6OIdFxINOno2NIPXTsnV8x08ejY0gddOxdXznT56NjSO137B1fSZMHpGNIdBzJ5AHpGFLzHXunV9bUCekYEh2HMnVCOobUesfe4RU2dUQ6htR4x97dFTdxRjqG1HbH3tWVN3FIOoZEx7FMHJKOITXdsXd0FUycko4h0XEsE6ekY0gtd+zdXBXXx6RjSA137F1cHdfnpGNI7XbsHVwtVwelY0jNduydWzVXJ6VjSHQczdVJ2+/44dXd3euE67AaHYdzedLmO37/1dePD3/4Wl6H9Vrt2Du2ii6P2nzH73738X++ey2vw3qNduzdWk2XZ22+453d9+THx9vbWzquos2OvVOr6+KwPXT84c0XSddhpSY79g6tsovTdtDx+y+/SLoOa9FxQOenbb/jh1evk67DanQc0Plpm+/4NGM6rqPFjr0zq+78uM13/M+7Hf6+uqYGO/aurL7z8zbf8arrkKW9jr0jM3B+YDqGRMcBnR+YjiE117F3YybOTkzHkFrr2LswG2dHpmNIdBzR2ZHpGFJjHXsHZuX0zHQMqa2OvfMyc3poOoZExzGdHJqOITXVsXdchk5OTceQWurYuy1LJ8emY0h0HNTx2HQMqaGOvcuydTw3HUOi46CO56ZjSO107B2WsePB6RhSMx17d2XteHI6htRKx95ZmTsenY4hNdKxd1UOXs5Ox5DoOKqXs9MxJDqO6uXsdAypjY69m3LxfHg6htREx95F+Xg+PR1DaqFj76CcPB+fjiHRcVyH49MxpAY69s7JzeH8dAyJjuM6nJ+OIcXv2LsmR08PgI4hhe/YuyVPT0+AjiHRcWBPT4COIUXv2DslX/tHQMeQgnfsHZKz/TOgY0ixO/buyNv+IdAxJDqObP8Q6BgSHUe2fwh0DCl0x94Zuds/BTqGFLlj74oC2D0GOoYUuGPvhiLYPQc6hhS3Y++EQtg9CDqGRMexPdIxEtBxbI90jARhO/YOKIhHOkaCqB179xMGHSNB0I6964mDjpEgZsfe8QRCx0hAx9HRMTQ6jo6OoYXs2DudUOgYWsSOvcuJhY6hBezYO5xg6BgaHUdHx9DidezdTTR0DC1cx97ZxEPHkOg4PDqGFK1j72gComNIdBweHUMK1rF3MxHRMaRYHXsnExIdQwrVsXcxMdExJDoOj44h0XF4dAwpUsfewQRFx5ACdezdS1Rv395bv4oV6NhVnI69cwmLjiGF6di7lrjoGFKUjr1jCexjxzMh0zEO6Dg8OoZEx+HRMaQgHXu3EhkdQ4rRsXcqodExpBAde5cSGx1DouPw6BhShI69Qwlu1/F0yHSMAzoOj44hBejYu5Po6BiSf8femYS373gyZDrGgXvH3pXER8eQ6Dg8Oobk3bF3JA146ngqZDrGgXPHzok0gY4h0XF4dAzJt2PnQtpAx5BcO3YOpBGHjidCpmMc0HF4dAzJs2PfPJpBx5AcO/atox10DImOw6NjSH4d+8bREDqG5NaxbxstoWNIdBzec8fXIdMxDrw6di2jLXQMyalj1zAaQ8eQfDp27aI1Lx1fhUzHOKDj8OgYkkvHnlW0h44h0XF4x44vQ6ZjHHh07NhEi+gYkkPHjkk0iY4h2XfsWESb6BgSHYd30vFFyHSMA/OO/XpoVXsd//zZL/6Sch2KoePw2uv48fGbm5tP/5FwHQqx7tgvh2addnwectyOd9+Tb25+lXAdijDu2C2GhrXZ8VPJn3yrr0MBth27tdCyNjv+4ebmlx8/X199uqbjKkw7dkuhaQ12/K8/39z8fvfFT1ffkOm4CjoOr72Of/7s+vP01HUoxrJjtxLa1l7HJa5DFjoOj44hGXbsFkLjzjo+C5mOcWDXsVcGzaNjSGYde1XQPjqGRMfh0TEkq469IujAecenIdMxDow6dkqgC3QMyaZjpwL6QMeQ6Di8i45PQqZjHJh07LP/XtAxJIuOfebfDTqGRMfh0TEkg4591t+Py46PIdMxDup37LL9ntAxpOodu0y/K3QMqXbHLsvvCx1DouPw6BgSHYdHx5Aqd+wy/M7QMaS6HbvsvjdXHb+ETMc4qNqxx+r7Q8eQanbsMfoO0TGkih17bL5HdAyJjsO77vg5ZDrGQb2OHRbfJzqGRMfh0TGkah07DL5TEx3fZ70KE3TsqlbH9nPvFh1DqtSx/dr7RceQ6Dg8OoZUp2P7sXdsquP75FdhhY5dVenYfOpdo2NINTo2X3rf6BhShY7Nh965yY7vk16FITp2Vb5j65l3j44h0XF4dAyJjsOjY0jFO7Zeef/oGFLpjq1HPgA6hlS4Y+uNj4COIZXt2HriQ6BjSHQc3nTH9+JVWKNjV0U7Nh74IOgYUsmOjfc9CjqGVLBj43kPY6bjezrGi3Id2457IHQMqVjHttseCR1DouPw6BhSqY5tpz2UuY7v6RjPCnVsOuzB0DGkMh2b7no0dAypSMemsx7ObMf3dIwDOg6PjiHRcXh0DKlEx5ajHhAdQyrQseWmR0THkLZ3bDnpIdExpM0dWy56TD10/PDH75Ouw0pbO7Yc9KA66Pjd3W/ouCo6Dq/9jr/79f/x/biujR1b7nlU7Xf88rn69vaWjqvY1rHlnIfVUcfyOqy0qWPLNY+LjiFt6dhyzAOjY0h0HB4dQ9rQseWWR0bHkNZ3bDnlofXQcf51yELH4dExJDoOj44h0XF4dAyJjsOjY0h0HB4dQ6Lj8OgYEh2HR8eQ6Dg8OoZEx+HRMSQ6Do+OIeV2fByP5ZSHRseQMju+p2NzdAwpr+PT8dgNeXB0DImOw6NjSFkdn43HbMejo2NIOR2fj8dqxsOjY0gZHV+Mx2jFoGNI6R1fjsdmxKBjaHQcHh1DSu74ajwmGwYdI0Fqx9fjsZgw7ukYCRI7nhiPwYKxQ8eQ0jqeGk/9AWNPd3xfeyUJ6NgVHYdHx5BSOp4ZD2zQMaSEjufGAxt0DEl3PDse2KBjSLLj+fHABh1DouPw6BiS6nhhPLCR0HGAkOnYleh4cTwwQceQljteHg9M0DGkxY7FeGAipWP/kOnY1VLHcjywQMeQFjrW44EFOoY033HCeGCBjiHNdjw7Kzo2ltSxe8h07Gqu4/lZ0bExOoZEx+HRMaTpjpdmRcfG0jr2DpmOXU12vDgrOjZGx5CmOl6eFR0bo2NIdBxeYsfOIdOxq4mOxazo2BgdQ7ruWM2Kjo3RMaTLjvWs6NgYHUO66DhhVnRsjI4h0XF4dAzpvOOUWdGxMTqGdNpx2qzo2Fhqx74h07Grk44TZ0XHxugYEh2HR8eQjh2nzoqOjSV37BoyHbt67jh9VnRsjI4hHTrOmBUdG6NjSE8d58yKjo2ld+wZMh27ouPw6BjSruO8WdGxMTqG9LHjzFnRsbGMjh1DpmNX8yNJHw+qomNIdBweHUOi4/DoGMrCSNLHg6roGMsWR5I+HlRFx1i2OJL08aCqnI79QqZjN8sjyRgPaqJjLFAjyRgPaqJjzJMjyRgPasrq2C1kOvYhR5IzHlREx5iTMJKc8aAiOsaMlJHkjAcV5XXsFTIdm0sbSdZ4UA8dY1LaSLLGg3oyO3YKmY6tJY4kbzyoho5xLXkkeeNBNXSMK+kjyRsPqsnt2CdkOraUMZLM8aAWOsaFnJFkjge10DHO5I0kczyohY5xKnMkmeNBLXSME7kjyRwPaqFjHGWPJHM8qCW7Y5eQ6djCipHkjgeV0DGerRhJ7nhQSX7HHiHTcX2rRpI9HtRBx9hbN5Ls8aAOOsaj+i1sdBzeio4dQqbjulaPJH88qIKOIX+ZIh2HR8fD2zKS/PGgijUd24dMxxVtGcmK8aAGOh7ctpGsGA9qoOOhbR3JivGgBjoe2eaRrBgPaqDjgW0fyYrxoIZVHZuHTMcVlBjJmvGgAjoeVZGRrBkPKljXsXXIdFxcmZGsGg/Ko+MhlRrJqvGgPDoeULmRrBoPylvZsXHIdFxSwZGsGw+Ko+PhlBzJuvGguLUd24ZMx8WUHcnK8aA0Oh5L4ZGsHA9KW92xach0XEbxkawdDwqj43FUGMna8aAwOh5GjZGsHQ8Ko+NB1BnJ2vGgMDoeQ6WRrB0PClvfsWXIdLxNtZHMoWNjdNy/iiOZQ8fGNnRsGDIdb1BzJHPo2Bgdd67uSObQsbEtHduFTMcr1R7JHDo2Rsc9qz6SOXRsbFPHZiHT8RoGI5lDx8bouFsWI5lDx8bouFM2I5lDx8bouEdmI5lDx8a2dWwVMh1nsRvJHDo2RsfdKb4ROo5vY8dGIdNxstID2aHj8Oi4K6Xn8YSOw9vasU3IdJyk8DZe0HF4mzs2CZmOE5Qdxik6Do+OO1F2FufoOLztHVuETMdC0U1coePw6LgDRRcxgY7Do+PmFd3DJDoOj47bVnQMc+g4vAIdG4RMx9NKLmEBHYdXouP6IdPxlIIzWEbH4dFxqwqOQKHj8Ip0XD1kOr5QbgAp6Di8Mh3XDpmOzxR7+4noODw6bk2xV5+OjsOj47YUe/E56Di8Qh1XDpmOn5R665noODw6bkWpN74CHYdXquO6IdNxode9Dh2HV6zjqiEP3nGZV70eHYdXruOaIY/ccZHXvA0dh1ew44ohD9txiVe8HR2HR8dxlXi/RdBxeHQcU4l3Wwwdh1ey43ohD9ZxgfdaFB2HV7TjaiGP1PH2d1ocHYdXtuNaIY/S8ebXWQcdh1e440ohj9Dx1jdZER2HR8chbH2NddFxeKU7rhNy1x1vfIMG6Di84h1XCbnXjre9OzN0HF75jmuU3GPHW96aMToOj44dbHlhHug4vBodlw+5n443vCo/dBxelY6Ll9xFx6vfkTs6Dq9Sx4VDbrzjtS8nCjoOr1bHZUNutuOVDy8YOg6vWsdFS26w49WPLSA6Dq9ixwVLbqjjDY8rLDoOr2rHxUJuoOMtjyk6Og6vbsf3hVIO3PHGx9MEOg6vesdFQo7X8dan0hQ6Dq9+xztbq1nq8/2Xd7/9MeG6Eko8iwbRcXg2Hd9vTHmhzw9vXj/+83f6uvUKPYCG0XF4Zh3vrU1poc/3//P948Mfv5fXZSh54i7QcXi2HT/JL2uhz4fPf3x8/9XXH7+6vb0N8HOnHr19m/tvPI+nxt1gwvwrivQqFvp899vnjpevw3p0HF7zHR+/Hy9fh/XoOLzmO67w52NcoOPwmu/4w5svKv99Neg4vOY7Nv358aDoOLz2O15zHbLQcXh0DImOw6NjSHQcHh1DouPw6BgSHYdHx5DoODw6hkTH4dExJDoOj44h0XF4dAyJjsOjY0h0HB4dQ6Lj8OgYEh2HR8eQ6Dg8OoZEx+HRMSQ6Do+OIdFxeHQMiY7D66zjVtx630BEPJQJXTyU3I6bcet9AxHxUCZ09VDoeAQ8lAldPRQ6HgEPZUJXD6W7joEB0THQPjoG2kfHQPvoGGhfNx0ffxnVu7u733zvfDdBHB/Khzd3v/7a+W6COD6Uh1fdLKWXjj+8eX345ZC7X/X68msix3Z8KI/fvd7/XnqcPJT3X378qpOH0kvHZ7+s+eSroR0fyu4r7B0fysPnPz6+/6qPTym9dHz+Tvh+vHd8KA+f/43P1U9OHwodR7P70Pj8Th5eMdm940N5ePV6P1ucPJT95+pOptJLx+f/be3lv7Ib9fmtZ6OTR/Hw6u6//9rHQ+ml4/M/H+/+Wgenfz7+Ex0fnC+ll7846KXjD2++OPyp+PQT9uCOD2X3HzY+V+8dH8puJb38TUovHR9+Krj/odMdf6VzcHwoH7/q5UelWx0fyru7u05+7NRPx8DA6BhoHx0D7aNjoH10DLSPjoH20THQPjoG2kfHQPvoGGgfHY/lh1/85fHnz37lfRsojI4H880vd/+gM3Q8mP//j//9z2+9bwKl0fFofrj5vfctoDg6Hs03N3ys7g8dD+anT/7+Gd+Qu0PHY/n5Y8Q/fcIfkHtDx2P55tN/PP7rz3yy7g0dA+2jY6B9dAy0j46B9tEx0D46BtpHx0D76BhoHx0D7aNjoH3/BkOEnErEzZBFAAAAAElFTkSuQmCC" /><!-- --></p>
<blockquote>
<p><strong>The difference is strong with this one.</strong></p>
</blockquote>
<p>Contrary to the <strong>HDI</strong>, for which all points within the interval have a higher probability density than points outside the interval, the <strong>ETI</strong> is <strong>equal-tailed</strong>. This means that a 90% interval has 5% of the distribution on either side of its limits. It indicates the 5th percentile and the 95th percentile. In symmetric distributions, the two methods of computing credible intervals, the ETI and the HDI, return similar results.</p>
<p>This is not the case for skewed distributions. Indeed, it is possible that parameter values in the ETI have lower credibility (are less probable) than parameter values outside the ETI. This property seems undesirable as a summary of the credible values in a distribution.</p>
<p>On the other hand, the ETI range does change when transformations are applied to the distribution (for instance, for a log odds scale to probabilities): the lower and higher bounds of the transformed distribution will correspond to the transformed lower and higher bounds of the original distribution. On the contrary, applying transformations to the distribution will change the resulting HDI.</p>
<h1 id="the-support-interval">The Support Interval</h1>
<p>Unlike the HDI and the ETI, which look at the posterior distribution, the <strong>Support Interval (SI)</strong> provides information regarding the change in the credability of values from the prior to the posterior - in other words, it indicates which values of a parameter are have gained support by the observed data by some factor greater or equal to <em>k</em> (Wagenmakers, Gronau, Dablander, &amp; Etz, 2018).</p>
<div class="sourceCode" id="cb3"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb3-1" title="1">prior &lt;-<span class="st"> </span><span class="kw">distribution_normal</span>(<span class="dv">1000</span>, <span class="dt">mean =</span> <span class="dv">0</span>, <span class="dt">sd =</span> <span class="dv">1</span>)</a>
<a class="sourceLine" id="cb3-2" title="2">posterior &lt;-<span class="st"> </span><span class="kw">distribution_normal</span>(<span class="dv">1000</span>, <span class="dt">mean =</span> <span class="fl">.5</span>, <span class="dt">sd =</span> <span class="fl">.3</span>)</a>
<a class="sourceLine" id="cb3-3" title="3"></a>
<a class="sourceLine" id="cb3-4" title="4">si_<span class="dv">1</span> &lt;-<span class="st"> </span><span class="kw">si</span>(posterior, prior, <span class="dt">BF =</span> <span class="dv">1</span>)</a>
<a class="sourceLine" id="cb3-5" title="5">si_<span class="dv">3</span> &lt;-<span class="st"> </span><span class="kw">si</span>(posterior, prior, <span class="dt">BF =</span> <span class="dv">3</span>)</a>
<a class="sourceLine" id="cb3-6" title="6"></a>
<a class="sourceLine" id="cb3-7" title="7"><span class="kw">ggplot</span>(<span class="dt">mapping =</span> <span class="kw">aes</span>(<span class="dt">x=</span>x, <span class="dt">y=</span>y)) <span class="op">+</span></a>
<a class="sourceLine" id="cb3-8" title="8"><span class="st">  </span><span class="kw">theme_classic</span>() <span class="op">+</span></a>
<a class="sourceLine" id="cb3-9" title="9"><span class="st">  </span><span class="co"># The posterior</span></a>
<a class="sourceLine" id="cb3-10" title="10"><span class="st">  </span><span class="kw">geom_area</span>(<span class="dt">fill =</span> <span class="st">&quot;orange&quot;</span>,</a>
<a class="sourceLine" id="cb3-11" title="11">            <span class="dt">data =</span> <span class="kw">estimate_density</span>(posterior, <span class="dt">extend =</span> <span class="ot">TRUE</span>)) <span class="op">+</span></a>
<a class="sourceLine" id="cb3-12" title="12"><span class="st">  </span><span class="co"># The prior</span></a>
<a class="sourceLine" id="cb3-13" title="13"><span class="st">  </span><span class="kw">geom_area</span>(<span class="dt">color =</span> <span class="st">&quot;black&quot;</span>, <span class="dt">fill =</span> <span class="ot">NA</span>, <span class="dt">size =</span> <span class="dv">1</span>, <span class="dt">linetype =</span> <span class="st">&quot;dashed&quot;</span>,</a>
<a class="sourceLine" id="cb3-14" title="14">            <span class="dt">data =</span> <span class="kw">estimate_density</span>(prior, <span class="dt">extend =</span> <span class="ot">TRUE</span>)) <span class="op">+</span></a>
<a class="sourceLine" id="cb3-15" title="15"><span class="st">  </span><span class="co"># BF = 1 SI in blue</span></a>
<a class="sourceLine" id="cb3-16" title="16"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>si_<span class="dv">1</span><span class="op">$</span>CI_low, <span class="dt">color=</span><span class="st">&quot;royalblue&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb3-17" title="17"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>si_<span class="dv">1</span><span class="op">$</span>CI_high, <span class="dt">color=</span><span class="st">&quot;royalblue&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb3-18" title="18"><span class="st">  </span><span class="co"># BF = 3 SI in red</span></a>
<a class="sourceLine" id="cb3-19" title="19"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>si_<span class="dv">3</span><span class="op">$</span>CI_low, <span class="dt">color=</span><span class="st">&quot;red&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>) <span class="op">+</span></a>
<a class="sourceLine" id="cb3-20" title="20"><span class="st">  </span><span class="kw">geom_vline</span>(<span class="dt">xintercept=</span>si_<span class="dv">3</span><span class="op">$</span>CI_high, <span class="dt">color=</span><span class="st">&quot;red&quot;</span>, <span class="dt">size=</span><span class="dv">1</span>)</a></code></pre></div>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA8cAAAJACAMAAACXE+S8AAAAq1BMVEUAAAAAADoAAGYAOpAAZrYzMzM6ADo6AGY6kNtBaeFNTU1NTW5NTY5NbqtNjshmAABmtv9uTU1uTW5uTY5ubo5ubqtuq+SOTU2OTW6OTY6Obk2OyP+QOgCQtpCQ2/+rbk2rbm6rbo6ryKur5OSr5P+2ZgC2/7a2///Ijk3I///bkDrb///kq27k////AAD/pQD/tmb/yI7/25D/5Kv//7b//8j//9v//+T///+r9GFFAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAU8klEQVR4nO3dC3vcRgGF4W3BQEzTFErSlhRSIIWmTaIGx2T//y/Da3vtvegyI81oZs75vqeF3VqKViu9O7o49mZLRK23Kf0CiGhxOCZqPxwTtR+OidoPx0Tth2Oi9ot3jHyZnjwZ//r793F/3uT0UwukueHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ni7ecXdTzPSRC6S54di4aMddNwoZx8XCsXGxjrtuHDKOi4Vj43AsE46Ni3TcdROQcVwsHBuHY5lwbFyc466bgozjYuHYOBzLhGPjohx33SRkHBcLx8bhWCYcG4djmXBsXIzjrpuGjONi4dg4HMuEY+NwLBOOjYtwfMq4FzKOi4Vj43AsE46Nw7FMODYu3PE54z7IOC4Wjo3DsUw4Ng7HMuHYOBzLhGPjgh33Me6BjONi4dg4HMuEY+NwLBOOjcOxTDg2LtRxP+NzyDguFo6Nw7FMODYOxzLh2Dgcy4Rj4wIdDzHGcT3h2Liljk8h47hYODYOxzLh2Dgcy4Rj43AsE46Nw7FMODYuzPEw41PIOC4Wjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo0LcjzG+AQyjouFY+NwLBOOjcOxTDg2Dscy4dg4HMuEY+NCHI8zPoaM42Lh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj4wIcTzE+gozjYuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYuCSOu6Pply2Q5oZj46YdBzDGcQ3h2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Lg0jruD6RcukOaGY+NwLBOOjcOxTDg2Dscy4dg4HMuEY+NwLBOOjUvkuHucfuECaW44Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjJlh1OG4mHBuXyvEeMo6LhWPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ni6Z43vIOC4Wjo0bZ9XhuJ1wbByOZcKxcTiWCcfG4VgmHBuHY5lwbByOZcKxcekc30HGcbFwbNwoqw7HDYVj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46NS+j4FjKOi4Vj48ZYdThuKRwbh2OZcGwcjmUaV3n19Zvb/7/+5vLpu6A5qKFwLNOoyg+XX9w6/vTq5fbtlyFzUEvhWKYxlT/98Z934/H1d28ehmYc64RjmYKOq6+ev9tef/v65tHFxQWOZUrpeAcZx8UKcvzh6d7x5BzUUCOsOhw3VeR4PDkHNRSOZQpyzPmxZjiWKcjxp1cvuF4tGI5lmna8+5f7x5LhWCa+n8s4HMuEY+NwLBOOjRtm1cU77nBcMBwbh2OZcGwcjmXCsXE4lgnHxuFYJhwbh2OZcGwcjmXCsXFpHXc4LheOjRtk1eG4sXBsHI5lwrFxOJYJx8bhWCYcG4djmXBsHI5lwrFxOJYJx8YldtzhuFg4Nm6I1aNLHDcSjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYuAFWBy5x3Eg4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYuH5Why5x3Eg4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPj0jvu5iyQlodj43pZHbvEcRvh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Ny+B4HDKOc4Vj4/pYnbrEcRPh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Ny+F4FDKOc4Vj43pYnbvEcQvh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Ny+J4DDKOc4Vj485Z9bnEcQPh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Ny+N4BDKOc4Vj485Y9bvEcf3h2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Ny+R4GDKOc4Vj405ZDbnEcfXh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Nw7FMODbuhNWgSxxXH46Ny+V4EDKOc4Vj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Lhsjocg4zhXODbumNWYSxxXHo6Nw7FMODYOxzLh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Lh8jgcg4zhXODbuiNW4SxzXHY6Nw7FMODYOxzLh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2LiMjvsh4zhXODbukNWUSxxXHY6Nw7FMODYOxzLh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Licjnsh4zhXhyo/Pvvs73FzUNMdsJp2ieOaO1b542bzm1+i5qCGw7FMpyo/Ptts/hA1BzUbjmU6V7mT/Pl/YuagRsOxTKcqf95sfntzfD1ydI1jmXAs05HK/32/2Xy1e/DryICMY5lwLNPx9eqx4+m+OajpcCwT94+Ny+q4DzKOc4Vj4x5ZhbjEccXh2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj43AsE46Ny+u4BzKOc4Vj4x5YhbnEcb3h2Dgcy4Rj43AsE46Nw7FMODYOxzLh2Dgcy4Rj4zI7PoeM41zh2Lg9q1CXOK42HBuHY5lwbByOZcKxcTiWCcfG4VgmHBuHY5lwbByOZcKxcTiWCcfG3bMKdonjasOxcbkdn0HGca5wbByOZcKxcTiWCcfG4VgmHBuHY5lwbByOZcKxcXesIlziuNZwbFx2x6eQcZwrHBuHY5lwbByOZcKxcTiWCcfG4VgmHBuHY5lwbByOZcKxcfkdd30LpPTh2LhbVnEucVxnODYOxzLh2Dgcy4Rj43AsE46Nw7FMYyqvv7l8+u720dvLy8sv3kzPQU2FY5lGVH569XL79svbhz+9DJqD2grHMo2ovP7uzfbq690o/OmH10FzUFvtWEW6jHXcnS2QcjSi8ur5u+31tzvBNwfYl5e3Q/LFxQWOZcKxTCMqPzzdO7768+vHMRnHMuFYpqDx+Lb9OTKOZcKxTEHnx7fhWC4cyzR6vfrF/fXq3RH2p39w30ktHMs0ff94NyS/vbz84+uAOaipcCwT389l3BqOu9MFUo5wbNwNq3iXOK4xHBuHY5lwbByOZcKxcTiWCcfG4VgmHBuHY5lwbNyTJzNcRk9/tMBi6yoejo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYuHUcd4cLLLey2uHYOBzLhGPjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjVnLcHSyw4NpKh2PjcCwTjo3DsUw4Ng7HMuHYOBzLhGPjcCwTjo3DsUw49q1by/EDZBznCse+4VgnHPuGY51w7BuOdcKxbzjWCce+4VgnHPuGY51wbFuHY51wbNuKjveQcZwrHHu1uW977/j+KY4bD8cO3cHdPjLeO354juO2w7F6j3BPOjiufnB8DBrHzYRj8YYU954fnwzOOG4mHIvXb3hX73WuI8g4biYcazYwBh/UDV6vzuD4HjKOc4VjwYaPpQ8adhziEsd1hWO9QhSHON68fz90RwrHlYVjuUIQbwPH4+F7yziuKhy7hmOlcKxT2Di8L8hxuGEcFw3HIoWdFD/WRToOGZbH3N8uFMe5wrFGsYxjHQcdYOO4WDhWKBbxdsZ4HCAZx8XCsULRiuMdhwzJOC4Wjk2Ld3z6tyhwXFE4brsZI/FdcxxPheNi4bjlZpwX39ctcTw0LI+63y0Vx7nCcbvNV7zM8byfH7JbKo5zheNmW8J4qeN+yjguFo6bbYHihcfVQzehcFwsHHu2zHHXf+0ax8XCsWeLHcdPv8VxvnDcXEtOi/d1ONYKx6216PLWvnSODw+vcVwsHLdVEsUJHR9d78JxsXDcVIkYpxyPg39O7hbH+cJxU6VRnPT8OPjn5G5xnC8cO9YldHxwbD0xPY7zhWPH0jrucFw8HDdSqiPq2xI7Dpwex/nCcROlur51H47VwnELJWaMY7lw3ECJFedyvJmaHsfZwnH95WKcYTwe/8E/OM4WjusvNeOcjkd/hBeOs4Vjw7I6HoGM42zh2LB8jsd/piaOs4Xjikt+QH1fTsdj4ThbOK635Ne37ntwtbbjDse5wnGtpb5n/NgajvuPrnGcKxxXWj7GazgeuNyF41zhuNKyKV5pPO6VjONc4divNRz3Q8ZxrnBs1yGrfI57IeM4Vzi2ayXHfbeScZwrHNdWvutb963m+Dwc5wrHlZXxOvV9OBYMx1WVXzGOJcNxVeVXvL7jg9NkHOcKxzW1wmh8wHg1xw+ScZwrHLu1uuMDyDjOFY7dWt/x47dp4jhXOHarhOM9ZBznCsd1tMYFrtsOda14vXrvuFtpNd3CcQ2tcbvpvkKOHxa40mq6heMKWpExjjXDcflWVIxj0XBcvFUZl3bM3pMnHBevGOMijtl9soRjr8o7RnKOcOxVace339zFLpQ8HHtV3HG38vUAk3BcrgL78zmr1R1zxTpHOC5VkeNLHIuG40KVObrEsWg4LlOhk8SaHLMjJQzHRaqCcSHH3cNbwK6ULBw7VZNjrlunDMdOVeWYETlhOHaqLsdATheO163ontvPanXHB1escZwoHK9a2RGoPseUKByvWeHjSBzLhuMVK306iGPZcLxetTEu5vgMMqfJS8OxTzU7RvKycOxTvY6BvDQc+1SxYyQvDMcrVMceOsyqBsfFrx60HY7zV8keWpFjrlinDsfZq+WIEcfC4ThztSjGsXQ4zlwtinsYF3Q8DLmaj722wrFL7ThmD4sPxy414hjJs8JxtirbG5txDOQZ4ThXle2M46xWdzwFea23RSUc56m6MaUlxxQdjrNUHWMca4fjHFWnuJdxUcdBkGt7F+sNxx416bi+j8Nqw7FHTTqu8PSk1nCctkr3u2lWs1zGTn+0wJDXzW9LDgzHKat2rwtgNcdl7PTRjjm2DgzHCauWccOOKSgcp6taxQOMCzsGcsJwnKx6GSs4rva9rSMcOxTEaobL2OnnO673lKWOcGxQGKsZLmOnP1lgzDoAeTQcJ6j2tySQVbzL2OkXOEbyaDheWv27VyireJex058uMG49qn+jC4bjhdXPWMYxDYfjRTWgGMcO4XhRDSgeZlze8SzILXx0rh6O5YtgFesydvpUjqF8Go5n1s6eFMEq1mXs9Ekc87cnesLxvNrZkWJYxbqMnf58gTNXqp23f6VwPKOmxoMoVpEuY6dP5phOwnF8KoxrcAzkNOE4vnYQb9UdN/SBmjccaxfLKs5l7PR9C1y0ely7vg/H2kWzinIZO31yx42d5OQLx8G1uL/Es4pyGTt97wIXr2SD2yV5OA6r0QO4GaxiXMZOn8cx4TiwNhVPMK7EcTLILW6hVOE4pCYRb80cN/pZmyYcCzflrhLHSSGb7p44Hq3t3WLKnZjjbbsHTovD8XCtf7xPuqvFMZe6FofjwRpXPM1Y2XHbmy4+HA/VuOKmHCeHbHeqjGPVAtzV4zj9iGwmGccnyWz9AHfKjrcCR1QR4fiwjc6neIi7ihxzqWtZOD5IBnEY46oc54Sss1kHw7FkYe6qcpwRstBh1lA4VhqF9wW6q8tx1kNrwY18lL1jyc/qQHdGjrfit5TdHesZ3gYzrs3xOte6FD+3fR3rbcnHgt3V5ngVyBvJIzBDx5Lb8aBwd9U5Xu3uk9z293MM42BW6zte+zayzL5g4/hhc0lstsGi3FXoeGXJMgfZFo5VNtZ0ce6qdLz+d3Y97hrt7iTKjh9HYBPFkYwrdVzuWzQbHp1FHbe7QRYUhy6AVSHHJb/Z+nC3aWkHknK8OR6BK36lGYojF8SqmOMq/tZEU6OzjOOW3vQcxYkLY1XOcU2S+55UV7uOjz8u636XsxfnLZRVScd1SH7sdHSubG8bezHX31w+fXfyqKDjPriVvZtlisMWwaqs4642ytuDm5Y9rIvuiiPL/vTq5fbtl8ePxudI8XpO3xvgjhcHI45VccddhZQfOj/oHtp113gxw1+6/u7N9urrN0ePxuc4+8MPin++2Zx+/XAant8R3j/f7/cxz3esxr6+cxnz501Of7PAua+3gvd74vnj/z528nwb/Hwb2cgMV8/fba+/fX306OLiInwRo6s08Pxs2OV5vudPnox//f37uD9vcvqbBda0/qs+v/9vgc+3kY3M8OHpXu/jo/E5qK12jsfauYxpcvqpBdLcIsfj8TmorXAsU9bzY6o7HMs0er36xcP16hdrXa+mFcOxTNP3j3cDcRX3jyl1OJap3e/nosXhWCYcG4djmXBsHI5lwrFxOJYJx8bhWCYcG4djmXBsHI5lwrFxOJYJx8bhWCYcG4djmXBsHI5lwrFxOJYJx8bhWCYcG4djmXBsHI5lwrFxOJYJx8bhWCYcG4djmWY43nexkY7VazmX1Zvv+KGL+bO2EKvXcm6rh+OhWL2Wc1s9HA/F6rWc2+px1Yqo/XBM1H44Jmo/HBO1H46J2m+J40+vXiZ7HdV19afLS9nVO/x9XYJJb7td5/KWOH4r/Gbtft3z1Z9fl34ZedrtBg+/P1Mv6W1327m8BY6v/vJXXccfdrv5T6Lrd/T7rPWS3na7euTNd/zph38pH1dv7z7XJbt6/k535e5SXr0+efMdv30hfX5899vbNfvwVN2x7rbb9sub5finy8svbz7TVR3vVm93LUh2V5Afj4W33e3WS3id6+3lLt236+pPmh9Su8TPj6W33YA87jv1J70r7I46ha9XS2+729Led1J2fPehp7p+2vePtbfdrrSOiaiOcEzUfjgmaj8cE7UfjonaD8dE7YdjovbDMVH74Zio/XBM1H44Nuznz/6+/fjsD6VfBqULx479+NvdP6QTjh377+/+9vv/lH4RlDAcW/bz5qvSL4FShmPLftxwWC0Vjh379fN/P2NAVgrHhn28Qfzr55wgC4Vjw378zS/b/33PkbVQOCZqPxwTtR+OidoPx0Tth2Oi9sMxUfvhmKj9cEzUfjgmaj8cE7Xf/wENKdZLe0nAwAAAAABJRU5ErkJggg==" /><!-- --></p>
<p>Between the blue lines are values the recived <em>some</em> support by the data (this is a BF = 1 SI), which received at least “moderate” support from the data.</p>
<p>From the presepctive of the Savage-Dickey Bayes factor, testing against a point null hypothesis for any value within the support interval will yeild a Bayes factor smaller than 1/<code>BF</code>.</p>
<h1 id="references">References</h1>
<div id="refs" class="references">

<div id="ref-kruschke2014doing">

<p>Kruschke, J. (2014). <em>Doing bayesian data analysis: A tutorial with r, jags, and stan</em>. Academic Press.</p>
</div>

<div id="ref-mcelreath2018statistical">

<p>McElreath, R. (2018). <em>Statistical rethinking: A bayesian course with examples in r and stan</em>. Chapman; Hall/CRC.</p>
</div>

<div id="ref-wagenmakers2018SI">

<p>Wagenmakers, E.-J., Gronau, Q. F., Dablander, F., &amp; Etz, A. (2018). <em>The support interval</em>. <a href="https://doi.org/10.31234/osf.io/zwnxb">https://doi.org/10.31234/osf.io/zwnxb</a></p>
</div>

</div>

</body>
</html>
back to top