VGAM (version 0.9-3)

aux.posbernoulli.t: Auxiliary Function for the Positive Bernoulli Family Function with Time Effects

Description

Returns behavioural effects indicator variables from a capture history matrix.

Usage

aux.posbernoulli.t(y, check.y = FALSE, rename = TRUE, name = "bei")

Arguments

y

Capture history matrix. Rows are animals, columns are sampling occasions, and values should be 0s and 1s only.

check.y

Logical, if TRUE then some basic checking is performed.

rename, name

If rename = TRUE then the behavioural effects indicator are named using the value of name as the prefix. If FALSE then use the same column names as y.

Value

• A list with the following components.

  [object Object],[object Object],[object Object],[object Object],[object Object]

Details

This function can help fit certain capture--recapture models (commonly known as $M_{tb}$ or $M_{tbh}$ (no prefix $h$ means it is an intercept-only model) in the literature). See posbernoulli.t for details.

See Also

posbernoulli.t, deermice.

Examples

# Fit a M_tbh model to the deermice data:
(pdata <- aux.posbernoulli.t(with(deermice, cbind(y1, y2, y3, y4, y5, y6))))

deermice <- data.frame(deermice,
                    bei = 0,  # Add this
                    pdata$cap.hist1)  # Incorporate these
head(deermice)  # Augmented with behavioural effect indicator variables
tail(deermice)