secr (version 1.4.0)

stoatDNA: Stoat DNA data

Description

Data of A. E. Byrom from a study of stoats (Mustela erminea) in New Zealand. Individuals were identified from DNA in hair samples.

Usage

data(stoatDNA)

source

Efford, M. G., Borchers D. L. and Byrom, A. E. (2009) Density estimation by spatially explicit capture-recapture: likelihood-based methods. In: D. L. Thomson, E. G. Cooch and M. J. Conroy (eds) Modeling Demographic Processes in Marked Populations. Springer, New York. Pp. 255--269.

Details

The data are from a pilot study of stoats in red beech (Nothofagus fusca) forest in the Matakitaki Valley, South Island, New Zealand. Sticky hair-sampling tubes (n = 94) were placed on a 3-km x 3-km grid with 500-m spacing between lines and 250-m spacing along lines. Tubes were baited with rabbit meat and checked daily for 7 days, starting on 15 December 2001. Stoat hair samples were identified to individual using DNA microsatellites amplified by PCR from follicular tissue (Gleeson et al. 2010). Six loci were amplified and the mean number of alleles was 7.3 per locus. Not all loci could be amplified in 27% of samples. A total of 40 hair samples were collected (Gleeson et al. 2010), but only 30 appear in this dataset; the rest presumably did not yield sufficient DNA for genotyping. The data are provided as a single-session capthist object 'stoatCH'. Hair tubes are treated as 'proximity' detectors which allow an individual to be detected at multiple detectors on one occasion (day), although there are no multiple detections in this dataset. Three pre-fitted models are included: stoat.model.HN, stoat.model.HZ, and stoat.model.EX (with halfnormal, hazard-rate and negative exponential detection functions, respectively).

References

Gleeson, D. M., Byrom, A. E. and Howitt, R. L. J. (2010) Non-invasive methods for genotyping of stoats (Mustela erminea) in New Zealand: potential for field applications. New Zealand Journal of Ecology in press. Available on-line at http://www.newzealandecology.org.

See Also

capthist, detection functions, secr.fit

Examples

data(stoatDNA)
summary(stoatCH)

stoat.model.HN <- secr.fit(stoatCH, buffer = 1000, detectfn = 0)
stoat.model.HZ <- secr.fit(stoatCH, buffer = 1000, detectfn = 1)
stoat.model.EX <- secr.fit(stoatCH, buffer = 1000, detectfn = 2)
confint(stoat.model.HN, 'D')
## Profile likelihood interval(s)...
##         lcl        ucl
## D 0.01381419 0.04563511

## plot fitted detection functions
xv <- seq(0,800,10)
plot(stoat.model.EX, xval = xv, ylim = c(0,0.12), limits = FALSE,
    lty = 2)
plot(stoat.model.HN, xval = xv, limits = FALSE, lty = 1, add = TRUE)
plot(stoat.model.HZ, xval = xv, limits = FALSE, lty = 3, add = TRUE)

## review density estimates
collate(stoat.model.HZ, stoat.model.HN, stoat.model.EX,
    realnames='D', perm=c(2,3,4,1))
model.average(stoat.model.HZ, stoat.model.HN, stoat.model.EX,
    realnames='D')