est_g0 {eoa3}

estimate g from fitted pk and cp models and search schedule

Description

Given a fitted pk model (from est_pk0), a fitted cp model (from est_cp0 or a survreg object), and a search schedule, estimate detection probability.

Usage

est_g0(pkmodel, cpmodel, days, a = NULL, v = NULL, ...)

Arguments

pkmodel

fitted pk model (from est_pk0)

cpmodel

fitted cp model (from est_cp0 or survreg object)

days

vector of days since searches begin (days[1] == 0)

a

fraction of carcasses arriving in the area searched

v

fraction of carcasses arriving in the period spanned by the monitoring

...

additional arguments (ignored)

Value

list of parameters for a beta distributions fit to the vectors of for \hat{g} for the searched area within the period monitored ($BabRaw), for the whole site within the period monitored ($Bab), and for the whole site extrapolated to the whole year ($BabAnn). In addition, the models and parameters that went into the estimate are included as well (pkmodel, cpmodel, a, v).

Examples

 pkmodel <- est_pk0(pkdata = pkdata0)
 cpmodel <- est_cp0(cpdata = cpdata0, dist = "weibull")
 ghat <- est_g0(pkmodel = pkmodel, cpmodel = cpmodel, days = days0, a = 0.4, v = 0.75)
 summary(ghat)