1,578 research outputs found

    Analysis of a heroin epidemic model with saturated treatment function

    Get PDF
    A mathematical model is developed that examines how heroin addiction spreads in society. The model is formulated to take into account the treatment of heroin users by incorporating a realistic functional form that "saturates" representing the limited availability of treatment. Bifurcation analysis reveals that the model has an intrinsic backward bifurcation whenever the saturation parameter is larger than a fixed threshold. We are particularly interested in studying the model's global stability. In the absence of backward bifurcations, Lyapunov functions can often be found and used to prove global stability. However, in the presence of backward bifurcations, such Lyapunov functions may not exist or may be difficult to construct. We make use of the geometric approach to global stability to derive a condition that ensures that the system is globally asymptotically stable. Numerical simulations are also presented to give a more complete representation of the model dynamics. Sensitivity analysis performed by Latin hypercube sampling (LHS) suggests that the effective contact rate in the population, the relapse rate of heroin users undergoing treatment, and the extent of saturation of heroin users are mechanisms fuelling heroin epidemic proliferation

    Optimal control of a heroin epidemic mathematical model

    Get PDF
    A heroin epidemic mathematical model with prevention information and treatment, as control interventions, is analyzed, assuming that an individual's behavioral response depends on the spreading of information about the effects of heroin. Such information creates awareness, which helps individuals to participate in preventive education and self-protective schemes with additional efforts. We prove that the basic reproduction number is the threshold of local stability of a drug-free and endemic equilibrium. Then, we formulate an optimal control problem to minimize the total number of drug users and the cost associated with prevention education measures and treatment. We prove existence of an optimal control and derive its characterization through Pontryagin's maximum principle. The resulting optimality system is solved numerically. We observe that among all possible strategies, the most effective and cost-less is to implement both control policies.publishe

    Global analysis of a multi-group SIR epidemic model with nonlinear incidence rates and distributed moving delays between patches

    Get PDF
    In this paper, applying Lyapunov functional approach, we establish sufficient conditions under which each equilibrium is globally asymptotically stable for a class of multi-group SIR epidemic models. The incidence rate is given by nonlinear incidence rates and distributed delays incorporating not only an exchange of individuals between patches through migration but also cross patch infection between different groups. We show that nonlinear incidence rates and distributed delays have no influence on the global stability, but patch structure has. Moreover, the present results generalize known results on the global stability of a heroin model with two delays considered in the recent literatures. We also offer new techniques to prove the boundedness of the solutions, the existence of the endemic equilibrium and permanence of the model

    A survey on Lyapunov functions for epidemic compartmental models

    Get PDF
    In this survey, we propose an overview on Lyapunov functions for a variety of compartmental models in epidemiology. We exhibit the most widely employed functions, and provide a commentary on their use. Our aim is to provide a comprehensive starting point to readers who are attempting to prove global stability of systems of ODEs. The focus is on mathematical epidemiology, however some of the functions and strategies presented in this paper can be adapted to a wider variety of models, such as prey–predator or rumor spreading

    Survival analysis and probability density function of switching heroin model

    Get PDF
    We study a switching heroin epidemic model in this paper, in which the switching of supply of heroin occurs due to the flowering period and fruiting period of opium poppy plants. Precisely, we give three equations to represent the dynamics of the susceptible, the dynamics of the untreated drug addicts and the dynamics of the drug addicts under treatment, respectively, within a local population, and the coefficients of each equation are functions of Markov chains taking values in a finite state space. The first concern is to prove the existence and uniqueness of a global positive solution to the switching model. Then, the survival dynamics including the extinction and persistence of the untreated drug addicts under some moderate conditions are derived. The corresponding numerical simulations reveal that the densities of sample paths depend on regime switching, and larger intensities of the white noises yield earlier times for extinction of the untreated drug addicts. Especially, when the switching model degenerates to the constant model, we show the existence of the positive equilibrium point under moderate conditions, and we give the expression of the probability density function around the positive equilibrium point

    Stability analysis of drinking epidemic models and investigation of optimal treatment strategy

    Get PDF
    In this research we investigate a class of drinking epidemic models, namely the SPARS type models. The basic reproduction number is derived, and the system dynamical behaviours are investigated for both drinking free equilibrium and drinking persistent equilibrium. The purpose is to determine the long term optimal treatment method and the optimal short period vaccination strategy for controlling the population of the periodic drinkers and alcoholics

    Uncertainty quantification in dynamical models. An application to cocaine consumption in Spain

    Full text link
    [EN] The present Ph.D. Thesis considers epidemiological mathematical models based on ordinary differential equations and shows its application to understand the cocaine consumption epidemic in Spain. Three mathematical models are presented to predict the evolution of the epidemic in the near future in order to select the model that best reflects the data. By the results obtained for the selected model, if there are not changes in cocaine consumption policies or in the economic environment, the cocaine consumption will increase in Spain over the next few years. Furthermore, we use different techniques to estimate 95% confidence intervals and, consequently, quantify the uncertainty in the predictions. In addition, using several techniques, we conducted a model sensitivity analysis to determine which parameters are those that most influence the cocaine consumption in Spain. These analysis reveal that prevention actions on cocaine consumer population can be the most effective strategy to control this trend.[ES] La presente Tesis considera modelos matemáticos epidemiológicos basados en ecuaciones diferenciales ordinarias y muestra su aplicación para entender la epidemia del consumo de cocaína en España. Se presentan tres modelos matemáticos para predecir la evolución de dicha epidemia en un futuro próximo, con el objetivo de seleccionar el modelo que mejor refleja los datos. Por los resultados obtenidos para el modelo seleccionado, si no hay cambios en las políticas del consumo de cocaína ni en el ámbito económico, el consumo de cocaína aumentará en los próximos años. Además, utilizamos diferentes técnicas para estimar los intervalos de confianza al 95% y, de esta forma, cuantificar la incertidumbre en las predicciones. Finalmente, utilizando diferentes técnicas, hemos realizado un análisis de sensibilidad para determinar qué parámetros son los que más influyen en el consumo de cocaína. Estos análisis revelan que las acciones de prevención sobre la población de consumidores de cocaína pueden ser la estrategia más efectiva para controlar esta tendencia.[CA] La present Tesi considera models matemàtics epidemiològics basats en equacions diferencials ordinàries i mostra la seua aplicació per a entendre l'epidèmia del consum de cocaïna en Espanya. Es presenten tres models matemàtics per a predir l'evolució d'aquesta epidèmia en un futur pròxim, amb l'objectiu de seleccionar el model que millor reflecteix les dades. Pels resultats obtinguts per al model seleccionat, si no hi ha canvis en les polítiques de consum de cocaïna ni en l'àmbit econòmic, el consum de cocaïna augmentarà en els pròxims anys. A més, utilitzem diferents tècniques per a estimar els intervals de confiança al 95% i, d'aquesta manera, quantificar la incertesa en les prediccions. Finalment, utilitzant diferents tècniques, hem realitzat un anàlisi de sensibilitat per a determinar quins paràmetres són els que més influencien el consum de cocaïna. Aquestos anàlisis revelen que les accions de prevenció en la població de consumidors de cocaïna poden ser l'estratègia més efectiva per a controlar aquesta tendència.Rubio Monzó, M. (2015). Uncertainty quantification in dynamical models. An application to cocaine consumption in Spain [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/55844TESI

    Cross immunity protection and antibody-dependent enhancement in a distributed delay dynamic model

    Get PDF
    Dengue fever is endemic in tropical and subtropical countries, and certain important features of the spread of dengue fever continue to pose challenges for mathematical modelling. Here we propose a system of integro-differential equations (IDE) to study the disease transmission dynamics that involve multi-serotypes and cross immunity. Our main objective is to incorporate and analyze the effect of a general time delay term describing acquired cross immunity protection and the effect of antibody-dependent enhancement (ADE), both characteristics of Dengue fever. We perform qualitative analysis of the model and obtain results to show the stability of the epidemiologically important steady solutions that are completely determined by the basic reproduction number and the invasion reproduction number. We establish the global dynamics by constructing a suitable Lyapunov functional. We also conduct some numerical experiments to illustrate bifurcation structures, indicating the occurrence of periodic oscillations for a specific range of values of a key parameter representing ADE.Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brazil (CAPES) - Finance Code 001 LIAM - Laboratory for Industrial and Applied Mathematics, Department of Mathematics and Statistics, York University-CA

    The 1950s “War on Narcotics”: Harry Anslinger, The Federal Bureau of Narcotics, and Senator Price Daniel’s Probe

    Get PDF
    corecore