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 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. |
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/
closedp
, closedp.mX
, closedp.h
data(hare) closedp.Mtb(hare) ## Example producing an unstable estimate data(mvole) period4<-mvole[,16:20] closedp.Mtb(period4)