The propagation of rumors is a typical form of social interaction that exerts a powerful influence on human affairs. Rumors represent unproven expositions about or interpretations of news, events, or problems that are of public interest. Because rumors are unconfirmed information, it is hard to determine whether they are true or false  . Traditionally, rumors are propagated by word of mouth    . Nowadays, with the emergence of the internet, rumors spread by instant messengers, emails or publishing blogging and social media that provide a faster speed of transmission   . Once a rumor spreads, the truth is at risk of being distorted in the public sphere. Many case studies indicate that most rumors are false news, which exert negative impacts on the public and on social security  . Rumour is a generally circulated story, report or statement without facts to confirm its truth; it is like a virus, once it is transmitted to an individual, an outbreak of events happens in a short period of time. It is thus that spreading a rumor is like a contagious disease that does not need more than a few infected individuals to infect the entire population. If a rumor starts in a population, sooner or later everyone will know it  .
In recent years, online social media is growing rapidly. Social media provides a convenient communication scheme for people. Meanwhile, the scheme enables unreliable sources to spread large amounts of unverified information among people      . Rumors are thus possible to spread more quickly and widely through online social media compared to traditional offline social communities. The wide spread of misinformation may bring disorder to people especially when they are facing religion or political crises. This indicates that it is crucial for social media to identify misinformation in time so as to limit the spread of rumors  .
At the beginning, mathematical models for the rumors were considered merely speculative and imprecise, but for the fact that rumor spreading is now seen like the transmission of disease  , it leads to a deeper understanding of the future spread of rumor. Rumor transmitted diseases are particularly well suited for modelling because the route of contact is clearly defined: it is through human activities   . In addition, there are so many questions that come to play when it comes to addressing the spread of rumor. Modeling is the only practical approach. The spreading of rumor is in many ways similar to the spreading of epidemic infection by the spreader or the infections to notify or infect the susceptible      . After notified or infected, that the spreader can become the stifler is similar as the infections can recover after sometimes  .
2. The Rumor Models
Rumor can be viewed as an infection of the mind     . Therefore, the rumor propagation problem can be studied with many modeling techniques used in the study of epidemics. Some of these approaches include deterministic models, stochastic models and complex network  .
Different possible behaviors in the dynamics of rumor spreading were studied by  . In their study, they proposed a deterministic propagation rumor model and detected relevant conditions, derived from local stability analysis of the rumor free equilibrium and rumor-endemic equilibrium.
Rumor propagation through different types of mathematical models was investigated by  . In their study, they consider rumor propagation with truth-spreading and determine the threshold which governs the dynamics of the
The dynamics of rumor model in complex heterogonous networks was analyzed by  . They first introduced a useful stochastic method that allows them to obtain meaningful time profiles for the quantities characterizing the propagation process.
3. Model Formulation
Following the classical assumption on the variables described in Table 1, we formulate a deterministic mathematical model to describe the propagation dynamics of rumors. In the model, the total population is divided into four epidemiological classes of ignorant , spreaders targeting community through media , spreaders targeting community through verbal communication and stiflers . In the formulation of the model, we wish to clearly state some essential assumptions that:
1) Recruitment into the ignorant compartment is in constant rate,
2) There is no movement between the classes of spreaders,
3) Ignorant individual become spreader after interaction with spreader through media or spreader through verbal communication,
Table 1. Variables and parameters of the model and their descriptions.
4) The probability of an ignorant to become a spreader through media or verbally are equal,
5) The spreading process are mutually exclusive event,
6) Spreaders through media and verbal communications become stifler at the rate and respectivelly,
7) Proportion of ignorant will interact with the spreaders at the rate and become stifler.
From the model diagram (Figure 1), we have ignorant class to be individuals who are not aware of the rumor. The model is structured such that ignorant individual can either be infected with rumor by individual in or . Note that, proportion of ignorant individual can easily progress to without spreading the rumor. The spreader class are individual who are aware of the rumor and are actively spreading it through media. The spreader class are individuals who have been informed about the rumor and they are actively spreading it orally to the community. The stifler , are individuals who are individuals who are aware of the rumor, spent time spreading it but are no longer spreading the rumor due to some reasons and they also discourage the spread of the rumor.
Thus, the differential equations for the deterministic model are as follow:
4. Basic Properties
Since the model (1.1) to (1.4) monitors human populations, all the variables and the associated parameters are non-negative at all time. It is important to show that the model variables of the model remain non negative for all non-negative initial conditions.
Lemma 1: The region is positively invariant and attract all solutions in
Adding all the equations from (1.1) to (1.4), gives the rate of change of the total human population
Figure 1. Model flow diagram.
Since whenever , then , implying is bounded by .
Thus, a standard comparison theorem by  can be used to show that
In particular, if . Thus, R is positively invariant (i.e. all solution in D remain in D for all time). Furthermore, if then either the solution enters R in finite time or approaches and the spreader variables M and G approaches zero. Hence, D is attracting (i.e. all solution in eventually approach, enter or stay in D). Therefore, the model is epidemiologically and mathematically well posed since all the variables remain non-negative for all . Hence it is sufficient to study the dynamics of the system (1.1) to (1.4) in D.
5. Analysis of the Model
In this section, we discuss the existence and uniqueness of Rumor Free Equilibrium (RFE) of the model and its analysis. The model Equations (1.1) to (1.4) has an RFE given by
The local stability of RFE given will be investigated using the next generation matrix method. We calculate the next generation matrix for the system of the question (1.1) to (1.4) by enumerating the number of ways that: 1) new spreaders arise 2) number of ways that individuals can move but only one way to create a spreader. So, let
F = rate of appearance of new spreaders into the compartment and,
V = rate of transfer into (out) of compartment
Hence the NGM with large domain is two dimensional and is given by
The dominant eigenvalue of (1.6) is equal to , therefore we evaluate the characteristic equation of (1.6) by using , which gives
The jacobian of (1.1) to (1.4) at the equilibrium point is
Now we try to calculate the eignvalues of (1.7) by finding the characteristic equation using the formula
Solving (1.8), we have
Theorem 1: The rumor-free equilibrium of the model equation (1.1) and (1.4) given by , is locally asymptotically stable if and unstable if . Thus this theorem 1 implies that for any given rumor in a population, it can be eliminated when .
Having and to be negative, we also see that and are both negative too when and . Since all the eigenvalues of (1.8) have negative real parts when and , we conclude that rumor-free equilibrium is locally asymptotically stable.
Stability Analysis of the Rumor Endemic Equilibrium
When rumor persists in a population (i.e. at least or ), the model question (1.1) to (1.4) has two equilibrium points denoted by
called rumor verbal-endemic equilibrium and rumor media-endemic equilibrium point. For the existence and uniqueness of and their coordinate has to satisfy the following , , and , , , respectively.
Equating (1.1) to (1.4) to zero, we obtained
From (1.10), we have
Substituting (1.13) and (1.14) in (1.9), we have
Using (1.13), (1.14) and (1.15) in (1.12), we obtained
Local Stability of rumor verbal-endemic equilibrium
We used the Jacobian Stability approach to prove the stability of the rumor verbal-endemic equilibrium.
The Jacobian of (1.9) to (1.12) at the equilibrium point is
If we evaluate (1.17) at and find the eigenvalue using characteristic equation , we will obtain
From the above questions, it is observed clearly that only and , therefore we conclude that rumor verbal-endemic equilibrium is unstable since some of the eigenvalue are greater than zero. Furthermore, once a rumor spreads, the truth is at risk of being distorted in the public sphere, therefore it must spread first before it decline naturally or by counter rumor.
Rumor media-endemic equilibrium
To obtain the rumor media endemic equilibrium we use the same approach as in rumor verbal-endemic equilibrium.
From (1.11), we have
Substituting (1.18) and (1.19) in (1.9), we have
Using (1.18), (1.19) and (1.20) in (1.12), we obtained
Local Stability of Rumor Media-endemic Equilibrium
We also use the Jacobin stability approach to prove the stability of the rumor media endemic. The Jacobian of (1.9) to (1.12) at equilibrium point is
If we evaluate (1.22) at and find the eigenvalue using characteristic equestion , we will obtain
Since some of the eigenvalues are not less than zero, we conclude that rumor through media-endemic equilibrium is unstable. Consequently, when rumor occurs in an ignorant population it must spread first before it declines either naturally or by counter rumor.
6. Numerical Solution and Discussion of Results
This section deals with the numerical studies of the developed model and discussion of the simulation result using estimated parameters and adopted some from other models.
We performed some numerical experiments using ode45 function from MATLAB to study the behavior of the systems on the ignorant, spreaders through media (spreader 1), spreader through verbal communication (spreader 2) and stiflers population. The parameter values used in the simulation of this model are presented in Table 2.
The results of the simulation of the model (1.1) to (1.4) with parameter value in Table 2 are presented in the figures below: (Figures 2-7).
Experiment 1: Investigating the impact of contact rates between the ignorant and spreaders.
Experiment 2: Impact of rate at which spreaders become stiflers.
Figure 2. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when and .
Figure 3. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when and .
Figure 4. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when , .
Figure 5. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population and .
Figure 6. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when , .
Figure 7. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when , .
Table 2. Estimated and adopted parameter values and variables used in the simulation of the model.
7. Discussion of Results
Here we discuss the results obtained in the analytical studies and the numerical experiments.
The modified model consists of a 4-dimensional system of ordinary differential equations. The Rumour Free Equilibrium was established for the system (1.1) to (1.4). We also obtained the reproduction number of the two spreader comportments using next generation matrix method and we established the stability of the rumor Free Equilibrium and Endemic Equilibrium of the modified model using linearization methods. We observed that all the eigenvalues have negative real parts, implying that the rumor free equilibrium is locally asymptotically stable and unstable otherwise. For the endemic equilibrium point, it is observed that some of the eigenvalues of polynomial (1.17) and (1.22) are positive and we conclude that rumor-endemic equilibrium is unstable (Figure 5-7). That is if there exists a rumor in population, it has to spread through first before it decline naturally or suppress by another rumor.
It is observed that the impact of is the result of the interaction of each parameters and it is easy to decrease by regulating and controlling a few parameters so also . The parameter with the most positive impact on is the rate of contact between the spreaders targeting the community through media and the ignorant. Similarly the most positive impact on is the rate of contact between the spreaders targeting the community through verbal communication and the ignorant. While the parameter that has moderate influence is the recruitment rate. In this study, we proposed two classes of spreaders and we observed that the propagation of rumor through media and verbal communication are mutually exclusive. Therefore, the appropriate increase in the propagation of rumor through media is conducive to the elimination of the propagation of rumor through verbal communication and vice versa (Figures 2 to 7). On the other hand, the propagation of rumor through media and verbal communication coexist for long time under certain condition (Figure 4). At this point, the improvement of the propagation of the rumor through media is also conducive to the propagation of rumor through verbal communication in a certain range (Figures 2 to 7).
Deterministic model on rumor propagation was analyzed to get insight into its dynamical features and to obtain associated spreaders thresholds. Some of the theoretical findings of the study are as follows:
1) rumor spreading has a life span which depends on the type of class of the spreader [see (Figures 2-4)].
2) the analysis carried out shows that the increase in stiflers population increased the rate at which spreaders convinced.
The model can be used to control rumor propagation by varying the parameters in the model, not only that it can also be modified to manage some other social negative behaviors like smoking, drug abuse, corruption, prostitution, gang, etc.
 Zhao, L.J., Wang, X.L., Wang, J.J., Qiu, X.Y. and Xie, W.L. (2014) Rumor-Propagation Model with Consideration of Refutation Mechanism in Homogeneouse Social Networks. Discrete Dynamics in Nature and Society, 2014, Article ID: 659273.
 Fang, J., Edward, D., Yang, C. and Naren, R. (2013) Epidemiological Modelling of News and Rumors on Twitter. Proceedings of the 7th Workshop on Social Network Mining and Analysis, Chicago, Illinois, 11 August 2013, 1-9.
 Auwal, A.M. and Farouk, A.U. (2016) Rumor Spreading Model: Vulnerability and its Psychological Effect in Study Salt and Water Ebola Virus Cure Rumor of August 2014 in Nigeria. Proceedings of 56th ISERD International Conference, Mecca, Saudi Arabia, 16-17 December 2016, 30-38.
 Deering, S. (1988) Multicast Routing in Internetworks and Extended LANs. SIGCOMM '88 Symposium Proceedings on Communications Architectures and Protocols, Stanford, CA, USA, 16-18 August 1988, 55-64.
 Nekovee, M., Moreno, Y., Bianconi, G. and Marsili, M. (2007) Theory of Rumor Spreading in a Complex Social Networks. Physica A: Statistical Mechanics and its Applications, 374, 457-470.
 Zhao, L., Wang, J., Chen, Y., Wang, Q. and Cui, H.X. (2012) SIHR Rumor Spreading Model in Social Networks. Physica A: Statistical Mechanics and Its Application, 2444-2453.
 Piqueira, J.R.C., Castana, M.C. and Monteiro, L.H.A. (2004) Modeling the Spreading of HIV in Homosexual Populations with Heterogeneous Preventive Attitude. Journal of Biological Systems, 12, 439-456.