American Journal of Applied Mathematics
Volume 5, Issue 1, February 2017, Pages: 31-38

Smoking as Epidemic: Modeling and Simulation Study

Sintayehu Agegnehu Matintu

School of Mathematical and Statistical Sciences, Hawassa University, Hawassa, Ethiopia

Sintayehu Agegnehu Matintu. Smoking as Epidemic: Modeling and Simulation Study. American Journal of Applied Mathematics. Vol. 5, No. 1, 2017, pp. 31-38. doi: 10.11648/j.ajam.20170501.14

Received: December 4, 2016; Accepted: December 23, 2016; Published: February 23, 2017

Abstract: In this paper we propose smoking epidemic model which analyzes the spread of smoking in a population. The model consists of five compartments corresponding to five population classes, namely, potential-moderate-heavy-temporarily recovered- permanently recovered class. The basic reproduction number R0 has been derived, and then the dynamical behaviors of both smoking free equilibrium and smoking persistent equilibrium are analyzed by the theory of differential equation, and Numerical simulation has been carried out and the results have confirmed the verification of analytical results. Sensitivity analysis of R0 identifies β1, the transmission coefficient from potential smokers to moderate smokers and β2, the transmission coefficient from potential smokers to heavy smokers, as the most useful parameters to target for the reduction of R0.

Keywords: Smoking Model, Reproduction Number, Equilibrium Value, Stability, Sensitivity Analysis, Numerical Simulation

1. Introduction

Smoking is a practice in which a substance is burned and the resulting smoke breathed in to be tasted and absorbed into the bloodstream. Most commonly the substance is the dried leaves of the tobacco plant which have been rolled into a small square of rice paper to create a small, round cylinder called a "cigarette" [1]. Epidemiology is concerned with the distribution of diseases in populations and the factors that influence the transmission of diseases [2]. One of an interesting area in epidemiology to study is smoking subject. There are a lot of studies that has been done on the epidemics of smoking [3-12]. In order to understand the dynamics of this disease, like many infectious diseases, mathematical models can be used. Most of these models are compartmental models, with the population being divided into compartments with the assumptions about the nature and time rate of transfer from one compartment to another. The basic concepts in these models are that all people in a community start as healthy. Healthy people may become infected by diseases. But infected people may become healthy again in community.

Castillo-Garsow et al. [3] presented a mathematical model for giving up smoking. In their model the total population is classified into three compartments namely; potential smokers, smokers and quitters (former smokers). Later on, Sharomi and Gumel [4] extended the basic model to account for variability in smoking frequency, by introducing two classes of smokers who temporarily and permanently quit smoking. Zaman [5] derived and analyzed a smoking model taking into account the occasional smokers compartment in the basic model. In [6], Z. Alkhudhari, S. Al-Sheikh and S. Al-Tuwairqi adopted the model developed and studied in [6] and considered the effect of peer pressure on temporary quitters. In 2014 [7], Z. Alkhudhari, S. Al-Sheikh and S. Al-Tuwairqi introduced a new model by dividing the smokers into two subclasses: occasional smokers and heavy smokers, and the impact of these two subclasses on the existence and stability of equilibrium points. They also studied the effect of occasional smokers on potential smokers, and later in 2015 [8] they studied the effect of heavy smokers on potential smokers.

In this paper, we develop a new model based on the idea given in [6, 7 and 8] by considering, basically, the effect of both moderate and heavy smokers on the potential smokers, and the effect of heavy smokers on temporarily quit smokers. Furthermore, it assesses the impact of this peer pressure on the existence and stability of steady state solutions, and the result is supported with numerical simulation.

2. Mathematical Model Formulation

Let the total population size at time t is denoted by . We divide the total population  into five subclasses. These are described as follows: potential smokers (non-smokers)  but likely to get infected (start smoking) in future, moderate (light) smokers , heavy (addicted) smokers , smokers who temporary quit smoking  and smokers who permanently quit smoking  such that .

2.1. Flow of People Among the Compartments

People will join the potential smoker’s compartment  at a constant recruitment rate of. This recruitment is due to the total natural birth rate . Some of these people will vacate this compartment at a constant death rate of , due to the total natural death rate . It is considered that the potential smokers start smoking due to peer influence of moderate and heavy smokers, will vacate the potential compartment with the rates  and . Here,  is the contact rate between potential and moderate smokers and  is the contact rate between potential and heavy smokers.

