nary/Bernoulli) can either be specified in the form <code>prob ~ ..., weights = N</code>, or in the more typical two-column matrix <code>cbind(successes,failures)~...</code> form.</p></li>
<li><p>Behavior of <code>REML=TRUE</code> for Gaussian responses matches <code><a href="https://rdrr.io/pkg/lme4/man/lmer.html" class="external-link">lme4::lmer</a></code>. It may also be useful in some cases with non-Gaussian responses (Millar 2011). Simulations should be done first to verify.</p></li>
<li><p>Because the