leslie_from_spip: Return a Leslie-like matrix from the spip parameters

In CKMRpop: Forward-in-Time Simulation and Tallying of Pairwise Relationships

Description

Here we take the survival rates for females and males and the sex ratio, as well as the annual new cohort size (assumed constant), and we make a leslie-like matrix to compute the stable age distribution.

Usage

leslie_from_spip(P, C)

Arguments

P	a named list of the spip parameters.
C	the constant size of the newborn cohort each year

Value

This function returns a list with the following components:

• stable_age_distro_fem: a vector of the expected number of females in each age group of 0 up to MaxAge-1 once the stable age distribution has been reached. Note that this corresponds to the PREKILL_CENSUS from spip. In this vector, the size of the MaxAge group is left out because this is how it is needed to be to insert into the --initial-males and --initial-females options in spip(). If you want the size of all age classes, use the output list component stable_age_distro_fem_with_max_age_class, described below.

• stable_age_distro_male: same as above, but for males.

• stable_age_distro_fem_with_max_age_class: The expected number of females from age 0 to MaxAge once the stable age distribution has been reached.

• stable_age_distro_male_with_max_age_class: same as above, but for females.

• female_leslie_matrix: The Leslie matrix implied by the spip parameters in P.

Examples

result <- leslie_from_spip(species_1_life_history, 300)

# print the result list:
result

CKMRpop documentation built on July 17, 2021, 5:07 p.m.