The moderate smokers, are increased when potential smokers start to smoke, with rates  and . Some other people will leave the compartment with the rates  and , in which,  is the rate at which moderate smokers becomes heavy smokers without peer influence and  is the contact rate between moderate and heavy smokers.

People will join into heavy smoker’s class  with the rates  and . Some others will leave with the rates and. Here,  and  are the contact rate between heavy smokes and those who temporarily quit smoking and the rate of quitting smoking, respectively.

People will enter into temporary quitter’s compartment  with a rate  and leave with the rates  and , where,  is the fraction of heavy smokers who temporarily quit smoking (at a rate ).

The permanent quitters,  are increased when heavy smokers begin to smoke, with the rate  and some people of this compartment will die with the rate , where,  is the remaining fraction of heavy smokers who permanently quit smoking (at a rate).

2.2. Flow Diagram of Smoking Epidemic Model

Here, we have given the flow diagram of the model. The compartments of the model are represented by rectangular colored boxes. The flow directions of the people among the compartments are represented by directed arrows.

Figure 1. Flow diagram of Smoking Epidemic Model.

2.3. Model Assumptions

Here, for the above flow diagram of smoking endemic model, the following assumptions are made.

The total population is divided into two diseases (smoker) and three non diseases (non smoker) compartments.

The recruitment rate  is different from the natural death rate  so that the total population is not constant.

The potential smokers start to smoke under the peer influence of both moderate and heavy smokers.

All the newly born people are assumed to join only potential smoker’s compartment.

The people of all the compartments will die naturally with the same rate of .

No other influence to smoke other than smokers.

Temporarily quitted smokers could only be allowed to relapse into heavy smokers’ class.

Legal age could be considered to start smoking, in the model.

No person is more popular than another, so that interactions between any two compartments are equally likely.

2.4. Model Equations

By adding the rates at which people enter the compartment and subtracting the rates at which people leave the compartment, we get an equation for the rate at which the population of each compartment changes over time. Thus, by considering the above assumptions, smoking epidemic model presented in figure  may be represented by the following system of first order nonlinear ordinary differential equations:

(1)

With the initial densities

(2)

By introducing  and , equation can be modified (scaled) as,

(3a)

(3b)

(3c)

(3d)

(3e)

Provided . The initial densities in equation  may also be modified as follows:

(4)

Since the variable  of equation  does not appear in the first four equations, we shall consider the equations  and we shall find the feasible (solution) region of equations so that the model makes biological sense, as follows:

(5)

Therefore, as it has been shown in [14], the solution of equation  can be proved that , where .

Thus, the feasible solution of human population for the models  is restricted to the region,

We denote the boundary and the interior of  by  and , respectively.

3. Analysis of the Model

In this section we are going to analyze the smoking endemic model in equations .

3.1. Existence of Smoking Free Equilibrium

Many epidemiological models have a disease free equilibrium (DFE) at which the population remains in the absence of disease [10]. In this model, the smoking free equilibrium  is a steady state solution in which the whole population is non-smokers. Mathematically, the smoking free equilibrium of the differential equations given in equation  is obtained by setting  provided that. As a result, the smoking free equilibrium point is

3.2. The Basic Reproduction Number

Like as it has been described in [9], analyzing the differential equations given in equation can also give insight into how smoking spreads and how the spread can be limited. An important tool for analyzing a model of this type is its reproductive number.

The basic reproduction number, denoted, is defined as the number of new infections produced by a typical infective individual in a susceptible population at a disease free equilibrium [10]. In model , taking the infected compartments to be moderate and heavy smokers, the basic reproduction number is then the sum of the transmission rate from potential smokers to moderate smokers divided by the sum of the natural death rate of the moderate smokers population and the proportion of individuals who enter into heavy smoker’s compartment, i.e.,

(6)

However, the moderate smokers  in equation can be computed from equation  as follows:

(7)

By plugging equation  into equation, the basic reproduction number at the disease free equilibrium  is given as,

(8)

As it has been briefly discussed in [10], the dimensionless constant  tells us that if , then the DFE in , is locally asymptotically stable, and the disease cannot invade the population, but if , then the endemic equilibrium exists in  and is the DFE is unstable and invasion is always possible.

3.3. Sensitivity Analysis of

It is critical to take various actions to control the system parameters so that  is remarkably below one. This section is intended to propose some effective measures for achieving this goal. To examine the sensitivity of  to each of its parameters, following Arriola and Hyman [13], the normalized forward sensitivity index with respect to each of the parameters is calculated as follows:

