Create summary table of demographic rates from survival and recruitment surveys

Create summary table of demographic rates from survival and recruitment surveys

Usage

estimateBayesianRates(
  surv_data,
  recruit_data,
  N0 = NA,
  disturbance = NULL,
  priors = NULL,
  shiny_progress = FALSE,
  return_mcmc = FALSE,
  i18n = NULL,
  niters = formals(bboutools::bb_fit_survival)$niters,
  nthin = formals(bboutools::bb_fit_survival)$nthin,
  ...
)

Arguments

surv_data

dataframe. Survival data in bboudata format

recruit_data

dataframe. Recruitment data in bboudata format

N0

dataframe. Optional. Initial population estimates, required columns are PopulationName and N0

disturbance

dataframe. Optional. If provided, fit a Beta model that includes disturbance covariates.

priors

list. Optional. If disturbance is NA, this should be list(priors_survival=c(...),priors_recruitment=c(...)); see bboutools::bb_priors_survival and bboutools::bb_priors_recruitment for details. If disturbance is not NA, see betaNationalPriors for details.

shiny_progress

logical. Should shiny progress bar be updated. Only set to TRUE if using in an app.

return_mcmc

boolean. If TRUE return fitted survival and recruitment models. Default FALSE.

niters

integer. The number of iterations per chain after thinning and burn-in.

nthin

integer. The number of the thinning rate.

...

Other parameters passed on to bboutools::bb_fit_survival and bboutools::bb_fit_recruitment.

Value

If return_mcmc is TRUE then a list with results and fitted models, if FALSE just the results summaries are returned.

See also

Caribou demography functions: bayesianScenariosWorkflow(), bayesianTrajectoryWorkflow(), betaNationalPriors(), caribouPopGrowth(), compareTrajectories(), compositionBiasCorrection(), convertTrajectories(), dataFromSheets(), demographicProjectionApp(), demographyDefaults(), disturbanceDefaults(), estimateNationalRate(), getNationalCoefficients(), getScenarioDefaults(), monitoringDefaults(), nationalTrajectoryDefaults(), plotCompareTrajectories(), plotSurvivalSeries(), plotTrajectories(), popGrowthTableJohnsonECCC, simulateObservations(), timeDefaults(), trajectoriesFromBayesian(), trajectoriesFromNational(), trajectoriesFromSummary(), trajectoriesFromSummaryForApp()

Examples

s_data <- rbind(bboudata::bbousurv_a, bboudata::bbousurv_b)
r_data <- rbind(bboudata::bbourecruit_a, bboudata::bbourecruit_b)
estimateBayesianRates(s_data, r_data, N0 = 500)
#>   PopulationName     R_bar       R_sd R_iv_mean R_iv_shape R_bar_lower
#> 1              A 0.1900485 0.08182184 0.3154416   25.83285   0.1662411
#> 2              B 0.2032712 0.10050771 0.3154416   25.83285   0.1722409
#>   R_bar_upper     S_bar     S_sd S_iv_mean S_iv_shape S_bar_lower S_bar_upper
#> 1   0.2146784 0.8809044 0.223287 0.4952266   13.08777   0.8282027   0.9202089
#> 2   0.2362427 0.9069938 0.307349 0.4952266   13.08777   0.8468059   0.9498290
#>    N0 nCollarYears nSurvYears nCowsAllYears nRecruitYears
#> 1 500          900         31          2353            27
#> 2 500          519         18          2001            15