Dear Murray
I hope you are well and thanks for your continuous support in running this fantastic forum!
I recently read your tiger, ovenbird and possum paper and found it an absolutely intriguing reference and model for home-range related camera-trap problems. I'm currently trying to figure out how to apply this to leopards and am re-running your simulations to get a grip on the analysis.
It seems I can get the ovenbird and possum data to work perfectly, however I am struggling to get the tiger data to run! Is the data actually included in the supplementary section of your Ecography paper? The lines of code below don't seem to allow me access to re-run the code:
## Load package ‘secr’ and data (previously prepared multi-session capthist
## object CHc28 and list of masks masks28)
library(secr)
load (’dataobjects.RData’)
## Fully reserve-specific parameters
## -- output is a list of fits, rather than an ‘secr’ object
fit.0c28 <- mapply(secr.fit, CHc28, masks28, SIMPLIFY = FALSE, trace =
FALSE, MoreArgs = list(model = list(D~1, lambda0~1, sigma~1),
detectfn = ’HHN’, verify = FALSE, start = c( -7.7, -4.2, 8.1),
binomN = 1))
## extract estimates from reserve-specific model fits
lapply(fit.0c28, predict)
## Fully reserve-specific, with k parameterization
fit.0c28k <- mapply(secr.fit, CHc28, masks28, SIMPLIFY = FALSE, trace =
FALSE, MoreArgs = list(model = list(D~1, lambda0~1, sigmak~1), detectfn =
’HHN’, verify = FALSE, start = c( -7.7, -4.2, 4.2), binomN = 1))
## extract estimates from reserve-specific model fits
lapply(fit.0c28k, predict)
## Multi-reserve model with constant sigmak
fit.0csl28 <- secr.fit( CHc28, mask = masks28, model = list(D~session,
lambda0~session, sigmak~1), detectfn = ’HHN’, verify = FALSE, binomN = 1)
predict(fit.0csl28)
Best
Alex