title: NaNs produced when using aodml beta-binomial model in aods3 package
body: <p>When trying to use the <code>aodml</code> funciton to fit a beta-binomial distribution GLM with my own data I get a warning message, stating that there are NaNs produced. When I tried to run the same function with the dataset given in the example of the <code>aods3</code> package I get the same warning messages.</p>
  <p>Code:</p>
  <pre><code>library(aods3)
  data(orob2)
  fm1&lt;-aodml(cbind(m, n-m)~seed, data=orob2, family=&quot;bb&quot;)
  
  Warning messages:
  1: In lbeta(a, b): NaNs produced
  </code></pre>
  <p>I have an inkling that the warning is caused by the <code>cbind(m, n-m)</code> but I don't know why or how? Any ideas, especially strange since this is happening with the example data? Am I missing something here?</p>