This function produces descriptive statistics for capture-recapture data.
descriptive(X, dfreq=FALSE, dtype=c("hist","nbcap"), t=NULL)# S3 method for descriptive
print(x, …)
# S3 method for descriptive
plot(x, main="Exploratory Heterogeneity Graph", …)
The matrix of the observed capture histories (see Rcapture-package for a description of the accepted formats).
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.
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).
Requested only if dtype="nbcap". A numeric specifying the total number
of capture occasions in the experiment. The value t=Inf
is accepted. It indicates that captures occur in continuous time.
(see Rcapture-package).
An object, produced by the descriptive function, to print or to plot.
A main title for the plot
Further arguments to be passed to methods (see print.default and plot.default).
The total number of captured units.
A table of basic descriptive statistics. For \(i=1,\ldots,t\) :
fi: the number of units captured \(i\) times
ui: the number of units captured for the first time on occasion \(i\)
vi: the number of units captured for the last time on occasion \(i\)
ni: the number of units captured on occasion \(i\)
If dtype="nbcap", only the the number of units captured \(i\) times, fi, can be computed.
Only if dtype="hist": A \(t\) by \(t\) matrix containing recapture statistics for units
released on occasion \(i\), for \(i=1,\ldots,t\) :
ni: the number of units released on occasion \(i\);
c2: the number of units released on occasion \(i\) and recaptured for the first time on occasion 2;
…
ct: the number of units released on occasion \(i\) and recaptured for the first time on occasion \(t\);
not recapt: the number of units released on occasion \(i\) and never recaptured again.
The m.array matrix is used in Mark (see White and Burnham 1999) and MSurge (see Choquet, Reboulet, Pradel, Gimenez,
and Lebreton 2004) to fit open population models to the data.
The function call (object of class "call").
The plot.descriptive function produces an exploratory heterogeneity graph.
In the absence of heterogeneity, the relation(s) presented in the graph should be almost linear.
Convex functions indicate heterogeneity (Baillargeon and Rivest (2007)).
Baillargeon, S. and Rivest, L.P. (2007) Rcapture: Loglinear models for capture-recapture in R. Journal of Statistical Software, 19(5), 10.18637/jss.v019.i05.
Choquet, R., Reboulet, A.M., Pradel, R., Gimenez, O. and Lebreton, J.D. (2004). M-Surge: New Software Specifically Designed for Multistate Capture-Recapture Models. Animal Biodiversity and Conservation, 27.1, 207--215.
Hoaglin, D. C. (1980) A Poissonness Plot, The American Statistician, 34, 146--149
Lindsay, B. G. (1986) Exponential family mixture models (with least-squares estimators). Annals of Statistics, 14, 124--137.
Rivest, L.P. (2008) Why a time effect often has a limited impact on capture-recapture estimates in closed populations. Canadian Journal of Statistics, 36(1), 75--84.
White, G. and Burnham, K.P. (1999) Program Mark: Survival Estimation from Populations of Marked Animals. Bird Study, 46 (Supplement), 120--138. The software can be downloaded here: http://www.phidot.org/software/mark/.
# NOT RUN {
# hare data set
hare.desc <- descriptive(hare)
hare.desc
plot(hare.desc)
# cbird data set
cbird.desc1 <- descriptive(cbird, dfreq = TRUE, dtype = "nbcap", t = 11)
plot(cbird.desc1)
# To illustrate the option t = Inf.
cbird.desc2 <- descriptive(cbird, dfreq = TRUE, dtype = "nbcap", t = Inf)
plot(cbird.desc2)
# The y coordinate has changed.
# }
Run the code above in your browser using DataLab