In this paper, we propose and analyze a fuzzy model for the health risk challenges associated with alcoholism. The fuzziness gets into the system by assuming uncertainty condition in the measure of influence of the risky individual and the additional death rate. Specifically, the fuzzy numbers are defined functions of the degree of peer influence of a susceptible individual into drinking behavior. The fuzzy basic risk reproduction number is computed by means of Next-Generation Matrix and analyzed. The analysis of reveals that health risk associated with alcoholism can be effectively controlled by raising the resistance of susceptible individuals and consequently reducing their chances of initiation of drinking behavior. When perceived respectable individuals in the communities are involved in health education campaign, the public awareness about prevailing risks increases rapidly. Consequently, a large population proportion will gain protection from initiation of drinks which would accelerate their health condition into more risky states. In a situation where peer influence is low, the health risks are likely to be reduced by natural factors that provide virtual protection from alcoholism. However, when the perceived most influential people in the community engage in alcoholism behavior, it implies an increase in the force of influence, and as such, the system will be endemic.

1. Introduction

Alcohol consumption that leads to health-related epidemics has been one of the leading causes of mortality of individuals worldwide. In Tanzania, just like in other developing countries, the situation gets worse as alcohol consumption becomes an important part of community social activities. However, available literature cite alcohol consumption behavior is one among the risk factors for various health challenges in Tanzanian population [19]. While communities identify themselves with alcohol consumption, a large part of the religious population forbid its uptake. Being religious therefore plays an important role in promoting health and molding an individual’s behaviors [10, 11]. Increased health defect cases such as malnutrition, chronic pancreatitis diseases, liver cirrhosis, and some types of cancer have been found to be associated to alcoholic behavior [1219]. While in some parts of the world individuals believed that a low level of alcohol consumption prevents thrombosis [20], recent studies challenge the belief. For example, Griswold et al. [13] established that any amount of alcohol consumption leads to health risks.

Mathematical modeling has been very useful in providing solutions to epidemiology-related problems for both autonomous and nonautonomous system models [12, 14, 15, 2129]. While autonomous systems literally mean the use of constant controls, a nonautonomous system occurs when a continuous variable control is used in the system as a function of time [22]. However, an application of fuzzy set theory in solving dynamical systems has recently been an interesting research area [3033].

In this paper, we use an application of fuzzy set theory approach to modify the basic SPLMAR model developed in system (1) into a fuzzy model. Two main distinguishing aspects of the model are as follows: social cultural beliefs as an integral part of the society [3436] and the staged process in which alcoholism behaviors take in the spread of health risks [12, 14, 15, 24, 26, 29] are considered.

2. Fuzzy Model Formulation

This paper presents crisp model (1) with six time-dependent state variables described as follows: susceptible population, —comprising of individuals susceptible to alcohol drinking habits; protected population, —comprising of individuals with strong religious beliefs that provide them with virtual protection from alcohol drinking; and low-risk population, —individuals at low risk who drink alcoholic beverages once in a while. Others include medium/moderate-risk population, —individuals at medium/moderate risk who drink alcoholic beverages regularly; alcoholics, —high-risk drinkers or alcohol addicts who have developed physical or psychological dependence on alcoholic beverages; and recovered population, —former drinkers who have voluntarily quit alcohol drinking due to various reasons. The model considers multirisk levels established under the drinking cultures of Tanzanian population with active religious beliefs by using constant model parameters defined in Table 1with , , , , , , and .

The model is developed on the fact that a nonalcoholic drinker acquires alcohol drinking habits through regular social contacts [12, 37] at a rate of defined as

We propose a fuzzy model using as the variable describing the degree of peer influence of a susceptible individual to initiate drinking behavior. The model development is made under the assumption that both and are the functions of . The physical meaning of parameters used in the model determined the choice of the two parameters. While is attributed to the spread of health risks associated with alcoholism in the community, translates the consequences of health risks by means of additional death rate. The fuzzy numbers and represent the likelihood that a susceptible individual will drink alcohol after prolonged contact with drinking individuals and the additional death rate induced by alcoholism, respectively. A fuzzy number may be defined as a generalization of a real number referring to a connected set of possible values, rather than referring to one single value, where each possible value is weighted between 0 and 1. Thus, the fuzzy system (3) of differential equations is established with nonnegative initial conditions of the state variables: , , , , and where

It is further assumed that the population is homogeneous and, at the initial level, drinking behavior is acquired by choice. Further, it is assumed that individuals in the recovered population do not develop a permanent immunity system against alcohol drinking. Similarly, individuals in the protected compartment acquire a nonpermanent virtual protection from alcohol drinking for their entire life in the compartment. On the other hand, a susceptible drinker acquires alcohol drinking habits through prolonged social contacts with drinkers in shared drinking venues. The force of peer influence is redefined into

