From 9211021151dc12c10267767f7b52e2fb10e9a5e7 Mon Sep 17 00:00:00 2001
From: Loic-Lenof <loic.lenof@gmail.com>
Date: Thu, 8 Sep 2022 10:52:02 +0200
Subject: [PATCH] Correction of errors
---
Stats/BBP-click-whistles_3models.R | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Stats/BBP-click-whistles_3models.R b/Stats/BBP-click-whistles_3models.R
index 8a05c3c..09843c4 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)
--
GitLab