I am currently in the process of modeling different gene drive scenarios on snails to prevent schistosomiasis.
For my model, I ensured the fecundity of snails depends on the population density and carrying capacity to ensure that the snail population does not go over the carrying capacity. I accomplished this by using the equation that I have attached to this document. I have derived this equation through the following website:
However, from my research, the intrinsic growth rate of snails are all below the value 1. This ultimately produces a negative value, causing an error in my code. I was wondering whether this methodology is correct and whether there is an alternative method to ensure the fecundity of snails depend on carrying capacity.
As shown in the equation, I must find the intrinsic growth rate of snails. Currently all the intrinsic growth rate I have identified of snails are below 1, which ultimately produces a negative value in the result of the equation.
Thank you for your time!