This function generates encounter histories from simulated closed population capture-mark-recapture data consisting of multiple non-invasive marks.
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"
)A list containing the following:
A matrix containing the observed encounter histories with rows corresponding to individuals and columns corresponding to sampling occasions.
A matrix containing the true (latent) encounter histories with rows corresponding to individuals and columns corresponding to sampling occasions.
True population size or abundance.
The number of sampling occasions.
Logit- or probit-scale intercept term(s) for capture probability (p). Must be a scaler or vector of length noccas.
Additive logit- or probit-scale behavioral effect term for recapture probability (c).
Logit- or probit-scale individual heterogeneity variance term.
Conditional probability of type 1 encounter, given detection.
Conditional probability of type 2 encounter, given detection.
Conditional probability of simultaneous type 1 and type 2 detection, given both types encountered. Only applies when data.type="sometimes".
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 function for detection probability. Must be "logit" or "probit". Note that multimarkClosed is currently implemented for the logit link only.
Brett T. McClintock
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.
processdata, multimarkClosed
#simulate data for data.type="sometimes" using defaults
data<-simdataClosed(data.type="sometimes")
Run the code above in your browser using DataLab