Metropolis-Hastings en Scala
Tengo la sensación de que un lenguaje funcional (como Scala) está particularmente bien adaptado al tipo de operaciones que exige MCMC.
Juzguen Vds.
Primero, genero datos en R:
datos <- rnorm(500, 0.7, 1)
writeLines(as.character(datos), "/tmp/datos.txt")
<div class="highlight"><pre tabindex="0" style="color:#f8f8f2;background-color:#272822;-moz-tab-size:4;-o-tab-size:4;tab-size:4;"><code class="language-R" data-lang="R"><span style="display:flex;"><span>Son de una normal con media <span style="color:#ae81ff">0.7</span>. En el modelo que vamos a crear, suponemos <span style="color:#a6e22e">conocida </span>(e igual a <span style="color:#ae81ff">1</span>) la varianza de la normal y trataremos de estimar la media suponiéndole una distribución a priori normal estándar. Y con Scala, así<span style="color:#f92672">:</span>
</span></span><span style="display:flex;"><span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"><</span>div class<span style="color:#f92672">=</span><span style="color:#e6db74">"highlight"</span><span style="color:#f92672">><</span>pre tabindex<span style="color:#f92672">=</span><span style="color:#e6db74">"0"</span> style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f8f8f2;background-color:#272822;-moz-tab-size:4;-o-tab-size:4;tab-size:4;"</span><span style="color:#f92672">><</span>code class<span style="color:#f92672">=</span><span style="color:#e6db74">"language-Scala"</span> data<span style="color:#f92672">-</span>lang<span style="color:#f92672">=</span><span style="color:#e6db74">"Scala"</span><span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>import<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> scala.io.Source
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>import<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> java.io.PrintWriter
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> read data <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>val<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> source <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>Source<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>fromFile<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#e6db74"</span><span style="color:#f92672">>&</span><span style="color:#75715e">#34;/tmp/datos.txt&#34;</span><span style="color:#f92672">,</span> <span style="color:#e6db74">&#34;UTF-8&#34;</span><span style="color:#f92672">)</span></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>val<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> data <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> source<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>mkString<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>split<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#e6db74"</span><span style="color:#f92672">>&</span><span style="color:#75715e">#34;\\s+&#34;</span><span style="color:#f92672">).</span>map<span style="color:#f92672">(</span><span style="color:#66d9ef">_</span><span style="color:#f92672">.</span>toDouble<span style="color:#f92672">).</span>toList</span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> we need random numbers <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>var<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> r <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> scala<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>util<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>Random<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> <span style="color:#a6e22e">almost </span>(standard) dnorm <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>def<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> dnorm<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>x<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> math<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>exp<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672">-</</span>span<span style="color:#f92672">></span>x<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>*</</span>span<span style="color:#f92672">></span>x<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>/</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">2</span><span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> chain transitions <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>def<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> rg<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>y<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> y <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>+</</span>span<span style="color:#f92672">></span> r<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>nextGaussian
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> <span style="color:#a6e22e">prior </span>(standard normal) times likelihood <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>def<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> dpost<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>x<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> dnorm<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>x<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>*</</span>span<span style="color:#f92672">></span> data<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>map<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>y <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=&</span>gt;<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> dnorm<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>x <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>-</</span>span<span style="color:#f92672">></span> y<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)).<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>reduce<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>_<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>*</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>_<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> metropolis<span style="color:#f92672">-</span>hastings acceptance ratio <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>def<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> acceptanceRatio<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>from<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>,<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> to<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> dpost<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>to<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>/</</span>span<span style="color:#f92672">></span> dpost<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>from<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>def<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> posterior<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>a<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>:</</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Stream<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>[<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>Double<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>]<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>{<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>lazy<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>val<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> elegidos <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>Stream<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>continually<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>a<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> map rg filter <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>x <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=&</span>gt;<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> r<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>nextDouble <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>&</span>lt;<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> acceptanceRatio<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>a<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>,<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> x<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>))<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span> a <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span><span style="color:#75715e">#</span><span style="color:#66d9ef">:</span><span style="color:#66d9ef">:</span> <span style="color:#66d9ef">posterior</span><span style="color:#f92672">(</span><span style="color:#66d9ef">elegidos.head</span><span style="color:#f92672">)</span></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>}<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> Get values and drop burning obs <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>val<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> res <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> posterior<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">0.4</span><span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>).<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>take<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">10000</span><span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>).<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>toList<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>drop<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">2000</span><span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> Thin <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>val<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> res_thin <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> res zip <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>Stream<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> from <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">1</span><span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> filter <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>_<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>_2 <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>%<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">5</span><span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>==</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">0</span><span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> map <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>_<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>_1<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#75715e"</span><span style="color:#f92672">>/*</span> Export <span style="color:#f92672">*/</</span>span<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>val<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> out <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#66d9ef"</span><span style="color:#f92672">></span>new<span style="color:#f92672"></</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>PrintWriter<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>(<span style="color:#f92672"></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#e6db74"</span><span style="color:#f92672">>&</span><span style="color:#75715e">#34;/tmp/posterior.txt&#34;</span><span style="color:#f92672">)</span></span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>res map out<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>println
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>out<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>.<<span style="color:#f92672">/</span>span<span style="color:#f92672">></span>close<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">></span>()<span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">></</span>code<span style="color:#f92672">></</span>pre<span style="color:#f92672">></</span>div<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span>
</span></span><span style="display:flex;"><span>Leo los datos de la posteriori en R y los represento así<span style="color:#f92672">:</span>
</span></span><span style="display:flex;"><span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"><</span>div class<span style="color:#f92672">=</span><span style="color:#e6db74">"highlight"</span><span style="color:#f92672">><</span>pre tabindex<span style="color:#f92672">=</span><span style="color:#e6db74">"0"</span> style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f8f8f2;background-color:#272822;-moz-tab-size:4;-o-tab-size:4;tab-size:4;"</span><span style="color:#f92672">><</span>code class<span style="color:#f92672">=</span><span style="color:#e6db74">"language-R"</span> data<span style="color:#f92672">-</span>lang<span style="color:#f92672">=</span><span style="color:#e6db74">"R"</span><span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">></span>posterior <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>&</span>lt;<span style="color:#f92672">-</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>as.double<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>(<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>readLines<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>(<span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#e6db74"</span><span style="color:#f92672">>&</span><span style="color:#75715e">#34;/tmp/posterior.txt&#34;</span>))</span>
</span></span><span style="display:flex;"><span><span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"display:flex;"</span><span style="color:#f92672">><</span>span<span style="color:#f92672">><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#a6e22e"</span><span style="color:#f92672">></span>hist<span style="color:#f92672"></</span>span<span style="color:#f92672">></span>(posterior, breaks <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#f92672"</span><span style="color:#f92672">>=</</span>span<span style="color:#f92672">></span> <span style="color:#f92672"><</span>span style<span style="color:#f92672">=</span><span style="color:#e6db74">"color:#ae81ff"</span><span style="color:#f92672">></span><span style="color:#ae81ff">50</span><span style="color:#f92672"></</span>span<span style="color:#f92672">></span>)<span style="color:#f92672"></</span>span<span style="color:#f92672">></</span>span<span style="color:#f92672">></</span>code<span style="color:#f92672">></</span>pre<span style="color:#f92672">></</span>div<span style="color:#f92672">></span>
</span></span><span style="display:flex;"><span>
</span></span><span style="display:flex;"><span><span style="color:#f92672">!</span>[posteriori_scala<span style="color:#a6e22e">]</span>(<span style="color:#f92672">/</span>wp<span style="color:#f92672">-</span>uploads<span style="color:#f92672">/</span><span style="color:#ae81ff">2016</span><span style="color:#f92672">/</span><span style="color:#ae81ff">06</span><span style="color:#f92672">/</span>posteriori_scala.png<span style="color:#75715e">#center)</span></span></span></code></pre></div>
