Skip to content
Snippets Groups Projects
Commit 92110211 authored by Loic-Lenof's avatar Loic-Lenof
Browse files

Correction of errors

parent 4d1516d8
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Stats/BBP-click-whistles_3models.R
+ 2
− 2
View file @ 92110211
...@@ -168,7 +168,7 @@ par(mfrow=c(1,1)) ...@@ -168,7 +168,7 @@ par(mfrow=c(1,1))
zero.whi <- zeroinfl(total_whistles_duration ~ zero.whi <- zeroinfl(total_whistles_duration ~
acoustic + fishing_net + behavior + offset(log(number)), acoustic + fishing_net + behavior + offset(log(number)),
data=acoustic.dta, dist='negbin') 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)), acoustic + fishing_net + behavior + offset(log(number)),
data=acoustic.dta) data=acoustic.dta)
# comparison ZINB VS NB model # comparison ZINB VS NB model
...@@ -184,7 +184,7 @@ mod.whi$df.null/mod.whi$df.residual ...@@ -184,7 +184,7 @@ mod.whi$df.null/mod.whi$df.residual
# No normality of residuals for GLM # No normality of residuals for GLM
# overdispersion with GLIM quasipoisson # overdispersion with GLIM quasipoisson
#try with glim NB: #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)), + offset(log(number)),
data=acoustic.dta) data=acoustic.dta)
car::Anova(mod.bbp, type=3) car::Anova(mod.bbp, type=3)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment