library("countreg")
data("CrabSatellites", package = "countreg")
cs <- CrabSatellites[, c("satellites", "width", "color")]
cs$color <- as.numeric(cs$color)
fm <- hurdle(satellites ~ ., data = cs, dist = "negbin", zero = "negbin")
if(require("car")) hurdletest(fm)Wald test for hurdle models
Restrictions:
count_((Intercept) - zero_(Intercept) = 0
count_width - zero_width = 0
count_color - zero_color = 0
Model 1: restricted model
Model 2: satellites ~ .
Res.Df Df Chisq Pr(>Chisq)
1 168
2 165 3 50.362 6.689e-11 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1