Rcapture (version 1.4-3)

descriptive: Descriptive Statistics

Description

This function produces descriptive statistics for capture-recapture data.

Usage

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", …)

Arguments

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. The value t=Inf is accepted. It indicates that captures occur in continuous time. (see Rcapture-package).

x

An object, produced by the descriptive function, to print or to plot.

main

A main title for the plot

Further arguments to be passed to methods (see print.default and plot.default).

Value

n

The total number of captured units.

base.freq

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.

m.array

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.

call

The function call (object of class "call").

Details

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)).

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/v19/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/.

Examples

Run this code
# 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 DataCamp Workspace