summary(arvicola)
head(traps(arvicola))
## for speed, pre-compute distance matrix
userd <- networkdistance (traps(arvicola), glymemask, glymemask)
## fit model
glymefit <- secr.fit(arvicola, mask = glymemask, trace = FALSE,
details = list(userdist = userd))
## estimates of 'real' parameters
predict(glymefit)
tmp <- secr.test(glymefit, nsim = 1000)
tmp
plot(tmp)
## More voles were caught only once than is predicted by the model.
## This is probably due to within-population variation in movement or
## capture probability.Run the code above in your browser using DataLab