pop.sim.gomp: Simulation of a population of adults with Gompertz distribution
Description
In many instances, it is useful to calculate with a population with
known parameters. To generate a population with realistic
characteristics is less obvious than it seems. We operate here
with the Gompertz distribution which provides a reasonable
approximation of human mortality for adult mortality, that is
for the ages >= 15 years. The user has to specify
either the parameter b or the modal age M. The modal age M is
particular useful as it provides an intuitive understanding of
the resulting age distribution. In both instances, the second
parameter a is generated by the regression formula found by
Sasaki and Kondo 2016. If neither is given, a population
with random parameters realistic for pre-modern times is generated.
Usage
pop.sim.gomp(n, b = NULL, M = NULL, start_age = 15, max_age = 100)
Value
A list of two data.frames with the following items:
First data.frame
N: Number of individuals.
b: Gompertz parameter controlling mortality.
M: Modal age.
a: Gompertz parameter controlling hazard of the
youngest age group.
Second data.frame
ind: ID of individuals.
age: Simulated absolute age.
Arguments
n
number of individuals to be simulated.
b
numeric, optional. Gompertz parameter controlling the
level of mortality.
M
numeric, optional. Modal age M.
start_age
numeric. Start age, default: 15 years.
max_age
numeric. Maximal age, to avoid unlikely centenaries,
default: 100 years.