Simulate Boreal Caribou data from survival rates, fecundity rates and sampling parameters.

bbs_simulate_caribou(
  survival,
  fecundity,
  nsims = 1,
  adult_females = 500,
  proportion_adult_female = 0.65,
  proportion_yearling_female = 0.5,
  probability_unsexed_adult_female = 0,
  probability_unsexed_adult_male = 0,
  month_composition = 9L,
  group_size = 5,
  group_coverage = 0.3,
  group_min_size = 2,
  group_max_proportion = 1,
  collared_adult_females = 30,
  month_collar = 1L,
  probability_uncertain_mortality = 0,
  probability_uncertain_survival = 0,
  population_name = "A"
)

Arguments

survival

A list of the Caribou survival rates (output of bbs_survival_caribou())

fecundity

A list of the Caribou fecundity rates (output of bbs_fecundity_caribou())

nsims

A whole positive number of the number of simulations.

adult_females

A number of the initial number of adult females in the population.

proportion_adult_female

A number between 0 and 1 of the proportion of adults that are female.

proportion_yearling_female

A number between 0 and 1 of the proportion of yearlings that are female.

probability_unsexed_adult_female

A number between 0 and 1 of the probability that an adult female is unsexed.

probability_unsexed_adult_male

A number between 0 and 1 of the probability that an adult male is unsexed.

month_composition

A whole number between 1 and 12 of the month that composition surveys take place, relative to the start of the biological year.

group_size

A whole number of the average group size. Group sizes are drawn from a poisson distribution with lambda of group_size.

group_coverage

A number between 0 and 1 of the proportion of groups sampled.

group_min_size

A whole positive number of the minimum group size.

group_max_proportion

A number between 0 and 1 of the maximum group size as proportion of the total population.

collared_adult_females

A whole positive number of the number of collared adult females. The number of collared adult females is 'topped up' each year at month_collar.

month_collar

A whole number between 1 and 12 of the collaring month.

probability_uncertain_mortality

A number between 0 and 1 of the probability of uncertain mortality.

probability_uncertain_survival

A number between 0 and 1 of the probability of uncertain survival.

population_name

A string of the population name. This does not affect simulation but can be used as a unique identifier.

Value

A list of three tibbles. The first named survival has columns Year, Month, StartTotal, MortalitiesCertain and MortalitiesUncertain. The second named recruitment has columns Year, Month, Cows, Bulls, UnknownAdults, Yearlings and Calves. And the final named abundance has columns Year, Month, Stage, Abundance.

Details

See bbs_survival_caribou() and bbs_fecundity_caribou() for generating survival and fecundity rates. See bbs_population_caribou() for details on how population is simulated. See bbs_population_groups_survey() for details on how groups are assigned for each composition survey. See bbs_survival_collared() for details on how survival of collared adult female is determined.

Survival and recrutiment data.frames generated are formatted to be used as input data to bboutools::bb_fit_survival() and bboutools::bb_fit_recruitment(), respectively.

Examples

survival <- bbs_survival_caribou(0.84)
fecundity <- bbs_fecundity_caribou(0.7)
caribou <- bbs_simulate_caribou(survival, fecundity = fecundity, nsims = 2)