closedp {Rcapture} | R Documentation |
The functions closedp.t
and closedp.0
fit various loglinear models for closed populations in capture-recapture experiments. For back compatibility, closedp.t
is also named closedp
.
closedp.t
fits more models than closedp.0
but for data set with more than 20 capture occasions, the function migth fail. However, closedp.0
works with fairly large data sets (see Details).
closedp(X, dfreq=FALSE, neg=TRUE, trace=FALSE) closedp.t(X, dfreq=FALSE, neg=TRUE, trace=FALSE) closedp.0(X, dfreq=FALSE, dtype=c("hist","nbcap"), t, t0=t, neg=TRUE, trace=FALSE) ## S3 method for class 'closedp': print(x, ...) ## S3 method for class 'closedp': boxplot(x, main="Boxplots of Pearson Residuals", ...) ## S3 method for class 'closedp': plot(x, main="Residual plots for some heterogeneity models", ...)
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. |
dtype |
A characters string, either "hist" or "nbcap", to specify the type of data. "hist", the default, means that X contains complete observed capture histories. "nbcap" means that X contains numbers of captures (see Rcapture-package for details on data formats). |
t |
Requested only if dtype="nbcap" . A numeric specifying the total number of capture occasions in the experiment. |
t0 |
A numeric. Models are fitted considering only the frequencies of units captured 1 to t0 times. By default t0=t . |
neg |
If this option is set to TRUE, negative eta parameters in Chao's models are set to zero. |
trace |
A logical, by default FALSE. If set to TRUE, a note is printed for each model while the function runs. It is useful to identify which model is associated to a warning. |
x |
An object, produced by a closedp function, to print or to plot. |
main |
A main title for the plot |
... |
Further arguments to be passed to methods (see print.default , boxplot.default and plot.default ). |
closedp.t
fits models M0, Mt, Mh Chao, Mh Poisson2, Mh Darroch, Mh Gamma3.5, Mth Chao, Mth Poisson2, Mth Darroch, Mth Gamma3.5, Mb and Mbh. closedp.0
fits only models M0, Mh Chao, Mh Poisson2, Mh Darroch and Mh Gamma3.5. However, closedp.0
can be used with larger data sets than closedp.t
. This is explained by the fact that closedp.t
fits models using the frequencies of the observable capture histories (vector of size 2^t-1), whereas closedp.0
uses the numbers of units captured i times, for i=1,...,t (vector of size t).
closedp.0
has an additional argument t0
which gives to the numbers of units caught more than t0
times their own parameters in the loglinear model. For example, the model for Mh Gamma3.5 has 3+t-t0 parameters. This means that closedp.0
fits models considering only the frequencies of units captured 1 to t0
times.
Chao's models estimate a lower bound for the abundance, both with a time effect (Mth Chao) and without (Mh Chao). The estimate obtained under Mh Chao is Chao's (1987) moment estimator. Rivest and Baillargeon (2007) exhibit a loglinear model underlying this estimator and provide a generalization to Mth. For these two models, a small deviance means that there is an heterogeneity in capture probabilities; it does not mean that the lower bound estimates are unbiased.
Other models for heterogeneity are defined as follows :
Model | Colomn for heterogeneity in the design matrix |
Poisson2 | 2^k-1 |
Darroch | k^2/2 |
Gamma3.5 | -log(3.5 + k) + log(3.5) |
where k is the number of captures. Poisson and Gamma models with parameter values different than 2 and 3.5 respectively can be fitted with the closedpCI.t
and closedpCI.0
functions.
Darroch's models for Mh and Mth are considered by Darroch et al. (1993) and Agresti (1994). Poisson and Gamma models are discussed in Rivest and Baillargeon (2007). Poisson models typically yield smaller corrections for heterogeneity than Darroch's model since the capture probabilities are bounded from below under these models. On the other hand, Gamma models can lead to very large estimators of abundance. We suggest considering this estimator in experiments where very small capture probabilities are likely.
When the variance of an abundance estimate is large, it is useful to use the closedpCI.t
or closedpCI.0
function to construct a profile condifence interval for this abundance.
The boxplot.closedp
function produces boxplots of the Pearson residuals of the fitted loglinear models that converged.
The plot.closedp
function produces scatterplots of the Pearson residuals in terms of fi (number of units captured i times) for the heterogeneous models Mh Poisson2, Mh Darroch and Mh Gamma3.5 if they converged.
To calculate bias corrected abundance estimates, use the closedp.bc
function.
n |
The number of captured units. | ||||||||||
t |
The total number of capture occasions in the experiment. | ||||||||||
results |
A table containing, for every fitted model, the estimated population size, the standard error of estimation, the deviance, the number of degrees of freedom and the Akaike criteria. | ||||||||||
converge |
A logical vector indicating whether or not the fitted models converged. | ||||||||||
glm |
A list of the 'glm' objects obtained from fitting models. | ||||||||||
parameters |
Capture-recapture parameters estimates. It contains N, the estimated population size, and p or p1 to pt defined as follows for the different models :
| ||||||||||
neg.eta |
The position of the eta parameters set to zero in the loglinear parameter vector of model MhC and MthC. | ||||||||||
X |
A copy of the data given as input in the function call. | ||||||||||
dfreq |
A copy of the dfreq argument given in the function call. | ||||||||||
t0 |
For closedp.0 only, a copy of the t0 argument given in the function call. |
This function uses the glm
function of the stats
package.
Sophie Baillargeon Sophie.Baillargeon@mat.ulaval.ca and
Louis-Paul Rivest Louis-Paul.Rivest@mat.ulaval.ca
Agresti, A. (1994) Simple capture-recapture models permitting unequal catchability and variable sampling effort. Biometrics, 50, 494–500.
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.
Chao, A. (1987) Estimating the population size for capture-recapture data with unequal catchabililty. Biometrics, 45, 427–438.
Darroch, S.E., Fienberg, G., Glonek, B. and Junker, B. (1993) A three sample multiple capture-recapture approach to the census population estimation with heterogeneous catchability. Journal of the American Statistical Association, 88, 1137–1148.
Rivest, L.P. and Levesque, T. (2001) Improved log-linear model estimators of abundance in capture-recapture experiments. Canadian Journal of Statistics, 29, 555–572.
Rivest, L.P. and Baillargeon, S. (2007) Applications and extensions of Chao's moment estimator for the size of a closed population. Biometrics, 63(4), 999–1006.
Seber, G.A.F. (1982) The Estimation of Animal Abundance and Related Parameters, 2nd edition, New York: Macmillan.
closedpCI.t
, closedpCI.0
, closedp.bc
, closedp.Mtb
, uifit
.
data(hare) hare.closedp<-closedp.t(hare) hare.closedp boxplot(hare.closedp) data(mvole) period3<-mvole[,11:15] closedp.t(period3) data(BBS2001) BBS.closedp<-closedp.0(BBS2001,dfreq=TRUE,dtype="nbcap",t=50,t0=20) BBS.closedp plot(BBS.closedp) # Seber (1982) p.107 # When there is 2 capture occasions, the heterogeneity models cannot be fitted X <- matrix(c(1,1,167,1,0,781,0,1,254),byrow=TRUE,ncol=3) closedp.t(X,dfreq=TRUE)