Last chance! 50% off unlimited learning
Sale ends in
housemouse
capthist
object. There
are two individual covariates: sex (factor levels `f', `m') and age
class (factor levels `j', `sa', `a'). The sex of two animals is not
available (NA); it is necessary to drop these records for analyses
using `sex' unless missing values are specifically allowed, as in hcov
.
The datasets were originally in the CAPTURE `xy complete' format which
for each detection gives the `column' and `row' numbers of the trap
(e.g. ` 9 5' for a capture in the trap at position (x=9, y=5) on the
grid). Trap identifiers have been recoded as strings with no spaces by
inserting zeros (e.g. `0905' in this example).
Sherman traps are designed to capture one animal at a time, but the data
include 30 double captures and one occasion when there were 4
individuals in a trap at one time. The true detector type therefore
falls between `single' and `multi'. Detector type is set to `multi' in
the distributed data objects.
Otis et al. (1978) report various analyses including a closure test on
the full data, and model selection and density estimation on data from
the mornings only. We include several secr models fitted to the
`morning' data (morning.0
, morning.b
etc.). Of these, a model
including individual heterogeneity in both g0 and sigma has the lowest
AIC.
plot(housemouse, title = paste("Coulombe (1965), Mus musculus,",
"California salt marsh"), border = 5, rad = 0.5,
gridlines = FALSE)
morning <- subset(housemouse, occ = c(1,3,5,7,9))
summary(morning)
## drop 2 unknown-sex mice
known.sex <- subset(housemouse, !is.na(covariates(housemouse)$sex))
## reveal multiple captures
table(trap(housemouse), occasion(housemouse))
AIC(morning.0, morning.b, morning.t, morning.h2, morning.0h2, morning.h2h2)
## assess need to distinguish morning and afternoon samples
housemouse.0 <- secr.fit (housemouse, buffer = 20)
housemouse.ampm <- secr.fit (housemouse, model = g0~tcov, buffer = 20,
timecov = c(0,1,0,1,0,1,0,1,0,1))
AIC(housemouse.0, housemouse.ampm)
Run the code above in your browser using DataLab