It can be seen that, among these six parameters,  is the most sensitive to the change in ,  and . An increase or decrease in  will cause an increase or decrease in  with the same proportion. Physically, however, decreasing the natural birth rate  is neither practical nor ethical. An increase (or decrease) in the value of  or  leads to a corresponding increase (or decrease) in . As opposed to this, the other two parameters  and  have an inversely proportional relationship with , an increase in  and  will bring about a decrease in  with a proportionally smaller size of decrease. Recall that  is the natural death rate of the population. It is clear that increase in  is neither ethical nor practical. Additionally, the parameter  may have directly or inversely proportional relationship with the reproduction number. Thus we choose to focus on one of three parameters: either , the transmission rate from potential smokers to moderate smokers or , the transmission rate from potential smokers to heavy smokers or , the rate of quitting smoking. Given ’s sensitivity to  and, it seems sensible to focus efforts on the reduction of  and . In other words, this sensitivity analysis tells us that prevention is better than quitting smoking.

3.4. Local Stability of Smoking Free Equilibrium

The local stability of  is then determined from the signs of the eigenvalues of the Jacobian matrix. At the smoking free equilibrium , the Jacobian matrix is given by:

(9)

The characteristic equation of this matrix is given by det , where  is a square identity matrix of order. The equation thus becomes

(10)

Where,

The roots of the characteristic equation are the eigenvalues of the Jacobian matrix. It is clear that the one eigenvalue is determined from , and is . All the other remaining eigenvalues are determined from the quadratic equation

(11)

Now by using Routh-Hurwitz criterion [15], it can be stated that the smoking free equilibrium  is locally asymptotically stable if and only if the following inequalities satisfied: and

Theorem: If, then  and  are greater than zero.

Proof: Since each parameter’s value is greater than zero, it is trivial that  if. Because . Here, if  then . This implies that  Hence, all roots of equation  have negative real parts if .

So all eigenvalues of the Jacobian matrix are negative if , and hence,  is locally asymptotically stable. If , then equation  may be modified as  which has the rootsand, and hence,  is locally stable. If , then  and it is clear that the quadratic equation has some positive roots. This leads to conclude that the smoking free equilibrium becomes unstable.

3.5. Existence and Local Stability of Endemic Equilibrium

To find the smoking present equilibrium of the system of equation setting  provided that at least one of the infected compartments’ is non zero. Thus, the model has smoking present equilibrium , where

,

,

,

,

and  satisfies the equation . Analytically, however, it is tricky to obtain  as in terms of parameters, and hence the existence of smoking present equilibrium will be verified with numerical simulation.

For the smoking present equilibrium , the Jacobian matrix is given as,

(12)

Where,

and .

The characteristic polynomial for  is given by det , where  is a square identity matrix of order. The polynomial becomes

(13)

Where,

,

,

,

.

Algebraically, the epidemic equilibrium value is asymptotically stable by Routh-Hurwitz criteria [15], if  and . These conditions ensure that all of the four eigenvalues found from  have negative real parts. However, the stability of the endemic equilibrium  will be illustrated using numerical simulations, by imposing the value of parameters and initial population of compartments.

4. Numerical Simulations

A number of numerical simulations have been carried out using MATLAB to illustrate the dynamics of the system. We use the following initial population such that:

It can be concluded that the model given in this paper can represent the dynamics of Tobacco use in real life if it is possible to find the reasonable values for parameters. The values of some parameters have been studied in some literatures while others don’t. Here, we consider the idea presented in [11, 12]. When a person first becomes a smoker it is not likely that she/he quits for several years since tobacco contains nicotine, which is shown to be an addictive drug. We assume  to be the average value of  and  years, that is,  years. Hence, the value of  is set to  days. The natural death rate  is per  per year (currently  in the US). There is strong evidence that the attitude towards smoking is starting from high (junior) school time. Therefore is set to be  days.

For the following numerical simulations we use the following parameters:
for , and  and the others are as it is for.

