# \donttest{
data(salamanders, package = "HZIP")
## Fit zero-inflated Poisson with random effects
fit.salamander <- hzip(y ~ mined | mined + spp + mined * spp,
data = salamanders)
summary(fit.salamander)
# }
Run the code above in your browser using DataLab