"""
Data for the EC process.
This data was obtained and used in the following study.

Alblas, Maaike et al. “Cost-effectiveness of prophylactic hysterectomy in first-degree female
relatives with Lynch syndrome of patients diagnosed with colorectal cancer in the United States:
a microsimulation study.” Cancer medicine vol. 10,19 (2021): 6835-6844. doi:10.1002/cam4.4080
"""

hazard = 30.

age_factors = {
    "without_atypia": [
        (0, 0.0),
        (30, 0.0000088 * 6.142857143),
        (40, 0.000455 * 6.142857143),
        (50, 0.000667 * 6.142857143),
        (60, 0.004254 * 6.142857143),
        (70, 0.000625 * 6.142857143),
        (80, 0.8 * 6.142857143),
        (100, 0.8 * 6.142857143),
    ],
    "with_atypia": [
        (0, 0.0),
        (30, 0.0000088),
        (40, 0.000455),
        (50, 0.000667),
        (60, 0.004254),
        (70, 0.000625),
        (80, 0.8),
        (100, 0.8),
    ]
}

dwell_hyperplasia = {
    "without_atypia": {"shape": .75, "mean": 114.40339},
    "with_atypia": {"shape": .75, "mean": 7.766915},
}

dwell_preclinical = {"shape": .75, "mean": 5.}

cure = [
    (0, 0.780),
    (20, 0.780),
    (20, 0.811),
    (25, 0.811),
    (25, 0.873),
    (30, 0.873),
    (30, 0.845),
    (35, 0.845),
    (35, 0.815),
    (40, 0.815),
    (40, 0.817),
    (45, 0.817),
    (45, 0.812),
    (50, 0.812),
    (50, 0.833),
    (55, 0.833),
    (55, 0.831),
    (60, 0.831),
    (60, 0.806),
    (65, 0.806),
    (65, 0.761),
    (70, 0.761),
    (70, 0.712),
    (75, 0.712),
    (75, 0.661),
    (80, 0.661),
    (80, 0.442),
    (100, 0.442),
]

time_to_death = [
    (0.000000000, 0),
    (0.309278351, 1),
    (0.515463918, 2),
    (0.659793814, 3),
    (0.752577320, 4),
    (0.804123711, 5),
    (0.840206186, 6),
    (0.865979381, 7),
    (0.891752577, 8),
    (0.912371134, 9),
    (0.927835052, 10),
    (0.927835052, 11),
    (0.927835052, 12),
    (0.927835052, 13),
    (0.938144330, 14),
    (0.948453608, 15),
    (0.963917526, 16),
    (0.963917526, 17),
    (0.963917526, 18),
    (0.974226804, 19),
    (1.000000000, 20),
]

# Note ``onset_method`` can only take values ``first`` or ``before_oc``.
# See the EC process tutorial for more information.
onset_method = "before_oc"