GenEst (version 1.2.0)

estg: Estimate all carcass-level detection rates and arrival intervals

Description

Estimate g values and arrival intervals for a set of carcasses from fitted pk and cp models and search data

Usage

estg(data_CO, COdate, data_SS, SSdate = NULL, model_SE, model_CP,
  sizeCol = NULL, unitCol = NULL, nsim = 1000, max_intervals = 8,
  seed_SE = NULL, seed_CP = NULL, seed_g = NULL)

Arguments

data_CO

Carcass Observation data

COdate

Column name for the date found data

data_SS

Search Schedule data

SSdate

Column name for the date searched data. Optional. If not provided, estg will try to find the SSdate among the columns in data_SS. See prepSS.

model_SE

Searcher Efficiency model (or list of models if there are multiple size classes)

model_CP

Carcass Persistence model (or list of models if there are multiple size classes)

sizeCol

Name of column in data_CO where the size classes are recorded. Optional. If not provided, no distinctions are made among sizes. sizeCol not only identifies what the name of the size segregating class

unitCol

Column name for the unit indicator

nsim

the number of simulation draws

max_intervals

maximum number of arrival interval intervals to consider for each carcass. Optional. Limiting the number of search intervals can greatly increase the speed of calculations with only a slight reduction in accuracy in most cases.

seed_SE

seed for random draws of the SE model

seed_CP

seed for random draws of the CP model

seed_g

seed for random draws of the gs

Value

list of [1] g estimates (ghat) and [2] arrival interval estimates (Aj) for each of the carcasses. The row names of the Aj matrix are the names of the units where each carcass was found.

Examples

# NOT RUN {
 data(mock)
 model_SE <- pkm(formula_p = p ~ HabitatType, formula_k = k ~ 1,
              data = mock$SE)
 model_CP <- cpm(formula_l = l ~ Visibility, formula_s = s ~ Visibility, 
               data = mock$CP, dist = "weibull",
               left = "LastPresentDecimalDays", 
               right = "FirstAbsentDecimalDays"
             )
 ghat <- estg(data_CO = mock$CO, COdate = "DateFound",  data_SS = mock$SS,
           model_SE = model_SE, model_CP = model_CP, unitCol = "Unit", nsim = 100)

# }