diff --git a/Stats/BBP-click-whistles_3models.R b/Stats/BBP-click-whistles_3models.R
index 8a05c3c6465cc47de8ecd404573a760bd692921d..09843c4d63bdd424c3e912f481f5dcca5482b9ab 100644
--- a/Stats/BBP-click-whistles_3models.R
+++ b/Stats/BBP-click-whistles_3models.R
@@ -168,7 +168,7 @@ par(mfrow=c(1,1))
 zero.whi <- zeroinfl(total_whistles_duration ~ 
                       acoustic + fishing_net + behavior + offset(log(number)), 
                     data=acoustic.dta, dist='negbin')
-nb.whi <- glim.nb(total_whistles_duration ~ 
+nb.whi <- glm.nb(total_whistles_duration ~ 
                      acoustic + fishing_net + behavior + offset(log(number)), 
                    data=acoustic.dta)
 # comparison ZINB VS NB model
@@ -184,7 +184,7 @@ mod.whi$df.null/mod.whi$df.residual
 # No normality of residuals for GLM
 # overdispersion with GLIM quasipoisson
 #try with glim NB:
-mod.bbp <- glim.nb(number_of_bbp ~ acoustic + fishing_net + behavior 
+mod.bbp <- glm.nb(number_of_bbp ~ acoustic + fishing_net + behavior 
                   + offset(log(number)),
                   data=acoustic.dta)
 car::Anova(mod.bbp, type=3)