3. Analysis of the Fuzzy System

Following Verma et al. [38] and Nandi et al. [32], we study the fuzzy SPLMAR model (3) and provide the model analysis and interpretations. The parameters and can be described through two fuzzy membership functions. For example, considering a situation where a drinking individual interacts with a susceptible member, the minimum amount of the degree of peer influence, , is required to have an impact on a susceptible member. That is to say, when , the impact of behavioral influence is considered negligible. The quantity is taken as a parameter whose exact value would depend upon both the attitude, public opinions towards the drinking behavior or the drinking individual, and willingness of a susceptible individual to conform with the peer pressures. As increases, the behavior inducement rate increases to the maximum which is equal to unity at . Further, it is also assumed that the degree of peer influence is bounded above where marks an upper bound. Therefore, the values of with an effect to the system lies in the interval of . The fuzzy membership function for the fuzzy number is given by

Similarly, we also assume that the addition death rate, , is a fuzzy number as it occurs due to the risk “transmission.” When , we have a negligible amount of risk transmissions implying that the additional death rate, , gets to the minimum level, say . As increases, the additional death rate increases and gets the highest value at . Since the additional death rate may not reach as its highest score due to several limitations, we let the maximum value for some real number such that . We therefore establish the fuzzy membership function of as follows:

The graphs of membership functions and are presented in Figures 1(a) and 1(b), respectively. For the model to be more realistic, we assume that the degree of peer influence of the studied group differs from among drinking individuals depending on their social influence in the community. Consequently, can be considered a linguistic variable with varying classifications. If we use and as, respectively, the central value and dispersion of each one of the fuzzy sets assumed by , we model each classification by using a triangular fuzzy number whose membership function is given in (8) and the graph of membership function is presented in Figure 2.

Let , , and be the family solutions of the fuzzy model system (3). These are the numbers of risky population proportions created as the result of social interactions between the susceptible members and risky individuals with social influence at time . Now, , , and are fuzzy numbers which lie in the interval .

3.1. Existence of Equilibrium Points

The equilibrium points of the fuzzy model (3) can be obtained by setting the RHS to be equal to zero and solve the system simultaneously. Thus, we have where , , , , , and are simplifying parameters.

The risk-free equilibrium point is obtained when the entire population is free from health risks associated with alcoholism, that is, implying that . Solving system (9) by maintaining risky classes at zero, we have with the positive constant .

Similarly, when , , , and consequently , the risk endemic equilibrium REE point (11) is established. where by , and are defined in and , , , and are the simplifying factors. Thus, the Theorem 1 is established.

Theorem 1. System (3) has a risk-free equilibrium and a unique risk endemic equilibrium .

3.2. Bifurcation and Fuzzy Basic Risk Reproduction Number

In this section, we first compute the basic health risk reproduction number, denoted as by using the Next-Generation Matrix method [39, 40]. The value of serves as the determinant to indicate whether the health risk epidemic is possible or not. Based on Diekmann et al.’s study [39], the risk transmission model consists of nonnegative initial conditions together with the following system of equations: where , , , and denotes transpose.

We then formulate the risk “transmissions” matrix and risk “transitions” matrix such that

Thus, by direct computation, we have

Now, we have the following Next-Generation Matrix

The basic risk reproduction number, , is given by the dominant eigenvalue of matrix Therefore,

With appropriate choices of , Theorem 2 is established

Theorem 2. The risk-free equilibrium is locally asymptotic stable when and unstable otherwise.

Proof. Consider the Jacobian matrix evaluated at the risk free equilibrium point bellow where , , , , , , , , , , , , , , , , , and
Using a trace determinant method, the risk-free equilibrium point is locally stable provided that and The trace and determinant of are, respectively, given by Since , it follows that is asymptotically stable whenever and unstable when provided that

The value of may vary significantly depending on different risk dynamics studied in the population or different populations involved in similar studies [33]. Similarly, the stability of the risk-free equilibrium also changes from unstable to stable when increases through 1. The system acquires a bifurcation at the risk-free equilibrium when . Suppose that the bifurcation value occurs at (see Figure 3) where is given as provided that and

The basic risk reproduction number presented in (21) is the function of , that is , but it can not be referred to as the fuzzy number since it is open to values exceeding unity. It is clear that is the function of the degree of social influence in the spread of the behavior. However, both and incline to their maximum level whenever , and as such, we have

