ry 10 iterations</span>
for ( iter in <span class="dv">1</span>:<span class="kw">length</span>(iterations) ) {
    j =<span class="st"> </span>iterations[iter]
    <span class="co"># Get parameters at iteration j</span>
    beta0_j =<span class="st"> </span>output$bias[j]
    beta_j =<span class="st"> </span>output$beta[j,]
    for ( i in <span class="dv">1</span>:<span class="kw">length</span>(yTest) ) {
        pihat_ij =<span class="st"> </span><span class="dv">1</span> /<span class="st"> </span>(<span c