Primarily intended as an internal function, get_nextEvent randomly simulates an individual's next life event given their current age, disease status, and relative-risk of disease.
get_nextEvent(current_age, disease_status, RV_status, hazard_rates, GRR,
carrier_prob, lambda_birth, birth_range, fert = 1)Numeric. The individual's current age.
Numeric. The individual's disease status, where disease_status = 1 if individual has experienced disease onset, otherwise disease_status = 0.
Numeric. RV_status = TRUE if the individual is a carrier of a rare variant that increases disease suseptibility, and FALSE otherwise.
An object of class hazard, created by hazard.
Numeric. The genetic relative-risk of disease, i.e. the relative-risk of disease for individuals who carry at least one copy of the causal variant. Note: When simulating diseases with multiple subtypes GRR must contain one entry for each simulated subtype. See details.
Numeric. The carrier probability for all causal variants with relative-risk of disease GRR. By default, carrier_prob = 0.002
Numeric. The individual's birth rate.
Numeric vector of length 2. The minimum and maximum allowable ages, in years, between which this individual may reproduce.
Numeric. A constant used to rescale the fertility rate after disease-onset. By default, fert = 1.
A named matrix. The number of years until the next life event, named by event type. See Details.
Given their current age, get_nextEvent randomly simulates an individual's next life event by generating waiting times to reproduction, onset, and death. The event with the shortest waiting time is chosen as the next life event.
We assume that, given an individual's current age, their time to disease onset is the waiting time in a non-homogeneous Poisson process with an age-specific hazard rate that follows a proportional hazards model. In this model, individuals who have NOT inherited the rare variant experience disease onset according to the baseline (or population) hazard rate of disease. On the other hand, individuals who have inherited the rare variant are assumed to have an increased risk of disease onset relative to those who have inherited it. The user is expected to supply the population age-specific hazard rate of disease, the relative-risk of disease for genetic cases, and the population allele frequency of the rare variant. We impose the restriction that individuals may only experience disease onset once, and remain affected from that point on.
We assume that, given an individual's current age, their time to death is the waiting time in a non-homogeneous Poisson process with age-specific hazard rate determined by their affection status. We assume that unaffected individuals experience death according to the age-specific hazard rate for death in the unaffected population. If the disease of interest is sufficiently rare, the user may instead choose to substitute the population age-specific hazard rate for death in the general population. We assume that affected individuals experience death according to the age-specific hazard rate for death in the affected population. The user is expected to supply both of these age-specific hazard rates.
We assume that, given an individual's current age, their time to reproduction is the waiting time in a homogeneous Poisson process. That is, we assume that individuals reproduce at uniform rate during their reproductive years. For example, one's reproductive years may span from age 18 to age 45. We do not allow for offspring to be produced outside of an individual's birth_range.
If get_nextEvent returns the waiting time to the next life event, named for event type. The possible event types are as follows:
"Child" a reproductive event, i.e. creation of offspring
"Onset" disease onset event,
"Death" death event