closedp.Mtb {Rcapture}R Documentation

Closed Population Capture-Recapture Model Mtb

Description

This function fits model Mtb for closed populations in capture-recapture experiments.

Usage

closedp.Mtb(X, dfreq=FALSE)

## S3 method for class 'closedp.Mtb':
print(x, ...)

Arguments

X The table of the observed capture histories in one of the two accepted formats. In the default format, it has one row per unit captured in the experiment. In this case, the number of columns in the table represents the number of capture occasions in the experiment (noted t). In the alternative format, it contains one row per capture history followed by its frequency. In that case, X has t+1 columns. The first t columns of X, identifying the capture histories, must contain only zeros and ones. The number one indicates a capture.
dfreq This argument specifies the format of the data matrix X. By default, it is set to FALSE, which means that X has one row per unit. If it is set to TRUE, then the matrix X contains frequencies in its last column.
x An object, produced by the closedp.Mtb function, to print.
... Further arguments passed to or from other methods.

Details

The Mtb model is non-linear. It is fitted with the optim function instead of the glm fonction. Therefore, the abundance estimate can be unstable.

For the model to be identifiable, the parameters are constrained in the following way: logit(ci)=logit(pi)+b for i in 2,...,l.

Value

n The number of captured units
results A table containing the estimated population size, the standard error of estimation, the deviance, the number of degrees of freedom and the Akaike criteria.
parMtb Capture-recapture parameters estimates for model Mtb : the abundance N, p1 to pt, the probabilities of first capture for each capture occasion, and c2 to ct, the recapture probabilities for each capture occasion.

Author(s)

Sophie Baillargeon Sophie.Baillargeon@mat.ulaval.ca and
Louis-Paul Rivest Louis-Paul.Rivest@mat.ulaval.ca

References

Baillargeon, S. and Rivest, L.P. (2007). Rcapture: Loglinear models for capture-recapture in R. Journal of Statistical Software, 19(5), http://www.jstatsoft.org/

See Also

closedp, closedp.mX, closedp.h

Examples

data(hare)
closedp.Mtb(hare)

## Example producing an unstable estimate
data(mvole)
period4<-mvole[,16:20]
closedp.Mtb(period4)

[Package Rcapture version 1.1 Index]