Function simulates time-to-detection occupancy design data under the model of Garrard et al. (2008), also see Bornand et al. (2014)
simOccttd(M = 250, mean.psi = 0.4, mean.lambda = 0.3,
beta1 = 1, alpha1 = -1, Tmax = 10, show.plot = TRUE, verbose = TRUE)
A list with the values of the arguments input and the following additional elements:
Simulated values of covariate A, a vector of length M
Simulated values of covariate B, a vector of length M
Probability of occurrence at each site, a vector of length M
Realized occurrence at each site, a 0/1 vector of length M
Uncensored simulated time-to-detection at each site, a vector of length M
Censored simulated time-to-detection at each site, a vector of length M
Censoring indicator, a 0/1 vector of length M
Total number of sites occupied
Total number of sites where the species was observed
Number of sites
intercept of occupancy probability
intercept of Poisson rate parameter
slope of continuous covariate B on logit(psi)
slope of continuous covariate A on log(lambda)
maximum search time (in arbitrary units, which are same as response), response will be censored at Tmax
choose whether to show plots or not. Set to FALSE when running simulations.
if TRUE, output will be written to the console.
Marc Kéry & Andy Royle
Garrard, G.E., Bekessy, S.A., McCarthy, M.A., & Wintle, B.A. (2008) When have we looked hard enough? A novel method for setting minimum survey effort protocols for flora surveys. Austral Ecology, 33, 986-998.
Bornand, C.N., Kéry, M., Bueche, L., & Fischer, M. (2014) Hide and seek in vegetation: time-to-detection is an efficient design for estimating detectability and occurrence. Methods in Ecology and Evolution, 5, 433-442.
Kéry, M. & Royle, J.A. (2016) Applied Hierarchical Modeling in Ecology AHM1 - 10.12.1.
# Generate data with the default arguments and look at the structure
tmp <- simOccttd()
str(tmp)
Run the code above in your browser using DataLab