nt><font face="Arial" color="#000100" size="3">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;for (j in 1 : J) { # loop around sizes<br>
</font><font face="Arial" color="#000100" size="1">&nbsp;&nbsp;&nbsp;
</font><font face="Arial" color="#000100" size="3">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;g[j, k] &lt;- gamma[j, k] - mean(gamma[, k]); # sum to zero constraint<br>
</font><font face="Arial" color="#000100" size="1">&nbsp;&nbsp;&nbsp;
</font><font face="Arial" color="#000100" size="3">&nbsp;&nbsp;&nb