diff --git a/Stats/BBP-click-whistles_3models.R b/Stats/BBP-click-whistles_3models.R
index 09843c4d63bdd424c3e912f481f5dcca5482b9ab..d00efb1da45548877a58c5bf3b93034eaabc0ecd 100644
--- a/Stats/BBP-click-whistles_3models.R
+++ b/Stats/BBP-click-whistles_3models.R
@@ -196,10 +196,9 @@ mod.bbp$deviance/mod.bbp$df.residual
 
 ### Model for clicks
 # Using NB model:
-mod.cli <- gilm.nb(number_of_clicks ~ acoustic + fishing_net + acoustic:fishing_net + offset(log(number)), 
+mod.cli <- glm.nb(number_of_clicks ~ acoustic + fishing_net + acoustic:fishing_net + offset(log(number)), 
                data=acoustic.dta)
 car::Anova(mod.cli, type=3)
-shapiro.test(residuals(mod.cli)) # H0 : normality -> cannot be rejected if p > 0.05
 dwtest(mod.cli) # H0 -> independent if p>0.05 (autocorrelation if p<0.05)
 bptest(mod.cli) # H0 -> homoscedasticity if p<0.05
 mod.cli$deviance/mod.cli$df.residual