EMOD-Hub / EMOD

Source files for building the IDM EMOD disease transmission model.
MIT License
94 stars 41 forks source link

Add Blood_Feed_Mortality #45

Closed Bridenbecker closed 3 weeks ago

Bridenbecker commented 4 weeks ago

Prashanth is working on an article with Nicholi and needs to add an ability to have some likelihood that certain vectors die due to making a blood feed. There are some gene drives that have a fitness cost where the vector is more likely to die having made a blood feed.

The plan is to add a new probability check to the vector life cycle as a post feed mortality. It will be after the current "died post feed" or included in it. The current "died post feed" involves post IRS killing but this is supposed to be after that.

We will add a new parameter to Vector_Species_Params

"Blood_Feed_Mortality" : {
    "Default_Probability": 0.0,
    "Genetic_Probabilities": [
        {
            "Allele_Combinations" : [
                [ "a1", "a1" ]
            ],
            "Probability": 1.0
        },
        {
            "Allele_Combinations" : [
                [ "a1", "a0" ]
            ],
            "Probability": 0.8
        }
    ]
}

I propose that in VectorInterventionsContainer we change pDiePostFeeding such that it includes both dying from IRS post feeding and from the blood meal. I'm assuming that dying from the blood meal happens after dying from the post IRS. That is, the mosquito feeds on person, lands on IRS treated wall to rest, leaves house, and dies from blood meal. Hence, we currently have:

pDiePostFeeding = not_housing_IRSprefeed_ADIH_net_indrep_dieduringfeeding * p_kill_IRSpostfeed_effective;

should change to:

GeneticProbability p_kill_post_feed = p_kill_IRSpostfeed_effective + (1.0 - p_kill_IRSpostfeed_effective) * p_die_blood_meal;
pDiePostFeeding = not_housing_IRSprefeed_ADIH_net_indrep_dieduringfeeding * p_kill_post_feed;

where p_die_blood_meal is defined by the JSON above. There will be other changes required to use p_kill_post_feed

NOTE: The GeneticProbabilities created for each species will need to be combined into one GeneticProbability that is owned by VectorParameters. This combined value will be used in VectorInterventionsContainer.