Model-Based Estimation of Expected Time to Cholera Extinction in Lusaka, Zambia.

The developing world has been facing a significant health issue due to cholera as an endemic communicable disease. Lusaka was Zambia's worst affected province, with 5414 reported cases of cholera during the outbreak from late October 2017 to May 12, 2018. To explore the epidemiological characteristics associated with the outbreak, we fitted weekly reported cholera cases with a compartmental disease model that incorporates two transmission routes, namely environment-to-human and human-to-human. Estimates of the basic reproduction number show that both transmission modes contributed almost equally during the first wave. In contrast, the environment-to-human transmission appears to be mostly dominating factor for the second wave. Our study finds that a massive abundance of environmental vibrio's with a huge reduction in water sanitation efficacy triggered the secondary wave. To estimate the expected time to extinction (ETE) of cholera, we formulate the stochastic version of our model and find that cholera can last up to 6.5-7 years in Lusaka if any further outbreak occurs at a later time. Results indicate that a considerable amount of attention is to be paid to sanitation and vaccination programs in order to reduce the severity of the disease and to eradicate cholera from the community in Lusaka.

A noble extended stochastic logistic model for cell proliferation with density-dependent parameters.

Cell proliferation often experiences a density-dependent intrinsic proliferation rate (IPR) and negative feedback from growth-inhibiting molecules in culture media. The lack of flexible models with explanatory parameters fails to capture such a proliferation mechanism. We propose an extended logistic growth law with the density-dependent IPR and additional negative feedback. The extended parameters of the proposed model can be interpreted as density-dependent cell-cell cooperation and negative feedback on cell proliferation. Moreover, we incorporate further density regulation for flexibility in the model through environmental resistance on cells. The proposed growth law has similarities with the strong Allee model and harvesting phenomenon. We also develop the stochastic analog of the deterministic model by representing possible heterogeneity in growth-inhibiting molecules and environmental perturbation of the culture setup as correlated multiplicative and additive noises. The model provides a conditional maximum sustainable stable cell density (MSSCD) and a new fitness measure for proliferative cells. The proposed model shows superiority to the logistic law after fitting to real cell culture datasets. We illustrate both conditional MSSCD and the new cell fitness for a range of parameters. The cell density distributions reveal the chance of overproliferation, underproliferation, or decay for different parameter sets under the deterministic and stochastic setups.

An extended stochastic Allee model with harvesting and the risk of extinction of the herring population.

Overexploitation of commercially beneficial fish is a serious ecological problem around the world. The growth profiles of most of the species are likely to follow density regulated theta-logistic model irrespective of any taxonomy group [Sibly et al., Science, 2005]. Rapid depletion of population size may cause reduced fitness, and the species is exposed to Allee phenomena. Here sustainability is addressed by modelling the herring population as a stochastic process and computing the probability of extinction and expected time to extinction. The models incorporate an Allee effect, crowding effect which reduce birth and death rates at large populations, and two possible choices of harvesting models viz. linear harvesting and nonlinear harvesting. A seminal attempt is made by Saha [Saha et al., Ecol. Model, 2013] for this economically beneficial fish, but ignored the vital phenomena of harvesting. Moreover, in this model, the demographic stochasticity is introduced through the white-noise term, which has certain limitations when harvesting is introduced into the system. White noise is appropriate for such a system where immigration and emigration are allowed, but a harvesting model is rational for a closed system. The demographic stochasticity is introduced by replacing an ordinary differential equation model with a stochastic differential equation model, where the instantaneous variance in the SDE is derived directly from the birth and death rates of a birth-death process. The modelling parameters are fit to data of the herring populations collected from Global Population Dynamics Database (GPDD), and the risk of extinction of each population is computed under different harvesting protocols. A threshold for handling times is computed beneath which the risk of extinction is high. This is proposed as a recommendation to management for sustainable harvesting.

Chance of extinction of populations in food chain model under demographic stochasticity.

The extinction of different species from the earth is increasing at an alarming rate. So, assessment of probability of extinction of different important species in our ecosystem could help us to take proper conservation policy for those population whose chance of extinction is high. In this paper a method is developed to find the probability of extinction of populations in a general n-trophic food chain model under demographic stochasticity. The birth-death process is used to incorporate the demographic stochasticity and the necessary mathematical expressions are obtained. The theoretical finding is validated by numerical simulation for a two dimensional predator-prey system.

Mathematical analysis of a power-law form time dependent vector-borne disease transmission model.

In the last few years, fractional order derivatives have been used in epidemiology to capture the memory phenomena. However, these models do not have proper biological justification in most of the cases and lack a derivation from a stochastic process. In this present manuscript, using theory of a stochastic process, we derived a general time dependent single strain vector borne disease model. It is shown that under certain choice of time dependent transmission kernel this model can be converted into the classical integer order system. When the time-dependent transmission follows a power law form, we showed that the model converted into a vector borne disease model with fractional order transmission. We explicitly derived the disease-free and endemic equilibrium of this new fractional order vector borne disease model. Using mathematical properties of nonlinear Volterra type integral equation it is shown that the unique disease-free state is globally asymptotically stable under certain condition. We define a threshold quantity which is epidemiologically known as the basic reproduction number (R0). It is shown that if R0 > 1, then the derived fractional order model has a unique endemic equilibrium. We analytically derived the condition for the local stability of the endemic equilibrium. To test the model capability to capture real epidemic, we calibrated our newly proposed model to weekly dengue incidence data of San Juan, Puerto Rico for the time period 30th April 1994 to 23rd April 1995. We estimated several parameters, including the order of the fractional derivative of the proposed model using aforesaid data. It is shown that our proposed fractional order model can nicely capture real epidemic.