Figure 2 show that every compartment approaches its smoking free equilibrium value  as time increases. The results show that the smoking free equilibrium value of the system is asymptotically stable, which is as expected from our theoretical analysis as . Additionally, the stability of smoking free equilibrium value is not affected by the initial values of each compartment as it is presented in figure 2. On the other hand, we have considered another case in which  as it has been shown in figure 3. The existence and stability of smoking present equilibrium value(s) has been simulated in figure 3 even though it is messy to obtain analytically. Here, the stability of smoking present equilibrium value(s) is/are independent of the initial values. Furthermore, the parameters  and  have no significant on the stability of both smoking free and smoking present equilibrium values as long as their stability depend on the reproduction number .

Figure 2. Time series plots of Potential smokers, Moderate smokers, Heavy smokers, Temporary quitter smokers and Permanent quitter smokers with different initial densities for .

Figure 3. Time series plots of Potential, Moderate, Heavy, Temporary and Permanent quitter smokers with different initial densities for .

5. Conclusions

In this paper, we have proposed a mathematical model describing the dynamics of smoking. This model consists of five differential equations describing the rate of change of the subpopulation  and . The model has smoking free equilibrium and smoking present equilibrium. The threshold  is derived and it is shown that if , then the smoking will disappear, while if , then the smoking persists with the population which means that after some period of time the smoking will become hazardous. Numerical simulation has been carried out and the results have confirmed the verification of analytical results. Sensitivity analysis of  identifies , the transmission coefficient from potential smokers to moderate smokers and , the transmission coefficient from potential smokers to heavy smokers, as the most useful parameters to target for the reduction of . In future work, the possible extensions are to consider that Smoking related death will be introduced in the model and the temporary quitters decide to quit smoking permanently.

References

1. http://en.wikipedia.org/wiki/Smoking (3th October, 2016)
2. https://en.wikipedia.org/wiki/Epidemiology (17th November, 2016)
3. C. Castillo-Garsow, G. Jordan-Salivia, and A. Rodriguez Herrera, Mathematical models for the dynamics of tobacco use, recovery, and relapse, Technical Report Series BU-1505-M, Cornell University, Ithaca, NY, USA, 1997.
4. O. Sharomi, A.B. Gumel, Curtailing smoking dynamics: A mathematical modeling approach, Applied Mathematics and Computation 195 (2008) 475–499.
5. G. Zaman, Qualitative behavior of giving up smoking model; Bulletin of the Malaysian Mathematical Sciences Society (2) (2011) 403-415.
6. Z. Alkhudhari, S. Al-Sheikh, S. Al-Tuwairqi, Global Dynamics of a Mathematical Model on Smoking, ISRN Applied Mathematics vol. (2014), Article ID 847075.
7. Z. Alkhudhari, S. Al-Sheikh, S. Al-Tuwairqi, The effect of occasional smokers on the dynamics of a smoking model, International Mathematical Forum, 9, No. 25 (2014), 1207- 1222.
8. Z. Alkhudhari, S. Al-Sheikh, S. Al-Tuwairqi, The effect of heavy smokers on the dynamics of a smoking model,International Journal of Differential Equations and Applications, 14, No. 4 (2015),1311-2872.
9. B. Benedict, Modeling Alcoholism as a Contagious Disease: How "Infected" Drinking Buddies Spread Problem Drinking. SIAM News, Vol. 40, No. 3, 2007.
10. P. van den Driessche, J. Watmough, Reproduction numbers and sub-thershold endemic equalibria for compartmental models of disease transmission, Math. Biosci. 180 (2002) 29–48.
11. G. Zaman, Y. H. Kang, and I. H. Jung, "Stability analysis and optimal vaccination of an SIR epidemic model," BioSystems, vol. 93, no. 3, pp. 240–249, 2008.
12. American National Institute of Drug Abuse, "Cigarettes and Other Nicotine Products," http://www.nida.nih.gov/pdf/ infofacts/Nicotine04.pdf.
13. L. Arriola, J. Hyman, Lecture notes, Forward and adjoint sensitivity analysis: with applications in Dynamical Systems, Linear Algebra and Optimization, Mathematical and Theoretical Biology Institute, Summer, 2005.
14. Z. Alkhudhari, S. Al-Sheikh, S. Al-Tuwairqi, Global Dynamics of a Mathematical Model on Smoking, ISRN Applied Mathematics vol. (2014), Article ID 847075.
15. L. Edelstein Keshet, Mathematical Models in Biology. Random House, New York, 1988. SIAM edition.

 Contents 1. 2. 2.1. 2.2. 2.3. 2.4. 3. 3.1. 3.2. 3.3. 3.4. 3.5. 4. 5.
Article Tools