title: compare models in pymer4
body: <p>Is there a way in <code>pymer4</code> (or similar python modules for linear mixed-effects regression models) to compare two different models?</p>
  <p>In R, the command <code>anova</code> of the package <code>lme4</code> does this as follows:</p>
  <pre><code>model1 = lmer(DV ~ IV1 + (1|RV), data = data)
  model2 = lmer(DV ~ IV1*IV2 + (1|RV), data = data)
  
  anova(model1, model2)
  </code></pre>
  <p>The <a href="http://eshinjolly.com/pymer4/api.html" rel="nofollow noreferrer">API</a> page of <code>pymer4</code> does not seem to mention this. (I could do this in R, but I would rather do it in python, as I'll have to run this test through different groups through a <code>for</code> loop, for which python provides for flexibility.)</p>