title: How to report Likelihood Ratio Test results
body: <p>I am using a Likelihood Ratio Test (in R) to look for main effects in my model with three fixed factors (site, year, habitat) like this:</p>
  
  <pre><code>model1&lt;-glm(tot.mass~hab, data=biom, family = Gamma(link = "log"))
  anova(model1, test="Chisq")
  model2&lt;-glm(tot.mass~year, data=biom, family = Gamma(link = "log"))
  anova(model2, test="Chisq")
  model3&lt;-glm(tot.mass~site, data=biom, family = Gamma(link = "log"))
  anova(model3, test="Chisq")
  </code></pre>
  
  <p>I was told that by using "anova(model3, test="Chisq")" I will find whether (in this example) adding in the term "site" significantly improves the model. From what I understand, since my p value is >0.5, site doesn't have a significant effect on my model (e.g. isn't a main effect for this response variable). My output looks like this:</p>
  
  <pre><code>     Df Deviance Resid. Df Resid. Dev Pr(&gt;Chi)
  NULL                   687     1681.2         
  site  2   16.543       685     1664.7   0.2884
  </code></pre>
  
  <p>My questions are 1) am I interpreting this correctly and 2) how would I report these results? Would I report the Deviance and pvalue ("Pr(>Chi)")?</p>
  
  <p>Thank you!</p>