Hunting: Hunting Data from Pico Basile, Bioko Island, Equatorial...
In mcglm: Multivariate Covariance Generalized Linear Models

Hunting

R Documentation

Hunting Data from Pico Basile, Bioko Island, Equatorial Guinea

Description

This dataset contains a case study analyzed in Bonat et al. (2017) regarding animals hunted in the village of Basile Fang, Bioko Norte Province, Bioko Island, Equatorial Guinea. Monthly counts of blue duikers and other small animals shot or snared were collected from a random sample of 52 commercial hunters between August 2010 and September 2013. For each animal caught, the species, sex, capture method, and altitude were recorded. The dataset contains 1216 observations.

Usage

data(Hunting)

Format

A data.frame with 1216 observations and 11 variables:

ALT: Factor with five levels indicating the altitude where the animal was caught.
SEX: Factor with levels Female and Male.
METHOD: Factor with levels Escopeta and Trampa indicating the method of capture.
OT: Monthly number of other small animals hunted.
BD: Monthly number of blue duikers hunted.
OFFSET: Monthly number of hunter days.
HUNTER: Hunter index.
MONTH: Month index.
MONTHCALENDAR: Month as calendar number (1 = January, ..., 12 = December).
YEAR: Calendar year (2010–2013).
HUNTER.MONTH: Index indicating observations taken for the same hunter and month.

Source

Bonat, W. H., et al. (2017). "Modelling the covariance structure in marginal multivariate count models: Hunting in Bioko Island." Journal of Agricultural, Biological and Environmental Statistics, 22(4):446–464.

Bonat, W. H. (2018). "Multiple Response Variables Regression Models in R: The mcglm Package." Journal of Statistical Software, 84(4):1–30.

Examples

library(mcglm)
library(Matrix)
data(Hunting, package = "mcglm")
formu <- OT ~ METHOD*ALT + SEX + ALT*poly(MONTH, 4)
Z0 <- mc_id(Hunting)
Z1 <- mc_mixed(~0 + HUNTER.MONTH, data = Hunting)
fit <- mcglm(linear_pred = c(formu),
             matrix_pred = list(c(Z0, Z1)),
             link = c("log"),
             variance = c("poisson_tweedie"),
             power_fixed = c(FALSE),
             control_algorithm = list(max_iter = 100),
             offset = list(log(Hunting$OFFSET)),
             data = Hunting)
summary(fit)
anova(fit)

mcglm documentation built on Jan. 9, 2026, 1:07 a.m.