\begin{Schunk}
\begin{Sinput}
 anova(myPois1, myNB1, test = "LR")
\end{Sinput}
\begin{Soutput}
Analysis of Deviance Table

Model 1: undercount ~ rural + perAA + equip
Model 2: undercount ~ rural + perAA + equip
  Resid. Df Resid. Dev Df Deviance Pr(>Chi)
1       152      77702                     
2       152        180  0    77522         
\end{Soutput}
\begin{Sinput}
 library(lmtest)
 lrtest(myPois1, myNB1)
\end{Sinput}
\begin{Soutput}
Likelihood ratio test

Model 1: undercount ~ rural + perAA + equip
Model 2: undercount ~ rural + perAA + equip
  #Df LogLik Df Chisq Pr(>Chisq)    
1   7 -39439                        
2   8  -1106  1 76666  < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
\end{Soutput}
\end{Schunk}