We introduce a positive number such that , with an appropriate choices of ; we have a fuzzy set whose fuzzy expected value, , can be well defined. Therefore, the fuzzy basic risk reproduction number , which can be defined as the average number of secondary risk cases caused by one infected node introduced into entirely susceptible nodes [31], is given by

We know from the definition of fuzzy expected value that where is a fuzzy measure. For the purpose of this study, the possibility measure is given by

From , for a monotonic function , with set given as an interval , we let be the solution of the following equation:

Thus, where and .

Now, can be determined by considering the linguistic variable in three classes: “weak ,” “medium ,” and “strong .” Each of these classification is a fuzzy number based on , , and as shown in Figure 4. The classification of the degree of social influences in the community can be explained in three different cases as follows:

Case 1. In this case, we consider the weak degree of social influence where we define that . Suppose that there exists such that ; if c, we have We establish that . This can be translated that the fuzzy basic risk reproduction number , and hence, the extinction of the health risk is associated with alcoholism.

Case 2. In this case, we consider the medium degree of social influence where we define that and giving.

For , is the continuous function giving and . This translates as the fixed point and . Since is a continuous monotonic (increasing) function, the Intermediate Value Theorem suggests that there exists with in which the values of and coincide such that . Further, the average number of fuzzy basic risk reproduction number is higher than the number of secondary risk cases due to the medium level of social influence implying that the health risk associated with alcoholism is endemic.

Case 3. In this case, we consider the strong degree of social influence by defining and .
The results in (35) can be established.

For any given , is a monotonically decreasing and continuous function with and .

Therefore, is established as a fixed point such that

By dividing by throughout, we have the following results:

Now, since , it translates into endemic health risks associated with alcoholism.

3.3. Risk Control in Fuzzy Epidemic System

In this section, we perform the control analysis of the risk estimation in the population by using the fuzzy basic risk threshold . The spread of health risk in the proposed fuzzy SPLMAR model (3) depends on the degree of the social influence as an input value of the transmission factor . The description of existence and stability of the risk in the system is case-wise presented hereunder. Since the proposed fuzzy system (3) represents a family of systems depending on the parameter , these family systems can be simplified by a unique system of equations with the same results. It is shown that the bifurcation occurs at , that is, . (1)Weak influence: in this case, we have giving suggesting the extinction of the health risks associated with alcoholism in the community.(2)Medium influence: in this case, three possibilities may arise as follows: (i)If , then suggesting the risk free community(ii)If , then an indication of risk bifurcation(iii)If , then implying the risk endemic in the community(3)Strong influence: in this case, we have giving

The spread of the health risks depends upon the parameter . Let be an improved value of ; we can establish three possibilities in which the spread of health risks takes as follows: (i)If , then suggesting the health risks would be cleared in the community(ii)If , then implying that the system passes through a bifurcation state(iii)If , then suggesting that the health risk problem would spread out in the system

4. Discussion and Conclusion

In this work, we propose and analyze the fuzzy SPLMAR model with the aim of using fuzzy set theory in the model system. Although all the parameters associated in the model system are important, in an uncertain environment, we consider the selected two most important parameters and representing the probability that a susceptible individual will drink alcohol after prolonged contact with drinking individuals and the additional death rate induced by alcoholism, respectively. In particular, the two parameters and are considered fuzzy numbers and functions of the degree of force of influence, whose membership functions are well defined. Later, the fuzzy basic risk reproduction number is computed.

Since the basic risk reproduction number is directly related to varied factors of drinking settings, obtained by considering the classical model system (1) with constant parameters and , it implied that the variation of the risky population is always positive. However, when is considered a fuzzy set, it provide us with an additional information regarding the health risk dynamics. The health risks associated with alcoholism can be effectively controlled by controlling the fuzzy basic risk reproduction number . The analysis of the fuzzy model suggests that can be reduced by increasing the value of . This may be done through provision of public health education which increases the resistance of susceptible individuals and prevent them from initiation of drinking behavior which would accelerate their health condition into more risky states. It is generally observed that if the amount of degree of peer of influence to an individual is low, then the alcohol-related health risks in the community may not be the most effective. The reason is natural factors in the community provide virtual immunity from deeply engaging into alcoholic behaviors. However, when the perceived most influential people in the community engage in alcoholism, it implies an increase in the degree of force of influence, and as such, the system will be endemic.

In a future work, some other parameters may be considered with different functionalities under the uncertainty environment. However, the present model can be applied to those types of diseases or conditions which spread through direct contact between susceptible and infected individuals.

Data Availability

No data is included.

Conflicts of Interest

The authors declare that they have no conflicts of interest.