Simulating damage: model
Structural parameters
N_CAT: The number of damage categories each person has N_MECH: The number of mechanisms in each category
MECH_WEIGHTm: The contribution of a mechanism to a category
Fitting parameters
BASAL_M: The mean basal damage rate BASAL_SD: The standard deviation of the basal damage rate
BASAL_H: The homogeneity of basal damage rate in a single person EXP_M: The mean exponential damage rate
EXP_SD: The standard deviation of the exponential damage rate
EXP_H: The homogeneity of exponential damage rate in a single person
FATAL_M: The mean yearly challenge FATAL_SD: The standard deviation of the yearly challenge
Values set for each person at initialisation:
PB: Basal rate for the person: lognorm(BASAL_M, BASAL_SD)
PE: Exponential rate for the person: lognorm(EXP_M, EXP_SD)
MBc,m:Basal rate for each mechanism: lognorm(BASAL_M, BASAL_SD)*(1-BASAL_H) + PB*BASAL_H
MEc,m: Exponential rate for each mechanism: lognorm(EXP_M, EXP_SD)*(1-EXP_H) + PE*EXP_H
D_Mc,m : Cumulative damage for each mechanism: 0 D_Cc : Cumulative damage for each category: 0
Variables updated for each person at each time step (year):
Total damage: PD(t) = [SUM c=1..N_CAT] D_Cc(t) Damage increment: DI_Mc,m(t) = MBc,m + MEc,m*PD(t-1)
Cumulative damage: D_Mc,m(t) = DI_Mc,m(t) + D_Mc,m(t-1)
Cumulative category damage: D_Cc(t) = [SUM m=1..N_MECH] DI_Mc,m(t)
Fatality challenge: FATAL(t) = |norm(FATAL_M, FATAL_SD)|
If D_Cc(t) > FATAL(t) for any c, the person dies at age t