simdataClosed: Simulate closed population capture-mark-recapture data...

In multimark: Capture-Mark-Recapture Analysis using Multiple Non-Invasive Marks

Simulate closed population capture-mark-recapture data arising from multiple non-invasive marks

Description

This function generates encounter histories from simulated closed population capture-mark-recapture data consisting of multiple non-invasive marks.

Usage

simdataClosed(
  N = 100,
  noccas = 5,
  pbeta = -0.4,
  tau = 0,
  sigma2_zp = 0,
  delta_1 = 0.4,
  delta_2 = 0.4,
  alpha = 0.5,
  data.type = "never",
  link = "logit"
)

Arguments

N	True population size or abundance.
noccas	The number of sampling occasions.
pbeta	Logit- or probit-scale intercept term(s) for capture probability (p). Must be a scaler or vector of length noccas.
tau	Additive logit- or probit-scale behavioral effect term for recapture probability (c).
sigma2_zp	Logit- or probit-scale individual heterogeneity variance term.
delta_1	Conditional probability of type 1 encounter, given detection.
delta_2	Conditional probability of type 2 encounter, given detection.
alpha	Conditional probability of simultaneous type 1 and type 2 detection, given both types encountered. Only applies when data.type="sometimes".
data.type	Specifies the encounter history data type. All data types include non-detections (type 0 encounter), type 1 encounter (e.g., left-side), and type 2 encounters (e.g., right-side). When both type 1 and type 2 encounters occur for the same individual within a sampling occasion, these can either be "non-simultaneous" (type 3 encounter) or "simultaneous" (type 4 encounter). Three data types are currently permitted: data.type="never" indicates both type 1 and type 2 encounters are never observed for the same individual within a sampling occasion, and observed encounter histories therefore include only type 1 or type 2 encounters (e.g., only left- and right-sided photographs were collected). Observed encounter histories can consist of non-detections (0), type 1 encounters (1), and type 2 encounters (2). See bobcat. Latent encounter histories consist of non-detections (0), type 1 encounters (1), type 2 encounters (2), and type 3 encounters (3). data.type="sometimes" indicates both type 1 and type 2 encounters are sometimes observed (e.g., both-sided photographs are sometimes obtained, but not necessarily for all individuals). Observed encounter histories can consist of non-detections (0), type 1 encounters (1), type 2 encounters (2), type 3 encounters (3), and type 4 encounters (4). Type 3 encounters can only be observed when an individual has at least one type 4 encounter. Latent encounter histories consist of non-detections (0), type 1 encounters (1), type 2 encounters (2), type 3 encounters (3), and type 4 encounters (4). data.type="always" indicates both type 1 and type 2 encounters are always observed, but some encounter histories may still include only type 1 or type 2 encounters. Observed encounter histories can consist of non-detections (0), type 1 encounters (1), type 2 encounters (2), and type 4 encounters (4). Latent encounter histories consist of non-detections (0), type 1 encounters (1), type 2 encounters (2), and type 4 encounters (4).
link	Link function for detection probability. Must be "logit" or "probit". Note that multimarkClosed is currently implemented for the logit link only.

Value

A list containing the following:

Enc.Mat	A matrix containing the observed encounter histories with rows corresponding to individuals and columns corresponding to sampling occasions.
trueEnc.Mat	A matrix containing the true (latent) encounter histories with rows corresponding to individuals and columns corresponding to sampling occasions.

Author(s)

Brett T. McClintock

References

Bonner, S. J., and Holmberg J. 2013. Mark-recapture with multiple, non-invasive marks. Biometrics 69: 766-775.

McClintock, B. T., Conn, P. B., Alonso, R. S., and Crooks, K. R. 2013. Integrated modeling of bilateral photo-identification data in mark-recapture analyses. Ecology 94: 1464-1471.

See Also

processdata, multimarkClosed

Examples

#simulate data for data.type="sometimes" using defaults
data<-simdataClosed(data.type="sometimes")

multimark documentation built on March 31, 2023, 9:33 p.m.