Liquor Habit Transmission Model

In this paper, we analyse population dynamics of liquor habit. Liquor free and Liquor endemic equilibrium are computed. The local and global stabilities of the proposed problem are established. The numerical simulation is given to validate the transmission of population indifferent compartment using state-space model.


Introduction
The Ayurvedas depicted that the alcohol behave as a medicine if it is taken for the purpose of meditation otherwise it behave as a poison if it is taken in addiction manner for the humans.The liquor habits reflect person's social status and cultural prestige.It is observed that transmission of liquoring spreads frequently these days due to availability of the liquor in the market.There is a similarity between spread of infectious disease and liquor habits.In other words, liquor habit can be treated as a virus which transmits among the compartments by social pressure like parties with friends, peers and executive meetings.
In this paper, we analyze quantitative model of a liquor habit transmission in a population similar to SEIRmodel.The notations are described in Section 2. The mathematical model and basic reproduction number are formulated in Section 3. The local and global stability are derived in Sections 3.2.1 and 3.2.2respectively.Numerical simulations are illustrated in Section 4 using State-Space model.Discussions and conclusions are given in Section 5.

Notations
The model is derived using following notations.

( ) ( ) ( ) ( )
: the portion of the susceptible who has never liquored; : the portion of the total population who liquor less than 2 pegs day; : the portion of the total population who liquor more than 2 pegs day; : the portion of the total population

Mathematical Model
Anybody in the population is susceptible to liquor.In general, one starts taking pegs and increases it gradually.The consumption of two pegs is considered as a normal which falls in E-compartment.The work stress, financial stress, and many more factors to go for more than two pegs which fall under I-compartment.Certain fraction of population from E-and I-compartments may be removed.It is taken as R-compartment.Thus, it resembles to-SEIR-model.
Let us call , , S E I and R the portion of total population of each class.Now, we make SEIR-model with some assumptions.The portion of liquored person's who increases liquor habit at a rate proportional to ( ) 1 and 0, 0, 0, 0 Since, an epidemic model occurs in a short time period, we ignore moving portion of removal from liquor.So, we will analyze the first three equations forming new reduced system Adding all these three equations, we have ( ) Therefore, the feasible region for ( 2) is Now, the basic reproduction number 0 R will be found by using the next generation matrix.It is easy to see that (2) always has a liquor free equilibrium , where dash denotes derivative.So that where ( ) denotes the rate of appearance of new liquored in compartment and ( ) X υ represents the rate of transfer of liquors, which is given as F and V are 2 2 × matrices defined as ( ) ( ) Hence, the basic reproduction number 0 Next, we need to discuss equilibrium of the liquor habit system.

Equilibrium
The liquor free equilibrium is locally asymptotically stable if all the eigenvalues of the matrix have positive real values (Al-Amoudi et al. (2014) [1]).Theorem 1 (Johnson (2004) [2]): Consider the liquor transmission model given by ( 2) with X ′ .If 0 X is a liquor free equilibrium of the model, then 0 X is locally asymptotically stable if 0 1 R > , and unstable if 0 1 R < , where 0 R is given by (4).Proof Let j F V = − .Since V is a non-singular matrix and F is non-negative, Thus j F V − = − + has the Z-sign Then, we have j − is a non-singular matrix ⇔ ( ) FV − is non-negative, all eigenvalues of 1 FV − have magnitude less than or equal to 0 R .Thus, ( ) Hence, ( ) 0 Similarly, it follows that ( ) Hence, ( ) 0 It follows that ( ) 0 S j > if and only if 0 1 R > .

Local Stability
The liquor free equilibrium is stable if all the eigenvalues of the Jacobian matrix of the system (1) have negative real parts.For this, the Jacobian of the system (1) at

Global Stability
The liquor free equilibrium is globally stable if ( ) .

Numerical Simulation
In this section, we perform numerical simulation of the system (1) using with the state-space model.

State-Space Model
State equation ( ) where u-input, y-output, x-state vector.
The general state-space description for a linear time invariant, continuous time dynamical system is

Discussion and
From Figure 2(a), it is observe that 11% of the population from the susceptible compartment start liquoring exponentially in first five weeks.Liquoring habit in the susceptible class is small initially but thereafter it increases exponentially.that more than two pegs are taken immediately after 3 rd week and increases exponentially.Almost 15% gets chain liquor (those who can't survive without liquor).The removal compartment (Figure 2(d)), 44% enters either into less than two pegs i.e.E-compartment or into never liquor that is S-compartment.
In this paper, a nonlinear mathematical model for liquor transmission is analysed.The local and global stability of the liquor free equilibrium point are established.It is proved that the free equilibrium is locally asymptotically stable when basic reproduction number 0 1 R < and global stability of liquor transform 1 1 0 FV − − > .This model can be extended by educating youth for non-liquoring via advertisement, rehabilitation centre etc.
the portion of liquored will decrease with the some rate.Here, β is called the effective liquor- ing rate.The portion of liquor person's starting liquoring habit with a rate R ρ (progression rate) The portion of excessive liquored with the rate c with 0 c > .The portion of removing liquor habit with the rate I γ with 0 γ > and E λ with 0 λ > .See the transfer diagram given in Figure 1.From Figure1, the model is formulated as following system of differential equations.

.
Using this state-space model, we find the solution of the system(1).Taking data as follows: We carry out the simulation.The results are shown in Figure2.