e="3"># log(phi[i ,j, k]) &lt;- alpha[k] + beta[i, k] + gamma[j, k]<br>
</font><font face="Arial" color="#000100" size="1">&nbsp;&nbsp;&nbsp;
</font><font face="Arial" color="#000100" size="3"># }<br><br>
</font><font face="Arial" color="#000100" size="1">&nbsp;&nbsp;&nbsp;
</font><font face="Arial" color="#000100" size="3"># Fit standard Poisson regressions relative to baseline<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;lambda[i, j] ~ dflat()&nbsp;&nbsp;&nbsp;# vague priors 