Research Article  Open Access
Dynamical Model about Rumor Spreading with Medium
Abstract
Rumor is a kind of social remark, that is untrue, and not be confirmed, and spreads on a large scale in a short time. Usually, it can induce a cloud of pressure, anxiety, and panic. Traditionally, it is propagated by word of mouth. Nowadays, with the emergence of the internet, rumors can be spread by instant messengers, emails, or publishing. With this new pattern of spreading, an ISRW dynamical model considering the medium as a subclass is established. Beside the dynamical analysis of the model, we mainly explore the mechanism of spreading of individualstoindividuals and mediumtoindividual. By numerical simulation, we find that if we want to control the rumor spreading, it will not only need to control the rate of change of the spreader subclass, but also need to control the change of the information about rumor in medium which has larger influence. Moreover, to control the effusion of rumor is more important than deleting existing information about rumor. On the one hand, government should enhance the management of internet. On the other hand, relevant legal institutions for punishing the rumor creator and spreader on internet who can be tracked should be established. Using this way, involved authorities can propose efficient measures to control the rumor spreading to keep the stabilization of society and development of economy.
1. Introduction
Rumor is a kind of social phenomenon that a remark spreads on a large scale in a short time through chains of communication and runs through the whole evolutionary history of mankind [1]. Narrowly speaking, rumor is defined as untrue information, but has not been confirmed. Usually, it is dispersed by some people in order to achieve the specific purpose: slandering others, manufacturing momentum, diverting attention, causing panic, and so on [2–4]. Most rumors induce panic psychology or economic loss in the accompanying unexpected events. Emergencies cause serious negative impacts on people’s life in several ways: not only the event itself might lead to financial loss or personal injuries, but also the rumor might lead to panic feelings and irrational behavior [5]. For example, the nuclear leakage in Japan caused an altbuying frenzy in China. With the rumors spreading, this frenzy swept and caused social panic in just a few days and the preternatural rising of the salt price, which has a negative effect on society and economy [6]. Recently, there appeared a hearsay that beef in Suzhou city was infected bacillus anthracis; but after a few days, the rumor cannot stand and is confirmed.
The spreading of rumor is in many ways similar to the spreading of epidemic infection, by the spreader or infectious to notify or infect the susceptible. After notified or infected, that the spreader can become the stifler is similar as the infectious can recover after a few time. Thus many epidemic models have been used to describe the spread of information in the form of rumors [6–11]. A classical model for the spread of rumor was introduced by Daley and Kendall [12] who first proposed the basic DK model of rumor spreading. Maki and Thompson [13] focused on the analysis of the rumor spreading model, which was based on mathematical theory and developed the MK model. The DK and MK models have been used extensively for quantitative studies of rumor spreading [14–16]. The literatures [5, 17, 18] give a detailed introduction for the models about rumor spreading.
Traditionally, rumors are propagated by word of mouth. In this propagation mode, the spreading speed is slower and can produce a small effect on society stability. Nowadays, with the emergence of the internet, rumors can spread rapidly by instant messengers, emails, or publishing bloggings that provide faster velocity of transmission and result in widespread of information [6]. Under this case, it is more difficult to control the rumor spreading. The paper about infective vector can be seen in [19, 20]. With this new pattern of rumor spreading, we establish an dynamical model including spreading between individuals and mediumtoindividuals to describe more accurately the actual pattern of transmission, which has not been studied in previous papers. Then the mechanism under spreading between individuals and by medium can be investigated by resorting to the model. Furthermore, we also give the main influence factors of transmission to government that can propose efficient measures to keep the stabilization of society and development of economy.
2. A Dynamical Model for Rumor Spreading
The populations at any time we consider are denoted by , which is classified into three subclasses: the susceptible or ignorant, spreader, and stiflers denoted by , and , respectively. Besides the rumor spreading from to , the susceptible can be notified by medium, such as website news, microblog, Baidu Post Bar, community sites, Renren, and mobile phone. The quantity of rumor in medium is denoted by . The mathematical model to study the rate of change of the populations is given by system (5). Our assumptions on the dynamical transmission of rumor among humans are demonstrated in the following flowchart (Figure 1).
When a susceptible contacts a spreader or sees rumor in medium, the susceptible can get to know about the rumor. Thus, the change over time in the size of the susceptible population only consists of the output, infected by the rumor . For the incidence rate under direct contact between individuals, we adopt the standard incidence rate. According to the indirect infection between the susceptible individual and the rumor in medium, the probability to be infected for a susceptible individual is proportional to the quantity of rumor in medium. So, the indirect transmission rate is . With regard to , it is very difficult to determine the quantity in medium. In addition, for the susceptible, after notified rumor, they do not completely become the rumor, but can doubt its credibility and become a stifler. Thus, the susceptible notified rumor enters the spreader at a constant rate per unit time. Consequently, the change over time in the size of the susceptible dairy cattle population is as follows: The interpretation for the change over time in the size of the spreader population is as follows. Once infected, the individuals will enter the spreader compartment at a constant rate per unit time. However, when two spreaders contact, both of them transmit the rumor, and the spreader can become bored and loses interest in the rumor, and consequently both become stiflers. In this case, the term needs 2 times. During a unite time, the number of the spreader also can realize truth with the development of the matter. When a spreader contacts a stifler and transmits the rumor, the stifler tries to remove it, and the spreader becomes the stifler at a rate of :
On the one hand, there are from the susceptible department to the stifler department. On the other hand, there are , , and from the spreader department to the stifler class. Thus, the change over time in the size of the stifler population is as follows:
For the rumor in medium, they are discharged by the spreader, and the quantity discharged per unite time is assumed as . However, they have the submerged rate of message . The quantity of rumor in medium is close to the follows of a person in network, so its unit also can be individual. Thus, the change over time in the size of the is as follows: Combining the above equations, the model we study in this paper to govern the temporal evolution of rumor during human population is presented in (5), whose parameters are described in Table 1. In the next section, we will discuss the dynamical behavior of model

