Extract Estimates From Fitted Models

Description

If simulations have been saved from run.scenarios as fitted secr models it is necessary to use one of these functions to extract estimates for later summarization.

Usage

1
2
3
4
5
6
7
8
9
## S3 method for class 'fittedmodels'
predict(object, ...)

## S3 method for class 'fittedmodels'
coef(object, ...)

derived.SL(object, ...)

regionN.SL(object, ...)

Arguments

object

fitted model simulation output from run.scenarios

...

other arguments passed to predict, coef, derived or region.N

Details

These functions are used when output from run.scenarios has been saved as fitted models. derived.SL and regionN.SL require a full fit (including the design object) whereas a trimmed model is sufficient for predict and coef.

derived.SL is used to compute the Horvitz-Thompson-like estimate of density when secr.fit has been used with CL = TRUE; it is roughly equivalent to predict.

regionN.SL predicts the realised number (R.N) or expected number (E.N) in a masked area. When detector layouts and/or sigma vary, the masked area will also vary (arbitrarily, depending on the buffer argument ‘xsigma’) unless a mask is provided by the user; this may be done either in run.scenarios or in regionN.SL.

Value

An object with class (‘estimatetables’, ‘secrdesign’, ‘list’) with appropriate outputtype (‘predicted’, ‘coef’, ‘derived’, ‘regionN’; sSee also run.scenarios).

See Also

run.scenarios

Examples

1
2
3
4
5
6
7
8
## using nrepl = 2 just for checking
scen1 <- make.scenarios(D = c(5,10), sigma = 25, g0 = 0.2)
traps1 <- make.grid()  ## default 6 x 6 grid of multi-catch traps
tmp1 <- run.scenarios(nrepl = 2, trapset = traps1, scenarios = scen1,
    fit = TRUE, extractfn = trim)
tmp2 <- predict(tmp1)
tmp3 <- select.stats(tmp2, 'D', c('estimate','RB','RSE'))
summary(tmp3)