length 2 giving the parameters of the beta prior for <code>P[k,r]</code>. That is, the prior is <code>P[k,r]</code> <span class="math inline">\(\sim\)</span> Beta(<code>a[1]</code>,<code>a[2]</code>). Because the Beta distribution has two parameters, there will be 2 parameters for each locus and each population. So the output of your function should be a 2 x ‘K’ x <code>R</code> array (where ‘K=2’ because we have a mixture of 2 populations).</p></li>
</ol>
</div>
<div id="answer-templates" class="se