3. Dynamical Behavior of Model (5)
Because , is independent of time and is a constant, denoted by . We can obtain the equilibrium of the system (5) , where ; that is, the rumor must disappear with time, and all and , that satisfy , represent stable situations. Next, we have a look at the stability of the equilibrium .
Theorem 1. The equilibrium is globally asymptotically stable.
Proof. In order to verify the stability of these points, the linear part of the vector field around them is given by the Jacobian at the equilibrium as follows:
With regard to the above matrix, we know that one zero eigenvalue corresponds to the fact that the order of the dynamical system is two. The other zero eigenvalue is related to the stable center manifold that is the straight line on the plane. Thus, the stability of the equilibrium depends on the signal of the remainder eigenvalues. The corresponding characteristic polynomial can be written as
Using and the Hurwitz criterion, we can get when , is locally asymptotically stable. In fact, applying the twodimensional figures, we can know that when is always less than . Here, we give figures about , and , seen in Figure 2.
From Figure 2, we can see when , the red curve, representing , is a constant 5000000, and always over the blue line, . At , the two curves are very closer, but the gap can be observed by observing the microdiagram round . When , the gap between them becomes larger. At , we have a look at the change of and in terms of the parameter . It can be observed that, at , two curves are the closet, but the red line is still over the blue one. In addition, we also consider the situations of and at and , which show that the situation at is the same as , and the has very small effect on and when and are too small. This also illustrates that the effect of and is the biggest. The situations of other parameters are not given in here for their influences are very little. So, the equilibrium is locally asymptotically stable.
Furthermore, by observing the equations in system (5), we find that as , is monotone decreasing and monotone increasing. The monotonicity of is uncertain, but they are linear. From that, it is impossible to show limit cycle, from which the global stability of can be obtained.
(a)
(b)
(c)
(d)
We know the monotonicity of is uncertain, depending of the parameter values, which have a certain meaning that will be interpreted later. At the end, the rumor will disappear by undergoing for a short time, and a portion of becomes . Hence, the proportion of all spreader in the process of rumor spreading that had spread rumor to others or made statement about rumor in network should attract the attention of the public, which is called as the final scale of spreader, an important indicator of the severity of rumor spreading. Now, we give the expression for the final scale of spreader, denoted by .
Theorem 2. The final scale of spreader is
Proof. In the initial stages, , so we assume and , and at the end of rumor spreading, . From the first equation in system (5), we have From the last equation in system (5), it can be obtained Combining (10) and (11) and boundary conditions applied at 0 and , the following expression can be obtained: Substituting the relation , we can have the expression of the final spreader scale: On the other hand, we can have directly another expression about the final scale from the last equation in system (5): Due to the complexity of equations in model (5), we cannot give the expression only concluding . Comparing (8) with (9), Equation (8) is simpler and more convenient to use to assess parameters sensitivity.
4. The Interpretation of Model Parameter Values
In this section, we give the basis of parameter values in Table 1. Having no real data during the rumor spreading, we assume their values according to individual experience. Firstly, we assume that the average number of acquaintances one person meets per day is every day. When a spreader meets a susceptible, the probability that the spreader talks about the rumor and the susceptible is infected by the rumor is every day. After influenced by the rumor, half of the susceptible believe this rumor and begin to spread it, but the other half become the stiflers, not to spread this rumor. Thus, . In other words, a third of people have no response, and a third of people begin to spread the rumor, and the rest, a third of people, become stiflers. So we get . When two spreaders contact, both of them discuss the rumor, and they lose interest in it for hearing it again and again. We assume the rate of losing interest is every day, so . When a spreader contacts a stifler, the spreader will transmit the rumor, but the stifler tries to remove it. Assuming the probability that the spreader becomes a stifler is , so . As time goes on, the spreader could be aware of the falsity content of the rumor by obtaining involved knowledge. As a result, the spreader also becomes a stifler at a rate of 0.1. The most important and key point is disseminating of the rumor by the medium. Because the quantity of message in medium cannot be measured accurately, we can assume of spreaders will disseminate the rumor into internet, and a person has 500 follower, then the rate that a spreader can disseminate the rumor through the medium is 5. Moreover, due to the update of message on internet, the rumor will be hidden after 12 hours. So the submerged rate of message is 2 every day. With regard to the initial value adopted in this paper, we consider the total number of people is 5000010, where there are 10 spreaders, and .
5. Result and Sensitivity Analysis
Taking no media message into account, we know that, when , will directly decrease to 0, then the rumor cannot spread to others. When , will first increase to its maximum at , then decreases to 0. gives the threshold to justify whether the rumor spread to others can be obtained, which can be visually seen in Figure 3.
(a)
(b)
(c)
However, taking the medium into consideration, the situation becomes more complex. cannot give the threshold of rumor spreading, which needs to add the condition . For example, when and , . However, , and the rumor can be disseminated out, seen in Figure 4. In fact, to control is very difficult. In other words, in this case, it is more difficult to control the spreading of rumor. Moreover, because medium can enlarge the effect of rumor, the introducing of the medium leads to the big change of the solutions behavior, and even the transmission rate is very small. We can draw the influence of medium on the rumor spreading in Figure 5. When that is very small, all people will know the rumor, and the peak of the proportion of spreader can arrive to about 0.15 from 0.05. Moreover, observing the axis, the transmission speed becomes quicker, and the time of duration of rumor can be shortened.
(a)
(b)
What we are more interesting in is to find the key factors to affect the rumor diffusing. Firstly, we carry out sensitivity analysis of the spreader density about some parameters in Figure 6. From Figures6(a)6(b), we can compare the influence of with . When increasing 5 times to parameters , the change of will arrive 1.5 times, but for it has very little change. From that, we can know that the role of the medium is more important during the transmission of rumor. Moreover, Figures 6(c)6(d) have shown that reducing or increasing can reduce the peak of the proportion of spreader and delay the arrival time of the peak. However, has larger influence than on the fraction of spreader. That means the prevention of rumor releasing is more important than the eliminating of existing rumor. So, the enhancing managing network is a vital measure. To look at the influence of and , we give a detailed discussion in Figure 6(e). With the decreasing of , the reduction rate is decreasing, and the effect of becomes smaller; that is, when is very small, the sensitivity and effect of becomes smaller. In turn, when the becomes bigger, the influence of also becomes smaller. Therefore, reducing or increasing is good measure to control rumor spreading, but the effects are not always obvious, which need detailed real data to confirm. Besides, it is remarkable that changing the , , , and can reduce the peak, but the width of cure becomes bigger. So, there are more questions that need indepth discussions in future work.
(a)
(b)
(c)
(d)
(e)
6. Discussion
With the new pattern of rumor spreading corresponding to the appearance of internet, we establish an dynamical model considering the medium as a subclass. Beside the dynamical analysis of the model, that is, the stability of the equilibrium and final spreader scale, we mainly explore the mechanism of spreading between individuals and from mediumtoindividual. By numerical simulation, we find that, due to the introduction of the medium, the mechanism becomes complex. Controlling the rumor spreading need to not only control the derivative of the spreader subclass, but also control change rate of the information about rumor in medium. However, it is difficult to control the information about rumor in medium. By comparing the transmission rate between individuals with the transmission rate of mediumtoindividuals, we find the effect of is larger than for can enlarge the effect of rumor, seen in Figure 6. From Figure 5, we also can see when is very small, its influence is bigger relatively. It is thus clear that consideration of medium during the rumor spreading is necessary. Besides, in Figures 6(c) and 6(d), we find has a larger influence than . That is to say, controlling the effusion of rumor is more important than deleting existing rumor. Thus, on the one hand, government should enhance the management of internet, such as the establishing of realname registration system. On the other hand, relevant legal institutions for the rumor creator and spreader on internet that can be tracked in network should be established. If obtaining all real data about rumor spreading, we can apply the model to quantify various measures, which is very helpful. Resorting to the methods proposed in this paper, the authorities can propose efficient measures to control the rumor spreading to keep the stabilization of society and development of economy.
The model we consider in this paper is deterministic, which means it is a fixed system. Once the same initial conditions and the same parameter values are given, the same trajectory is always obtained. Obviously, there exists a certain gap with the real situation. Clearly, the result is just a situation which appears with a certain probability. So some papers adopted the stochastic model to investigate the transmission process [21–24]. During the rumor spreading, the behavior of individuals may have a certain randomness, that is, the number of contact with others and medium every day and the frequency of being on the net are not constant. Therefore, that the parameters , and are added with a noise (randomness) is more rational. In this case, the steady state is not fixed and is an oscillating curve. So the deterministic model is a mean situation, and the stochastic model is more detailed. The more detailed discussion will be our future work.
Besides, our model is built based on mean and rough conditions. In fact, during rumor spreading there are more human behavior, such as forgetremember mechanism [16] which is not considered in our model. Moreover, the contact number and followers of every people should be different, but it is assumed to be constant in this paper. In emergency situation, authority plays a role of manager [5] that usually publishs some statement to clarify rumors. In this paper, we do not consider it in the way that can confirm rumor and protect people from the rumor infection. These are the directions of research in our future work. In a word, this paper is just an initial working, and lots of work need more prefect models to study.
References
 K. Kawachi, M. Seki, H. Yoshida, Y. Otake, K. Warashina, and H. Ueda, “A rumor transmission model with various contact interactions,” Journal of Theoretical Biology, vol. 253, no. 1, pp. 55–60, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 J. Kostka, Y. A. Oswald, and R. Wattenhofer, “Word of mouth: rumor dissemination in social networks,” in Structural Information and Communication Complexity, vol. 5058 of Lecture Notes in Computer Science, pp. 185–196, Springer, Berlin, Germany, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 Z. L. Zhang and Z. Q. Zhang, “An interplay modal for rumour spreadingand emergency development,” Physica A, vol. 388, no. 19, pp. 4159–4166, 2009. View at: Publisher Site  Google Scholar
 M. Kosfeld, “Rumours and markets,” Journal of Mathematical Economics, vol. 41, no. 6, pp. 646–664, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. J. Zhao, Q. Wang, J. J. Cheng et al., “The impact of authorities' media and rumor dissemination on theevolution of emergency,” Physica A, vol. 391, no. 15, pp. 3978–3987, 2012. View at: Publisher Site  Google Scholar
 L. J. Zhao, J. J. Wang, Y. C. Chen et al., “SIHR rumor spreading model in socialnetworks,” Physica A, vol. 391, no. 7, pp. 2444–2453, 2012. View at: Publisher Site  Google Scholar
 A. Sudbury, “The proportion of the population never hearing a rumour,” Journal of Applied Probability, vol. 22, no. 2, pp. 443–446, 1985. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. Kawachi, “Deterministic models for rumor transmission,” Nonlinear Analysis: Real World Applications, vol. 9, no. 5, pp. 1989–2028, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Rapoport, “Spread of information through a population with sociostructural bias. I. Assumption of transitivity,” Bulletin of Mathematical Biophysics, vol. 15, pp. 523–533, 1953. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. Rapoport, “Spread of information through a population with sociostructural bias. II. Various models with partial transitivity,” Bulletin of Mathematical Biophysics, vol. 15, pp. 535–546, 1953. View at: Google Scholar  MathSciNet
 A. Rapoport and L. I. Rebhun, “On the mathematical theory of rumor spread,” Bulletin of Mathematical Biophysics, vol. 14, pp. 375–383, 1952. View at: Google Scholar  MathSciNet
 D. J. Daley and D. G. Kendall, “Epidemics and rumours,” Nature, vol. 204, article 1118, 1964. View at: Publisher Site  Google Scholar
 D. P. Maki and M. Thompson, Mathematical Models and Applications, PrenticeHall, Englewood Cliffs, NJ, USA, 1973. View at: MathSciNet
 B. Pittel, “On a DaleyKendall model of random rumours,” Journal of Applied Probability, vol. 27, no. 1, pp. 14–27, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 C. Lefèvre and P. Picard, “Distribution of the final extent of a rumour process,” Journal of Applied Probability, vol. 31, no. 1, pp. 244–249, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Gu, W. Li, and X. Cai, “The effect of the forgetremember mechanism onspreading,” European Physical Journal B, vol. 62, no. 2, pp. 247–255, 2008. View at: Publisher Site  Google Scholar
 L. A. Huo, P. Q. Huang, and X. Fang, “An interplay model for authorities' actions and rumor spreading in emergency event,” Physica A, vol. 390, no. 20, pp. 3267–3274, 2011. View at: Publisher Site  Google Scholar
 L. J. Zhao, Q. Wang, J. J. Cheng, Y. C. Chen, J. J. Wang, and W. Huang, “Rumor spreading model with consideration of forgetting mechanism: a case of online blogging live journal,” Physica A, vol. 390, no. 13, pp. 2619–2625, 2011. View at: Publisher Site  Google Scholar
 K. L. Cooke, “Stability analysis for a vector disease model,” The Rocky Mountain Journal of Mathematics, vol. 9, no. 1, pp. 31–42, 1979. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Wang, Z. Jin, Z. M. Yang, Z. K. Zhang, T. Zhou, and G. Q. Sun, “Global analysis of an SIS model with an infective vector on complex networks,” Nonlinear Analysis: Real World Applications, vol. 13, no. 2, pp. 543–557, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Q. Sun, L. Li, Z. Jin, and B. L. Li, “Effect of noise on the pattern formation in an epidemic model,” Numerical Methods for Partial Differential Equations, vol. 26, no. 5, pp. 1168–1179, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Krstić, “The effect of stochastic perturbation on a nonlinear delay malaria epidemic model,” Mathematics and Computers in Simulation, vol. 82, no. 4, pp. 558–569, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Y. Ji, D. Q. Jiang, Q. S. Yang, and N. Z. Shi, “Dynamics of a multigroup SIR epidemic model with stochastic perturbation,” Automatica, vol. 48, no. 1, pp. 121–131, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Mendez, D. Campos, and W. Horsthemke, “Stochastic fluctuationsof the transmission rate in the susceptibleinfectedsusceptible epidemicmodel,” Physical Review E, vol. 86, no. 1, Article ID 011919, 8 pages, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Xiaxia Zhao and Jianzhong Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.