![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
A new social epidemic model to depict alcoholism with media coverage is proposed in this paper. Some fundamental properties of the model including existence and positivity as well as boundedness of equilibria are investigated. Stability of all equilibria are studied. The existence of the optimal control pair and mathematical expressions of optimal control are obtained by Pontryagin's maximum principle. Numerical simulations are also performed to illustrate our results. Our results show that media coverage is an effective measure to quit drinking.
1. Introduction
The mathematical modelling of topics involving alcoholism is becoming an important area [Citation9, Citation15, Citation22]. There is a strong evidence that alcoholism is associated with many health problems and the global burden of disease [Citation22]. Alcohol consumption has negative effects on all organs of the body. Approximately, worldwide 3.6% of cancers derive from chronic alcohol drinking, including those of the upper aerodigestive tract, the liver, the colorectum and the breast [Citation9]. It has observed a significant indirect relationship between the degree of working memory decline induced by alcohol and adverse consequences of alcohol use [Citation15]. There are many valid ways to quit drinking, such as gradually reducing the amount of alcohol or taking tablets and so on. A number of mathematical models for alcoholism has been established and analysed [Citation19, Citation24, Citation26]. As early as 1953, Stratus and Bacon began to research classical problem of college students drinking [Citation26]. Sanchez et al. [Citation24] divided the people into three classes, namely susceptible, alcoholics, alcoholics undergoing treatment. They assumed that each of the three parts had the same death rate and took into account relapse, but they ignored the death rate due to heavy alcohol consumption. Mulone and Straughan [Citation19] studied a model for binge drinking taking into account admitting and non-admitting drinkers. For the other mathematical models about drinking, smoking or epidemics, we refer to [Citation3–5, Citation10, Citation11, Citation20, Citation25, Citation29–31] and the references cited therein.
Nowadays, media can play an important role in the spread and control of infectious disease. Media usually can cut down the opportunity and probability of contact transmission of disease. Liu et al. [Citation17] constructed a new model to explain the mechanism for outbreaks of emerging infectious diseases due to the psychological effect of the reported numbers of infectious and hospitalized individuals. Samanta et al. [Citation23] built a mathematical model to study the impact of awareness programs by media on the emergence of infectious diseases. Liu and Cui [Citation16] set up a model to investigate the impact of media coverage and controlling of infectious disease in a given region. Ma et al. [Citation18] introduced a dynamic alcohol consumption model with awareness programs and time delay and captured the effects of awareness programs and time delay in controlling the alcohol problems.
Different control strategies, such as education campaigns, resource allocation, vaccination and treatment, are used in many epidemic models [Citation1, Citation2, Citation6, Citation13]. Zaman [Citation32] considered two possible control variables in the form of education and treatment campaigns oriented to decrease smoking. Wang et al. [Citation28] formulated an alcohol quitting model in which they consider the impact of distributed time delay between contact and infection process by characterizing dynamic nature of alcoholism behaviours, and they generalized the infection rate to the general case. Furthermore they considered two control strategies,that is education and treatment.
The structure of this paper is organized as follows. In Section 2, a more realistic alcoholism model with media coverage is constructed. Motivated by Wang et al. [Citation27], we use saturation function (Holling type-II) to reflect the influence of media on the alcoholics. In Section 3, the basic reproduction number and existence of alcoholism equilibrium are derived. In Section 4, stability analysis of equilibria of model is investigated by using Lyapunov function and Routh–Hurwitz criterion. In Section 5, optimal control strategies by using the classical method of Pontryagin's maximum principle are discussed. In Section 6, some sensitivity analysis and numerical simulations are presented. Some conclusions and discussions are also given.
2. Model formulation
The total population is divided into three compartments in our model, susceptible individuals who do not drink or drink only moderately, denoted by ; alcoholics who binge drink and affect the physical health seriously, denoted by
, the quitting drinkers who permanently quit drinking, denoted by
. It is supposed that the cumulative density of susceptible individuals and quitting drinkers driven by media coverage at time t is
. The growth rate of
is proportional to the number of mortalities induced by alcoholics. The total number of population at time t is given by
(1)
(1)
The population flow among those compartments is shown in Figure .
The mathematical representation of the model consists of a system of nonlinear differential equations with four state variables. The complete model is as follows:
(2)
(2) where b is the recruitment rate of susceptible individuals who do not drink or drink only moderately, β is the transmission coefficient form the susceptible individuals to the alcoholics. Large numbers of media coverage causes less interaction between susceptible aware and alcoholic populations, a mathematical form of this assumption follows the Holling type-II functional form with half-saturating constant
(see [Citation27]). Similarly, we use the Holling type-II functional form with half-saturating constant
to reflect the effect of media coverage on alcoholics. q is the rate coefficient of the individuals who permanently quit drinking, μ is the natural death rate of the population, d is the death rate due to heavy alcohol consumption, m is the implementation rate of media coverage, γ is the depletion rate of media coverage due to ineffectiveness, social problems and similar factors.
Lemma 2.1
If
the solutions
of Equation (Equation2
(2)
(2) ) are positive for all t>0.
Proof.
If the conclusion dose not satisfy, then at least one of ,
,
,
is not positive. Thus, we have one of the following four cases.
There exists a first time
such that
There exists a first time
such that
There exists a first time
such that
There exists a first time
such that
In cases (1), we have
which is contradiction to
.
In cases (2), we have
which is contradiction to
.
In cases (3), we have
which is contradiction to
.
In cases (4), we have
which is contradiction to
.
Thus, the solutions ,
,
,
of Equation (Equation2
(2)
(2) ) are positive for all t>0.
Lemma 2.2
All feasible solutions
of Equation (Equation2
(2)
(2) ) are bounded by the region
(3)
(3)
Proof.
Parameters of Equation (Equation2(2)
(2) ) are assumed to be positive. Using the fact
, the system (Equation2
(2)
(2) ) reduced to the following system:
(4)
(4) From the third equation in (Equation3
(3)
(3) ) we can obtain
(5)
(5) It follows that
(6)
(6) where
is the initial value of total number of people. Thus,
(7)
(7) then
(8)
(8) Form the fourth equation in (Equation3
(3)
(3) ), we have
(9)
(9) it follows that
(10)
(10) where
represents the initial value of cumulative density media coverage. Thus,
(11)
(11) Hence, for the analysis of model (Equation2
(2)
(2) ), we get the region which is given by the set:
(12)
(12) which is a positively invariant set for (Equation2
(2)
(2) ). So we only need to consider dynamics of system on the set Ω.
3. Basic reproduction number and alcohol free equilibrium
It is easy to know that system (Equation2(2)
(2) ) always has an alcohol free equilibrium
. In the following, the basic reproduction number of system (Equation2
(2)
(2) ) will be obtained by the next generation matrix method formulated in [Citation7]. Let
, system (Equation2
(2)
(2) ) can be written as
where
The Jacobian matrices of
and
at the alcohol free equilibrium
are, respectively,
where
The basic reproduction number is
By the Theorem 2 of [Citation7], we have the following theorem.
Theorem 3.1
For system (Equation2(2)
(2) ), the alcohol free equilibrium
is locally asymptotically stable if
and unstable for
.
Theorem 3.2
If the alcohol free equilibrium
is globally asymptotically stable.
Proof.
Constructing Lyapunov function as follows:
(13)
(13) The derivative of
along with (Equation2
(2)
(2) ) is given by
(14)
(14)
Obviously, , and if
, we can get that
, hence, we obtain
. Furthermore,
holds if and only if
, A=0, Q=0. By LaSalle's invariance principle [Citation14], we can derive that the alcohol free equilibrium
is globally asymptotically stable. The proof is complete.
4. Existence and stability of alcohol present equilibrium
The alcohol present equilibrium of system (Equation2
(2)
(2) ) is determined by equations
(15)
(15) From the fourth equation in (Equation15
(15)
(15) ), we have
(16)
(16) From the third equation in (Equation15
(15)
(15) ), and together with (Equation16
(16)
(16) ), we can get
(17)
(17) From the second equation of (Equation15
(15)
(15) ), and together with (Equation16
(16)
(16) ), we can get
(18)
(18) Then substituting Equations (Equation16
(16)
(16) ) and (Equation18
(18)
(18) ) into the first equation in (Equation15
(15)
(15) ), we have a quadratic equation of
as follows:
(19)
(19) where
(20)
(20)
It is easy to see that is always positive and
is always negative if
. Hence, Equation (Equation19
(19)
(19) ) has an unique positive root
.
Through the above conclusion, the following result holds.
Theorem 4.1
When system (Equation2
(2)
(2) ) has an unique positive alcohol present equilibrium
.
Next we consider the stability of the alcohol present equilibrium.
Theorem 4.2
For the system (Equation2(2)
(2) ), the alcohol present equilibrium
is locally asymptotically stable if
.
Proof.
For the alcohol present equilibrium , the matrix at
can be written as
(21)
(21) where
(22)
(22) So the characteristic equation at the alcohol present equilibrium
is
(23)
(23) where
(24)
(24)
If , it's obvious that
,
,
are all negative and
,
,
,
,
are all positive. Here, we only calculate
,
(25)
(25) Similarly, we can easily get the coefficients
all are positive. According to the Routh–Hurwitz criterion [Citation12], we can get that the eigenvalues of the above characteristic equation have negative real parts under the conditions
(26)
(26) This shows that the alcohol present equilibrium
is locally asymptotically stable if
.
Conjecture
If , the alcohol present equilibrium
is globally asymptotically stable.
Remark
It's difficult to demonstrate the global stability, we only present simulation. How to prove the global stability of alcohol present equilibrium have still to be resolved.
5. Optimal control problem
5.1. The existence of optimal control
In this section, in order to achieve the goal of minimized alcoholics, we reconsider the model (Equation2(2)
(2) ) and formulate an optimal control problem (Equation28
(28)
(28) ) to reduce the numbers of alcoholics. Hence, we use two control variables
and
to reduce the number of alcoholics. The control
represents efforts intended to prevent the susceptible from contacting with the alcoholics, for example, the media coverage and educational programs are included. The control
represents the effort to encourage alcoholics to quit through proper treatment, such as taking pills or seeking other medical help. However, it is advisable to balance between the costs and the alcohol effects. In view of this, the optimal control problems to minimize the objective functional is given by
(27)
(27) subjects to
(28)
(28) with initial conditions
(29)
(29) Here,
is measurable and
, for all
, where
is the final time, the constant
are relative weight coefficients that balance the cost and the number of alcoholics. The optimal control problem (Equation28
(28)
(28) )
Next, we will study andanalyse the existence of the optimal control. To determine existence of optimal control, we use the result of Fleming and Rishel [Citation8].
Theorem 5.1
There exists an optimal control such that
(30)
(30) subjects to the control system (Equation28
(28)
(28) ) with initial conditions (Equation29
(29)
(29) ).
Proof.
To prove the existence of an optimal control, the following conditions must be satisfied.
The controls and corresponding state variables are non-negative values.
The control set
is convex and closed.
The right-hand side of the state system is bounded by a linear function in the state and control variables.
The integrand of the objective functional is concave on
.
There exist constants
and
such that the integrant
of the objective functional satisfies
(31)
(31) conditions 1, 2 and 4 are obvious satisfied. As for condition 3, we have already proved that four state variable are bounded, therefore,
(32)
(32) The last condition, we can easily see that there exist
, and for all
,
satisfies the equation
, which completes the existence of an optimal control.
5.2. The characterization of the optimal control
Since the existence of the optimal control is established. The optimal control are obtained from the necessary conditions which come from the Pontryagin's Maximum Principle in [Citation21]. To solve control problem actually, we should define Hamiltonian for the optimal control problems as follows:
(33)
(33)
where
,
, are the penalty multipliers satisfying
at optimal control
and
.
By applying Pontryagin's Maximum Principle we obtain the following theorem.
Theorem 5.2
There exist optimal control pairs and state variable of
that minimize the objective functional
over U. Given these optimal controls, there exist adjoint variable
satisfying the equations
(34)
(34) with the terminal conditions
(35)
(35) The optimal controls can be obtained from the following equations:
(36) (n34)
(36) (n34)
Proof.
By differentiating the Hamiltonian, we obtain the adjoint system:
(37)
(37) Therefore, the adjoint variable
satisfy Equation (Equation34
(34)
(34) ). By differentiating the Hamiltonian with respect to the controls, we have the following optimality conditions:
(38)
(38)
To determined an explicit expression for the optimal control without ,
. We consider the following three cases.
For
, we have
. Hence the optimal control is
(39)
(39)
For
, we have
. Hence
(40)
(40) Therefore,
implies that
(41)
(41)
For
, we have
. Hence,
(42)
(42) This implies that
(43)
(43) Combining these three cases, the optimal control is characterized as
(44)
(44) Similar to control variable
, we get
(45)
(45) This completes the proof.
6. Numerical simulations and discussions
In order to demonstrate the analytic results we have obtained, numerical experiments are carried out using Matlab. We choose a set of values of parameters in Table . We assume that the initial condition of (Equation2(2)
(2) ) is
. For the above set of parameter values, we obtain an unique positive alcohol present equilibrium
.
Table 1. The parameters description of the alcoholism model.
Now, we choose , it is easy to get
. From Figure , we know that the alcohol free equilibrium
is globally asymptotically stable. Then, we choose
, it can be calculated that
. From Figure (a, b), we see that the alcohol present equilibrium
is globally asymptotically stable.
In order to evaluate the effect of media coverage on the dynamics of alcoholics, we choose different values of m, γ (see Figures and ). In Figure , we compare the density of alcoholic individuals with different values of m (0, 0.1, 0.3, 0.9), the simulation shows that the more implementation rate of media coverage, the lower the final alcoholics' density. On the contrary, in Figure , the more the depletion rate of media coverage, the higher the final alcoholics' density. So we know that the number of alcoholics will reduce substantially with the media coverage increasing.
Now, we explore numerically an optimal control for the system (Equation28(28)
(28) ). To determine it, first, we choose
,
. Figure (a, b) indicate that the changes of population
and
with the different values of
,
. From Figure (a), we can easily get that the system with control is better than that without control. When
,
,
tends to a lowest value. When
is better than
, and comparing with
,
is better. From Figure (b), we know that the density of
decreases with the values of
and
increasing.
If we change weight coefficients into ,
, we know that there is little difference between Figures (a, b) and (a, b). So we can draw the same conclusion.
Figure shows the process of control measure with different weight coefficients. In Figure (a), when we choose ,
, the simulation shows that the control
starts from around 0.73, then it continues the slowly declining trend and finally toward a level of 0. But the control
is relatively stable approaching to 0. In Figure(b), we choose
,
, it's obvious that the control
reaches peak about 0.51 at the time about 10, then it gradually decrease to 0 with the increase of time. As for the control
, it gradually reduces from 0.2 to 0.
In this paper, we proposed an alcoholism model with media coverage. Unlike some of other previous model, we have taken into account the impact of media on alcoholics. First, we get the basic reproduction number by next generation matrix method. Second, we obtain that the alcohol free equilibrium of the model is globally asymptotically stable if
. For
, the alcohol present equilibrium is locally asymptotically stable. Next, we establish the optimality system and obtain the necessary condition of optimality by the Pontryagin's maximum principle. At last, we use numerical method to simulate outcomes which we have been proved.
Our results show that media coverage has substantial influence on the dynamics of alcoholics and can greatly influence the spread of the drinking, thus, it is crucial to educate people by media coverage to quit drinking.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- N.G. Becker and D.N. Starczak, Optimal vaccination strategies for a community of households, Math. Biosci. 139 (1997), pp. 117–132.
- M.L. Brandeau, G.S. Zaric, and A. Richter, Resource allocation for control of infectious diseases in multiple independent populations: Beyond cost-effectiveness analysis, J. Health Econ. 22 (2003), pp. 575–598.
- Y. Cai, Y. Kang, M. Banerjee, and W. Wang, A stochastic SIRS epidemic model with infectious force under intervention strategies, J. Differ. Equ. 259 (2015), pp. 7463–7502.
- Y. Cai, Y. Kang, M. Banerjee, and W. Wang, A stochastic epidemic model incorporating media coverage, Commun. Math. Sci. 14 (2016), pp. 892–910.
- J.F. Cao, W.T. Li, and F.Y. Yang, Dynamics of a nonlocal SIS epidemic model with free boundary, Discrete Contin. Dyn. Syst. Ser. B 22 (2017), pp. 247–266.
- C. Castilho, Optimal control of an epidemic through educational campaign, Electron. J. Diff. Equ. 2006 (2006), pp. 1–11.
- P. van den Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci. 180 (2002), pp. 29–48.
- W.H. Fleming and R.W. Rishel, Deterministic and Stochastic Optimal Control, Springer, Berlin, Germany, 1975.
- S. Ghosh, S. Guria, and M. Das, Alcohol as a risk factor for cancer burden: A review, Proc. Zoolog. Soc. 69 (2016), pp. 32–37.
- H.F. Huo and L.X. Feng, Global stability for an HIV/AIDS epidemic model with different latent stages and treatment, Appl. Math. Model. 37 (2013), pp. 1480–1489.
- H.F. Huo and X.M. Zhang, Complex dynamics in an alcoholism model with the impact of Twitter, Math. Biosci. 281 (2016), pp. 24–35.
- M. Kot, Elements of Mathematical Biology, Cambridge University Press, Cambridge, 2001.
- H. Laarabi, A. Abta, and K. Hattaf, Optimal control of a delayed SIRS epidemic model with vaccination and treatment, Acta Biotheor. 63 (2015), pp. 87–97.
- J.P. LaSalle, The Stability of Dynamical Systems, Regional Conference Series in Applied Mathematics, Society for Industrial and Applied Mathematics, Berlin, NJ, 1976.
- W. Lechner, A.M. Day, J. Metrik, A.M. Leventhal, and C.W. Kahler, Effects of alcohol-induced working memory decline on alcohol consumption and adverse consequences of use, Psychopharmacology 233 (2016), pp. 83–88.
- Y. Liu and J. Cui, The impact of media converage on the dynamics of infectious diseases, Int. J. Biomath. 1 (2008), pp. 65–74.
- R.S. Liu, J.H. Wu, and H.P. Zhu, Media/psychological impact on multiple outbreaks of emerging infectious disease, Comput. Math. Methods Med. 8 (2007), pp. 153–164.
- S.H. Ma, H.F. Huo and X.Y. Meng, Modelling alcoholism as a contagious disease: A mathematical model with awareness programs and time delay, Discrete Dynamics in Nature and Society, 2015, Article ID 260195, 13 pages.
- G. Mulone and B. Straughan, Modeling binge drinking, Int. J. Biomath. 5 (2012), pp. 57–70.
- S. Mushayabasa and C.P. Bhunu, Modelling the effects of heavy alcohol consumption on the transmission dynamics of gonorrhea, Nonlinear Dynam. 66 (2011), pp. 695–706.
- L.S. Pontryagin, V.G. Boltyanskii, R.V. Gamkrelidze, and E.F. Mishchenko, The Mathematical Theory of Optimal Processes, Interscience, New York, USA, 1965.
- J. Rehn, C. Mathers, S. Popova, M. Thavorncharoensap, Y. Teerawattananon, and J. Patra, Global burden of disease and injury and economic cost attributable to alcohol use and alcohol-use disorders, Lancet 373 (2009), pp. 2223–2233.
- S. Samanta, S. Rana, A. Sharma, A.K. Misra, and J. Chattopadhyay, Effect of awareness programs by media on the epidemic outbreaks: A mathematical model, Appl. Math. Comput. 219 (2013), pp. 6965–6977.
- F. Sanchez, X. Wang, C. Castillo-Chavez, D.M. Gorman, and P.J. Gruenewald, Drinking as an epidemica simple mathematical model with recovery and relapse, in Therapists Guide to Evidence-Based Relapse Prevention: Practical Resources for the Mental Health Professional, K.A. Witkiewitz and G.A. Marlatt, eds., Academic Press, Burlington, Canada, 2007, pp. 353–368.
- O. Sharomi and A.B. Gumel, Curtailing smoking dynamics: A mathematical modeling approach, Appl. Math. Comput. 195 (2008), pp. 475–499.
- R. Straus and S.D. Bacon, Drinking in College, Greenwood Press, New Haven, 1971.
- Y. Wang, J.D. Cao, Z. Jin, H.F. Zhang, and G.Q. Sun, Impact of media coverage on epidemic spreading in complex networks, Statist. Mech. Appl. 392 (2013), pp. 5824–5835.
- X.Y. Wang, K. Hattaf, H.F. Huo, and H. Xiang, Stability analysis of a delayed social epidemics model with general contact rate and its optimal control, J. Ind. Manag. Optim. 12 (2016), pp. 1267–1285.
- H. Xiang, L.X. Feng, and H.F. Huo, Stability of the virus dynamics model with Beddington-DeAngelis functional response and delays, Appl. Math. Model. 37 (2013), pp. 5414–5423.
- H. Xiang, Y.L. Tang, and H.F. Huo, A viral model with intracellular delay and humoral immunity, Bull. Malays. Math. Sci. Soc. (2016). doi: 10.1007/s40840-016-0326-2.
- G. Zaman, Qualitative behavior of giving up smoking models, Bull. Malays. Math. Sci. Soc. 34 (2009), pp. 403–415.
- G. Zaman, Optimal campaign in the smoking dynamics, Comput. Math. Methods Med. (2011), Article ID 163834, 9 pages.