## Simple example: generate and summarise trapping data at two densities
## result inherits from 'selectedstatistics'
scen1 <- make.scenarios(D = c(5,10), sigma = 25, g0 = 0.2, noccasions = 5)
traps1 <- make.grid() ## default 6 x 6 trap grid
tmp1 <- run.scenarios(nrepl = 5, trapset = traps1, scenarios = scen1,
fit = FALSE, seed = 123)
tmp2 <- run.scenarios(nrepl = 15, trapset = traps1, scenarios = scen1,
fit = FALSE, seed = 127)
summary(rbind(tmp1,tmp2))
summary(c(tmp1,tmp2))
Run the code above in your browser using DataLab