gibbscap2 {bayess}

Gibbs sampling for the Arnason-Schwarz capture-recapture model

Description

In the Arnason-Schwarz capture-recapture model (see Chapter 5), individual histories are observed and missing steps can be inferred upon. For the dataset eurodip, the moves between regions can be reconstituted. This is the first instance of a hidden Markov model met in the book.

Usage

gibbscap2(nsimu, z)

Arguments

nsimu

numbed of simulation steps in the Gibbs sampler

z

data, capture history of each individual, with 0 coding non-capture

Value

p

Gibbs sample of capture probabilities across time

phi

Gibbs sample of survival probabilities across time and locations

psi

Gibbs sample of interstata movement probabilities across time and locations

late

Gibbs averages of completed histories

Examples

data(eurodip)
res=gibbscap2(10,eurodip[1:100,])
plot(res$p,type="l",col="steelblue3",xlab="iterations",ylab="p")