closedp.Mtb {Rcapture} | R Documentation |
This function fits model Mtb for closed populations in capture-recapture experiments.
closedp.Mtb(X, dfreq=FALSE) ## S3 method for class 'closedp.Mtb': print(x, ...)
X |
The matrix of the observed capture histories (see Rcapture-package for a description of the accepted formats). |
dfreq |
A logical. By default FALSE, which means that X has one row per unit. If TRUE, it indicates that the matrix X contains frequencies in its last column. |
x |
An object, produced by the closedp.Mtb function, to print. |
... |
Further arguments to be passed to print.default . |
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.
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. |
Sophie Baillargeon Sophie.Baillargeon@mat.ulaval.ca and
Louis-Paul Rivest Louis-Paul.Rivest@mat.ulaval.ca
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/v19/i05.
data(hare) closedp.Mtb(hare) ## Example producing an unstable estimate data(mvole) period4<-mvole[,16:20] closedp.Mtb(period4)