# A Rumor Spreading Model considering the Cumulative Effects of Memory

Authors: Yi Zhang; Jiuping Xu
Journal: Discrete Dynamics in Nature and Society
 (2015)

---

## Abstract

This paper proposes a rumor spreading model which examines how the memory effects rate changes over time in artificial network and a real social network. This model emphasizes a special rumor spreading characteristic called “the cumulative effects of memory.” A functionp(t) reflecting the cumulative memory effects is established, which replaces the constant rate of memory effects in the traditional model. Further, rumor spreading model simulations are conducted with different parameters in three artificial networks. The results show that all the parameters but the initial memory rate of memory effects function have a significant impact on rumor spreading. At the same time, the simulation results show that the final size of the stiflers is sensitive to the average degree k when it is small but is not sensitive to k when the average degree is greater than a certain degree. Finally, through investigations on the Sina Microblog network, the numerical solutions show that the peak value and final size of the rumor spreading are much larger under a variable memory effects rate than under a constant rate.


---

## Body

## 1. Introduction

An old saying goes that rumors come true after being repeated a thousand times. In real life, if people are unable to distinguish authenticity, many rumors are deemed to be true after a large number of repetitions. When rumors are widely propagated, people tend to believe the rumor, especially if they lack timely real information. Because of the increased presence of online social networks, rumors are no longer spread by word of mouth over a small area but are spread amongst strangers in different regions and different countries, meaning that rumors are being spread faster and wider than ever before. This sustained and rapid spreading of rumors deepens people’s impression about the veracity of the rumor and thus improves the credibility. Rumor spreading, therefore, has the ability to shape public opinion and lead to social panic and instability [1]. For example, the 2011 Tohoku nuclear leakage accidents caused a number of rumors in China. Rumors said that taking materials containing iodine could help ward off nuclear radiation, which led to the fact that many people rushed to purchase iodized salt.Rumor spreading has attracted significant attention from researchers and, as a result, a great deal of research has been done on rumor spreading models. In the early stages, scholars borrowed from epidemic models to describe rumor spreading process [2–4]. A classical model is the SIR model, which studies the dynamic behavior of rumor spreading using an epidemic dynamics SIR model [5]. Based on this research, there have been many more rumor spreading models, most of which have been a variant of the SIR model [6–8]. With the development of network technology, many novel models have appeared inspired by empirical discoveries about network topology [7, 9–11]. At the same time, researchers have started to consider the specific features of rumor spreading in their models. Dodds and Watts [12] studied the effects of limited memory on contagion. Zhao et al. [13, 14] proposed a rumor spreading model which considered remembering mechanisms in homogeneous and inhomogeneous networks.However, most previous studies have discussed memory effects as a constant parameter in the model, but, in reality, people hear rumors many times and so have an accumulation of impressions about the rumors, which changes the probability as to when people become rumor spreaders. Therefore, memory effects have a strong time-dependency. Further, the remembering mechanisms can indicate repeatability, which affects the spreading characteristics of the rumor [13]. Even a small amount of memory can affect the rumor spread in small network sizes [15]. Lü et al. [16] proposed a model considered memory effects on information spread but did not construct a function of time and did not build spreading dynamics equations.In the existing research on rumor propagation, there have been two main classification methods for the population. In the first, the population is divided into the three groups borrowed from the epidemic spreading classifications [5], those that are susceptible (have never heard the rumor), those who are infected (are spreading the rumor), and those who are recovered (have heard the rumor but do not spread it). These classifications have been adopted in most rumor spreading models. In the second classification, however, the population is divided into four groups. This classification adds a new group to the first classification, but this group has been interpreted differently in different studies. In some research, this additional group, which arises from the spreaders, is called the hibernators, those that have heard but forgotten the information [13]. In Lü’s research [16], the additional group is called the known, those who are aware of the rumor but are not willing to transmit it as they are suspicious of its authenticity. An incubation class was suggested in [17], where it was proposed that a susceptible individual first goes through a latent period after being infected before becoming a spreader or a stifler, which is a similar idea to that of Lü. In this paper, we follow a four-group classification including incubation class, the unaware, the spreaders, the lurkers, and the stiflers.In this paper, we study a rumor spreading model with variable rate which considers cumulative memory effects as a function of time, which is more in line with reality. The remainder of this paper is organized as follows. In Section2, first, the rumor spreading model with a variable rate varying over time is described. Then, the mean-field equations are derived, and the dynamic analysis of the model is conducted. In Section 3, numerical simulations demonstrating the dynamics of the established model are conducted on regular networks, random networks, BA networks, and online social networks. We compare our model using different parameters and analyze the model’s impact at different average degree of networks. Then, we make comparisons between a constant versus a variable rate on Sina Microblog. In Section 4, further discussion is presented. Finally, in Section 5, we conclude this paper and provide several avenues for further research.


## 2. A Rumor Spreading Model with Variable Memory Rate

In this model, we construct a function of timet which reflects the accumulation of memory, the speed of memory change, and the importance of an event which triggers rumors. Then, with this function added as a parameter to a differential dynamic spread model, a new rumor propagation model is proposed.


### 2.1. Memory Effects Function

Consider a network withN nodes and E links representing the individuals and their interactions. At each time step, each individual is in one of the following four states:(1)
the unaware: this individual has not yet heard the rumor;(2)
the lurkers: this individual knows the rumor but is not willing to spread it because they require an active effort to discern the truth or falseness of the rumor;(3)
the spreaders: this individual knows this rumor and transmits it to all their contacts;(4)
the stiflers: this individual neither trusts the rumor nor transmits it.People generally hear a rumor after many times, and therefore they get an accumulated impression about the rumor, which means that the probability that people become a spreader changes from “will never believe” to “believes.” This can be described as the cumulative effect of memory, which affects the probability that an individual becomes a spreader from a lurker in the rumor spreading process. In information spreading theory, a function was established which reflected the probability that a person would approve the information at timet after having received the news m times [16]. This function is P(m)=(λ-T)e-b(m-1)+T, where λ=P(1) is the approving probability of the first receipt of the information and T∈(0,1] is the upper bound of the probability indicating maximal approval probability. In fact, as m is a function of time, we can transform this function into a function of time t; that is, the probability that a person approves the rumor at time t can be denoted by a function of time t. From the following analysis, we determine a specific form for the function of timet.First, we analyze the changing process of the lurkers. From a microcosmic point of view, lurkers do not automatically change their states at time stept. Some may become a stifler or a spreader, while others remain lurkers and may become stiflers or spreaders at a later time. We assume that the new lurkers at each time step have a part of the residuals which last until the end of the rumor spreading. This corresponds with the fact that there are always some people who take a long time to change their state in real life. From the above analysis, lurkers become spreaders at a variable probability, denoted by p(t), and become stiflers at the rate of p2, so we can determine the lurkers’ process of change at each time step, as shown in Figure 1. In the following paragraph, we explain this function p(t).Figure 1
The change of the lurker at each step.Ei1 represents new lurkers at the ith time step and Eij (j≠1) represents thejth residual of Ei1. The total number of lurkers is E(i)=Ei1+Ei-1,2+Ei-2,3+⋯+E2,i-1+E1i at the ith time step. For example, E(3)=E31+E22+E13.As the number of times the rumor is received, the probability that a individual agrees to the truth of the rumor grows and infinitely approaches a constant [16]. Thus, as time passes, the number of times the rumor is received for the residual lurkers gradually increases. Because the probability p(t) that an individual becomes a spreader from a lurker is a level that reflects the transformation probability of all lurkers, including the residual old lurkers and the new joined lurkers in each time step, as time passes, the probability increases gradually because of the cumulative effect of memory and infinitely approaches a constant when the accumulated memories achieve a certain degree.In the following, we show how probabilityp(t) becomes larger and infinitely approaches a constant. Figure 2 shows the probabilities that each part of the lurkers becomes spreaders. Firstly, we prove a lemma.Figure 2
The probabilities that each part of the lurkers becomes spreaders;pjh represents the probability that lurkers Ejh become spreaders.Lemma 1.
In the rumor spreading process, there is a time pointt=m, such that p1,k≤p2,k≤p3,k≤⋯≤pi,k  (i,k=1,2,3,…,i+k≤m) and pj,1≤pj,2≤pj,3≤⋯≤pj,h  (j,h=1,2,3,…,j+h≤m) before t=m.Proof.
According to rumor spreading theory, at the beginning of the process, there is an increase in the number of spreaders as the rumor propagates, and, with the further spreading of the rumor, the number of spreaders reaches a peak and then decreases as time goes on. Finally, the number of spreaders is zero and the rumor spreading terminates.
Therefore, with an increase in the spreaders, the probability of the lurkers receiving the rumor increases gradually before the number of spreaders reaches a peak at least. We assume this time point to bet=m. That is, the probability that the lurkers become spreaders increases before t=m. Therefore, there is a time point t=m, such that p1,k≤p2,k≤p3,k≤⋯≤pi,k  (i,k=1,2,3,…,i+k≤m) and pj,1≤pj,2≤pj,3≤⋯≤pj,h  (j,h=1,2,3,…,j+h≤m) before t=m.Becausep(t) is a level that reflects the transformation probability of all lurkers, including the residual old lurkers and the new joined lurkers in each time step, we prove that it has the following properties.Theorem 2.
With an increase in time, the probabilityp(t) becomes increasingly larger and infinitely approaches a constant.Proof.
First, the probabilityp(t) becomes increasingly larger. We use mathematical induction to prove this theorem. According to the lemma.
Fori=1,
(1)p2=p12E12+p21E21E12+E21>p11E12+p11E21E12+E21=p1E12+p1E21E12+E21=p1.
For i=2,
(2)p3=p13E13+p22E22+p(31)E31E13+E22+E31>p13E13+p22E22+p11E31E13+E22+E31>p12E13+p12E22+p1E31E13+E22+E31=p12E13+E22+p1E31E13+E22+E31=p2.
For i=3,
(3)p4=p14E14+p23E23+p32E32+p41E41E14+E23+E32+E41>p14E14+p23E23+p32E32+p11E41E14+E23+E32+E41>p13E14+p13E23+p22E32+p1E41(E14+E23)+E32+E41>p13E14+E23+p22E32+p1E41E14+E23+E32+E41=p3.
Assuming the result is true for i=k, that is,
(4)pk=p1,kE1,k+p2,k-1E2,k-1+p3,k-2E3,k-2⋯+pk,1Ek,1E1,k+E2,k-1+E3,k-2+⋯+Ek,1>pk-1.
To establish the truth ofp(k+1)>p(k), we need to show that
(5)p(k+1)=p1,k+1E1,k+1+p2,kE2,k+p3,k-1E3,k-1⋯+pk+1,1Ek+1,1E1,k+1+E2,k+E3,k-1+⋯+Ek+1,1>p1,k(E1,k+1+E2,k)+p2,k-1E3,k-1⋯+p(1)Ek+1,1(E1,k+1+E2,k)+E3,k-1+⋯+Ek+1,1=pk.
Consequently, with an increase in the time stept, the probability p(t) that a lurker becomes a spreader increases gradually because of the cumulative effect of memory. The probability is that the number of contact times infinitely approaches a constant [16]. Therefore, p(t) infinitely approaches a constant when accumulated memories achieve a certain degree.In light of above analysis, the probabilityp(t) that lurkers become spreaders affected by memory accumulation attth time step is
(6)pt=p-qe-c(t-1)+q,
where p, q, and c are parameters. These three parameters reflect the characteristics of the variable memory effects rate. p is the initial value of the memory effects function at t=1. The parameter p reflects the importance of an event triggering rumors in the spreading process, and it is the initial probability that an individual becomes a spreader. A larger value for p means that the spreaders more easily remember the rumor because the event is probably more important. q≤1-p2, q∈(0,1), is the maximal transformation probability. As time passes, p(t) infinitely approaches q. The parameter c can be regarded as the memory speed; namely, c captures how quickly p(t) reaches the maximum value q. The memory effects rate p(t) is a probability varying over time t. Here, we do not consider interest decay and assume that the time scale for the rumor spreading is much faster than the memory decay. Figure 3 shows the transformation probability p(t) as a function of t, given a different c and a fixed p2=0.1, p(1)=0.05, as shown in Figure 2.Figure 3
Memory effects function whenp2=0.1, p(1)=0.1.


### 2.2. Rumor Spreading Model

Denote byS(t), E(t), I(t), and R(t) the density of the unaware, lurkers, spreaders, and stiflers at time t, S(t)+E(t)+I(t)+R(t)=1. As shown in Figure 4, the rumor spreading rules can be summarized as follows.(1)
Everyone needs time to determine the authenticity of rumor, so an unaware becomes a lurker with a probability 1 when an unaware individual contacts a spreader. The contact probabilityk is decided by the specific network topology. Therefore, the reduced speed of the unaware dS(t)/dt is proportional to the number of the unawareS(t) and the spreadersI(t), so we can derive the differential equation:
(7)dS(t)dt=-kStIt.(2)
A lurker becomes a spreader at the rate ofp(t) and becomes a stifler at the rate of p2, which depends on cognition. For example, some unaware individuals have strong knowledge structures and logical reasoning abilities, so they may have little interest in rumors; thus the probability p2 is large. Because an unaware individual becomes a lurker with a probability 1 when an unaware contacts a spreader, the increased speed of the lurkers dE(t)/dt is equal to
(8)dE(t)dt=kStIt-ptEt-p2Et.(3)
When two spreaders contact each other, both may find the two pieces of information inconsistent, so they stop the spread. When a spreader contacts a stifler, the spreader tries to stop the spread, as the stifler shows no interest in the rumor or denies its veracity. We suppose that the above cases occur at the same probabilityp3 [17]. Therefore, the reduced speed of the spreaders dI(t)/dt is proportional to the number of I(t) and R(t)+I(t), Additionally, a lurker becomes a stifler at the rate of p(t), so we have
(9)dItdt=ptEt-kp3ItIt+Rt.(4)
The increasing speed of the stiflersdR(t)/dt is proportional to the number of existing I(t) and R(t)+I(t) from (3), and a lurker becomes a stifler at the rate of p2 from (2), so we get
(10)dR(t)dt=kp3ItIt+Rt+p2Et.Figure 4
Structure of rumor spreading process.Based on the previous discussion, the model is described using the following system of differential equations:(I)dStdt=-kStIt,dEtdt=kS(t)I(t)-p(t)E(t)-p2Et,dItdt=p(t)E(t)-kp3ItIt+Rt,dRtdt=kp3ItIt+Rt+p2Et,S(0)=S0,E(0)=0,I(0)=1-S0>0,R(0)=0.


### 2.3. Steady-State Analysis of Model

At the beginning of the rumor spreading process, the system has only the unaware and the spreaders and at the end has only the unaware and stiflers as an equilibrium state is reached. So here, we analyze the system’s steady-state.SubstitutingI(t)+R(t)=1-S(t)-E(t) into differential equations (9), as the differential equations (7), (8), and (9) are not related to the R(t), we only consider these three equations. Then the reduced limiting dynamic system is given by
(II)dStdt=-kStIt,dEtdt=kStIt-ptEt-p2Et,dItdt=ptEt-kp3It1-Et-St,S(0)=S0,E(0)=0,I(0)=1-S0>0.From differential equations theory, systems(I) and (II) are homogeneous, which means that analyzing the properties of system (II) is equal to analyzing the properties of system (I). First, some definitions and the required theorem are introduced.Our model is a nonautonomous differential dynamic system, the general form of which is(III)dxdt=ft,x,f(t,0)=0,x∈Rn.Suppose thatI=[t0,+∞], U={x∣x≤h}, and V(t,x) is a continuous differentiable function defined in I×U. W(x) is a continuous differentiable function defined in U.Definition 3 (see [18]).
IfW(x) is a positive (negative) definite function, such that V(t,x)≥W(x)(Vt,x≤-Wx) is true in I×U and V(t,0)=0, V(t,x) is called a positive (negative) definite function in I×U. If V(t,x)≥0  (V(t,x)≤0), V(t,x) is called a positive (negative) semidefinite function.Definition 4 (see [18]).
IfW1(x) is a positive (negative) definite function, such that Vt,x≤W1(x), it is called V(t,x) and has an infinitesimally small upper bound. If W2(x) is an infinite positive definite function, such that V(t,x)≥W2(x), it is called V(t,x) and has an infinite lower bound.Theorem 5 (see [18]).
Suppose thatV(t,x) is a positive definite function in I×Rn, which has an infinitesimally small upper bound and an infinite lower bound and V˙(t,x) is a negative semi-definite, so the zero solutions for system (III) are global uniformly asymptotically stable.Theorem 6.
P
*
=
(
S
*
,
0
,
0
)  (0≤S*<1) is the equilibrium of the system (II).Proof.
Let the right side of each of the differential equations (12) be equal to zero in the system which gives the equation
(11)-kS(t)I(t)=0,(12)kS(t)I(t)-p1E(t)-p2E(t)=0,(13)p1E(t)-kp3I(t)(1-E(t)-S(t))=0.
The feasible region for the equations isR3, and so we study the equations in a closed setA = {(S,E,I)∈R3∣S+E+I≤1,S,E,I≥0}.
From (11) and (12), we can get (p(t)+p2)E(t)=0, so E=0. Substituting E=0 into (13), we have kp3I(t)(1-S(t))=0. Adding (12), we can get I(t)=0. From (13), the system has the equilibrium P*=(S,E,I)=(S*,0,0)  (0≤S*<1).Therefore, we can obtain the equilibrium for the system(I)  Q*=(S*,0,0,R*), where S*+R*=1; that is, the rumor must disappear with time, and all S* and R* that satisfy S*+R*=1 represent stable situations. Next, we look at the stability of the equilibrium P*=(S*,0,0).Theorem 7.
The equilibriumP* is globally uniformly asymptotically stable.Proof.
Let the Liapunov function beV(S,E,I)=F(S)+G(E)+H(I), and taking the derivative of V versus t along the solution for the equations, we have
(14)V′S,E,I=-F′(S)kS(t)I(t)+G′E×kStIt-ptEt-p2Et+H′IptEt-kp3It1-Et-St.
In order to ascertainV(S,E,I) and V′(S,E,I), and taking the values of F′(S)=1, G′(E)=1, and H′(I)=1, we get V(S,E,I)=S+E+I, and
(15)V′S,E,I=-kS(t)I(t)+kS(t)I(t)-p(t)E(t)-p2E(t)+p(t)E(t)-kp3It1-Et-St=-p2E(t)-kp3It1-Et-St≤0.
So we getV(S,E,I), which is positive definite, 0≤V(S,E,I)≤1+S; that is, it has an infinitesimally small upper bound and an infinite lower bound, and V′(S,E,I) is negative semidefinite in the feasible region A={(S,E,I)∈R3∣S+E+I≤1,S,E,I≥0}.
From Theorem5, the equilibrium P* is globally uniformly asymptotically stable.


## 2.1. Memory Effects Function

Consider a network withN nodes and E links representing the individuals and their interactions. At each time step, each individual is in one of the following four states:(1)
the unaware: this individual has not yet heard the rumor;(2)
the lurkers: this individual knows the rumor but is not willing to spread it because they require an active effort to discern the truth or falseness of the rumor;(3)
the spreaders: this individual knows this rumor and transmits it to all their contacts;(4)
the stiflers: this individual neither trusts the rumor nor transmits it.People generally hear a rumor after many times, and therefore they get an accumulated impression about the rumor, which means that the probability that people become a spreader changes from “will never believe” to “believes.” This can be described as the cumulative effect of memory, which affects the probability that an individual becomes a spreader from a lurker in the rumor spreading process. In information spreading theory, a function was established which reflected the probability that a person would approve the information at timet after having received the news m times [16]. This function is P(m)=(λ-T)e-b(m-1)+T, where λ=P(1) is the approving probability of the first receipt of the information and T∈(0,1] is the upper bound of the probability indicating maximal approval probability. In fact, as m is a function of time, we can transform this function into a function of time t; that is, the probability that a person approves the rumor at time t can be denoted by a function of time t. From the following analysis, we determine a specific form for the function of timet.First, we analyze the changing process of the lurkers. From a microcosmic point of view, lurkers do not automatically change their states at time stept. Some may become a stifler or a spreader, while others remain lurkers and may become stiflers or spreaders at a later time. We assume that the new lurkers at each time step have a part of the residuals which last until the end of the rumor spreading. This corresponds with the fact that there are always some people who take a long time to change their state in real life. From the above analysis, lurkers become spreaders at a variable probability, denoted by p(t), and become stiflers at the rate of p2, so we can determine the lurkers’ process of change at each time step, as shown in Figure 1. In the following paragraph, we explain this function p(t).Figure 1
The change of the lurker at each step.Ei1 represents new lurkers at the ith time step and Eij (j≠1) represents thejth residual of Ei1. The total number of lurkers is E(i)=Ei1+Ei-1,2+Ei-2,3+⋯+E2,i-1+E1i at the ith time step. For example, E(3)=E31+E22+E13.As the number of times the rumor is received, the probability that a individual agrees to the truth of the rumor grows and infinitely approaches a constant [16]. Thus, as time passes, the number of times the rumor is received for the residual lurkers gradually increases. Because the probability p(t) that an individual becomes a spreader from a lurker is a level that reflects the transformation probability of all lurkers, including the residual old lurkers and the new joined lurkers in each time step, as time passes, the probability increases gradually because of the cumulative effect of memory and infinitely approaches a constant when the accumulated memories achieve a certain degree.In the following, we show how probabilityp(t) becomes larger and infinitely approaches a constant. Figure 2 shows the probabilities that each part of the lurkers becomes spreaders. Firstly, we prove a lemma.Figure 2
The probabilities that each part of the lurkers becomes spreaders;pjh represents the probability that lurkers Ejh become spreaders.Lemma 1.
In the rumor spreading process, there is a time pointt=m, such that p1,k≤p2,k≤p3,k≤⋯≤pi,k  (i,k=1,2,3,…,i+k≤m) and pj,1≤pj,2≤pj,3≤⋯≤pj,h  (j,h=1,2,3,…,j+h≤m) before t=m.Proof.
According to rumor spreading theory, at the beginning of the process, there is an increase in the number of spreaders as the rumor propagates, and, with the further spreading of the rumor, the number of spreaders reaches a peak and then decreases as time goes on. Finally, the number of spreaders is zero and the rumor spreading terminates.
Therefore, with an increase in the spreaders, the probability of the lurkers receiving the rumor increases gradually before the number of spreaders reaches a peak at least. We assume this time point to bet=m. That is, the probability that the lurkers become spreaders increases before t=m. Therefore, there is a time point t=m, such that p1,k≤p2,k≤p3,k≤⋯≤pi,k  (i,k=1,2,3,…,i+k≤m) and pj,1≤pj,2≤pj,3≤⋯≤pj,h  (j,h=1,2,3,…,j+h≤m) before t=m.Becausep(t) is a level that reflects the transformation probability of all lurkers, including the residual old lurkers and the new joined lurkers in each time step, we prove that it has the following properties.Theorem 2.
With an increase in time, the probabilityp(t) becomes increasingly larger and infinitely approaches a constant.Proof.
First, the probabilityp(t) becomes increasingly larger. We use mathematical induction to prove this theorem. According to the lemma.
Fori=1,
(1)p2=p12E12+p21E21E12+E21>p11E12+p11E21E12+E21=p1E12+p1E21E12+E21=p1.
For i=2,
(2)p3=p13E13+p22E22+p(31)E31E13+E22+E31>p13E13+p22E22+p11E31E13+E22+E31>p12E13+p12E22+p1E31E13+E22+E31=p12E13+E22+p1E31E13+E22+E31=p2.
For i=3,
(3)p4=p14E14+p23E23+p32E32+p41E41E14+E23+E32+E41>p14E14+p23E23+p32E32+p11E41E14+E23+E32+E41>p13E14+p13E23+p22E32+p1E41(E14+E23)+E32+E41>p13E14+E23+p22E32+p1E41E14+E23+E32+E41=p3.
Assuming the result is true for i=k, that is,
(4)pk=p1,kE1,k+p2,k-1E2,k-1+p3,k-2E3,k-2⋯+pk,1Ek,1E1,k+E2,k-1+E3,k-2+⋯+Ek,1>pk-1.
To establish the truth ofp(k+1)>p(k), we need to show that
(5)p(k+1)=p1,k+1E1,k+1+p2,kE2,k+p3,k-1E3,k-1⋯+pk+1,1Ek+1,1E1,k+1+E2,k+E3,k-1+⋯+Ek+1,1>p1,k(E1,k+1+E2,k)+p2,k-1E3,k-1⋯+p(1)Ek+1,1(E1,k+1+E2,k)+E3,k-1+⋯+Ek+1,1=pk.
Consequently, with an increase in the time stept, the probability p(t) that a lurker becomes a spreader increases gradually because of the cumulative effect of memory. The probability is that the number of contact times infinitely approaches a constant [16]. Therefore, p(t) infinitely approaches a constant when accumulated memories achieve a certain degree.In light of above analysis, the probabilityp(t) that lurkers become spreaders affected by memory accumulation attth time step is
(6)pt=p-qe-c(t-1)+q,
where p, q, and c are parameters. These three parameters reflect the characteristics of the variable memory effects rate. p is the initial value of the memory effects function at t=1. The parameter p reflects the importance of an event triggering rumors in the spreading process, and it is the initial probability that an individual becomes a spreader. A larger value for p means that the spreaders more easily remember the rumor because the event is probably more important. q≤1-p2, q∈(0,1), is the maximal transformation probability. As time passes, p(t) infinitely approaches q. The parameter c can be regarded as the memory speed; namely, c captures how quickly p(t) reaches the maximum value q. The memory effects rate p(t) is a probability varying over time t. Here, we do not consider interest decay and assume that the time scale for the rumor spreading is much faster than the memory decay. Figure 3 shows the transformation probability p(t) as a function of t, given a different c and a fixed p2=0.1, p(1)=0.05, as shown in Figure 2.Figure 3
Memory effects function whenp2=0.1, p(1)=0.1.


## 2.2. Rumor Spreading Model

Denote byS(t), E(t), I(t), and R(t) the density of the unaware, lurkers, spreaders, and stiflers at time t, S(t)+E(t)+I(t)+R(t)=1. As shown in Figure 4, the rumor spreading rules can be summarized as follows.(1)
Everyone needs time to determine the authenticity of rumor, so an unaware becomes a lurker with a probability 1 when an unaware individual contacts a spreader. The contact probabilityk is decided by the specific network topology. Therefore, the reduced speed of the unaware dS(t)/dt is proportional to the number of the unawareS(t) and the spreadersI(t), so we can derive the differential equation:
(7)dS(t)dt=-kStIt.(2)
A lurker becomes a spreader at the rate ofp(t) and becomes a stifler at the rate of p2, which depends on cognition. For example, some unaware individuals have strong knowledge structures and logical reasoning abilities, so they may have little interest in rumors; thus the probability p2 is large. Because an unaware individual becomes a lurker with a probability 1 when an unaware contacts a spreader, the increased speed of the lurkers dE(t)/dt is equal to
(8)dE(t)dt=kStIt-ptEt-p2Et.(3)
When two spreaders contact each other, both may find the two pieces of information inconsistent, so they stop the spread. When a spreader contacts a stifler, the spreader tries to stop the spread, as the stifler shows no interest in the rumor or denies its veracity. We suppose that the above cases occur at the same probabilityp3 [17]. Therefore, the reduced speed of the spreaders dI(t)/dt is proportional to the number of I(t) and R(t)+I(t), Additionally, a lurker becomes a stifler at the rate of p(t), so we have
(9)dItdt=ptEt-kp3ItIt+Rt.(4)
The increasing speed of the stiflersdR(t)/dt is proportional to the number of existing I(t) and R(t)+I(t) from (3), and a lurker becomes a stifler at the rate of p2 from (2), so we get
(10)dR(t)dt=kp3ItIt+Rt+p2Et.Figure 4
Structure of rumor spreading process.Based on the previous discussion, the model is described using the following system of differential equations:(I)dStdt=-kStIt,dEtdt=kS(t)I(t)-p(t)E(t)-p2Et,dItdt=p(t)E(t)-kp3ItIt+Rt,dRtdt=kp3ItIt+Rt+p2Et,S(0)=S0,E(0)=0,I(0)=1-S0>0,R(0)=0.


## 2.3. Steady-State Analysis of Model

At the beginning of the rumor spreading process, the system has only the unaware and the spreaders and at the end has only the unaware and stiflers as an equilibrium state is reached. So here, we analyze the system’s steady-state.SubstitutingI(t)+R(t)=1-S(t)-E(t) into differential equations (9), as the differential equations (7), (8), and (9) are not related to the R(t), we only consider these three equations. Then the reduced limiting dynamic system is given by
(II)dStdt=-kStIt,dEtdt=kStIt-ptEt-p2Et,dItdt=ptEt-kp3It1-Et-St,S(0)=S0,E(0)=0,I(0)=1-S0>0.From differential equations theory, systems(I) and (II) are homogeneous, which means that analyzing the properties of system (II) is equal to analyzing the properties of system (I). First, some definitions and the required theorem are introduced.Our model is a nonautonomous differential dynamic system, the general form of which is(III)dxdt=ft,x,f(t,0)=0,x∈Rn.Suppose thatI=[t0,+∞], U={x∣x≤h}, and V(t,x) is a continuous differentiable function defined in I×U. W(x) is a continuous differentiable function defined in U.Definition 3 (see [18]).
IfW(x) is a positive (negative) definite function, such that V(t,x)≥W(x)(Vt,x≤-Wx) is true in I×U and V(t,0)=0, V(t,x) is called a positive (negative) definite function in I×U. If V(t,x)≥0  (V(t,x)≤0), V(t,x) is called a positive (negative) semidefinite function.Definition 4 (see [18]).
IfW1(x) is a positive (negative) definite function, such that Vt,x≤W1(x), it is called V(t,x) and has an infinitesimally small upper bound. If W2(x) is an infinite positive definite function, such that V(t,x)≥W2(x), it is called V(t,x) and has an infinite lower bound.Theorem 5 (see [18]).
Suppose thatV(t,x) is a positive definite function in I×Rn, which has an infinitesimally small upper bound and an infinite lower bound and V˙(t,x) is a negative semi-definite, so the zero solutions for system (III) are global uniformly asymptotically stable.Theorem 6.
P
*
=
(
S
*
,
0
,
0
)  (0≤S*<1) is the equilibrium of the system (II).Proof.
Let the right side of each of the differential equations (12) be equal to zero in the system which gives the equation
(11)-kS(t)I(t)=0,(12)kS(t)I(t)-p1E(t)-p2E(t)=0,(13)p1E(t)-kp3I(t)(1-E(t)-S(t))=0.
The feasible region for the equations isR3, and so we study the equations in a closed setA = {(S,E,I)∈R3∣S+E+I≤1,S,E,I≥0}.
From (11) and (12), we can get (p(t)+p2)E(t)=0, so E=0. Substituting E=0 into (13), we have kp3I(t)(1-S(t))=0. Adding (12), we can get I(t)=0. From (13), the system has the equilibrium P*=(S,E,I)=(S*,0,0)  (0≤S*<1).Therefore, we can obtain the equilibrium for the system(I)  Q*=(S*,0,0,R*), where S*+R*=1; that is, the rumor must disappear with time, and all S* and R* that satisfy S*+R*=1 represent stable situations. Next, we look at the stability of the equilibrium P*=(S*,0,0).Theorem 7.
The equilibriumP* is globally uniformly asymptotically stable.Proof.
Let the Liapunov function beV(S,E,I)=F(S)+G(E)+H(I), and taking the derivative of V versus t along the solution for the equations, we have
(14)V′S,E,I=-F′(S)kS(t)I(t)+G′E×kStIt-ptEt-p2Et+H′IptEt-kp3It1-Et-St.
In order to ascertainV(S,E,I) and V′(S,E,I), and taking the values of F′(S)=1, G′(E)=1, and H′(I)=1, we get V(S,E,I)=S+E+I, and
(15)V′S,E,I=-kS(t)I(t)+kS(t)I(t)-p(t)E(t)-p2E(t)+p(t)E(t)-kp3It1-Et-St=-p2E(t)-kp3It1-Et-St≤0.
So we getV(S,E,I), which is positive definite, 0≤V(S,E,I)≤1+S; that is, it has an infinitesimally small upper bound and an infinite lower bound, and V′(S,E,I) is negative semidefinite in the feasible region A={(S,E,I)∈R3∣S+E+I≤1,S,E,I≥0}.
From Theorem5, the equilibrium P* is globally uniformly asymptotically stable.


## 3. Numerical Analysis

In the following sections, simulations are conducted using the Runge-Kutta method and MATLAB in artificial networks and online social networks, respectively. First, we examine how the proportion of spreaders and stiflers changes over time by varying the function parameters. Then, we discuss the effect of the average degree of networks on rumor spreading. Finally, numerical solutions for the established model are conducted on the Sina Microbolg online social network, and we make comparisons between a constant versus a variable memory rate on the Sina Microbolg.


### 3.1. Sensitivity Analysis in Artificial Networks

To compare the results, the proposed model was performed on two artificial networks with different values for parameterc. We chose regular networks and ER random networks with the same network size N=10000 and the same average degree k=16.(1)
Regular networks [19]: this is a nearest-neighbor coupled network, in which every node is of the same degree.(2)
Random networks [20]: two different nodes are connected with a probability of P=0.0016, so we get random networks with N nodes and PN(N-1)/2 edges. ER random networks degree distributions are an approximate Poisson distribution, k=P(N-1)≈PN. This distribution reaches its peak value at the average degree k [21].(3)
BA network [22]: the network exhibits a power-law degree distribution p(k)~k-γ with γ=3.Given that the other parameters are fixed, we compared the rumor spreading processes on the three networks with different values forc. The smaller the value of c, the slower the speed of increase for I(t). So here, we set c=0.05, c=0.1, and c=0.5. Figure 5 illustrates how the density of spreaders and stiflers changes over time for the different values of c in a regular network. From a macroscopic perspective, we found that the greater parameter c is, the higher the peak for I(t) is, and, as parameter c increases, the size of the stiflers increases, as larger c indicates a higher speed for p(t) to reach the maximum value q. Figure 5(a) describes how the density of spreaders changes with changes in parameter c. The green line indicates a scenario in which the memory speed is very fast. We have carefully checked that the results are not sensitive to c when c≥0.5. It can be seen that the higher parameter c is, the larger the spreader peak value is and the slower the rumor terminates. Figure 5(b) describes how the density of stiflers changes with c over time. The final value for the stifler density R(t) is greater, which indicates the number of people affected by a rumor when c is larger. Clearly, the larger the value of c when other parameters are fixed, the broader the rumor’s influence. A larger c indicates that the p(t) speed is higher and reaches the maximum value, less spreaders change into stiflers, and therefore the influence of the rumor increases.Density of spreaders and stiflers on regular networks.
(a)
(b)Figures6 and 7 illustrate how the density of spreaders and stiflers changes over time with different values for c in ER random networks and BA networks. Similar to the regular network, it can be seen that rumors spread more broadly when c is larger. However, there is a most significant impact from the importance of events on rumor spreading in random networks. That is to say, rumor spreading is sensitive to parameter c, which is related to the topology of random networks. From Figure 7(a), in the same situation, if c is larger, the rumors spread most quickly in BA networks than in the other two networks. And the number of I(t) follows the relation BA>random>regular. These results are consistent with the traditional understanding of epidemic and information spread in that the speed and range of the spread obey the relationship [8]. Figure 7(b) describes how the density of stiflers changes as c changes over time. When parameter c is fixed, we can see that the rumor spreads faster and the final size of the stiflers, R(t), is largest in BA networks. In other words, rumors spread faster and more broadly on BA networks than the other networks. For all the networks, the bigger the value of c, when other parameters are fixed, the bigger the rumors influence.Density of spreaders and stiflers on random networks.
(a)
(b)Density of spreaders and stiflers on BA networks.
(a)
(b)Because BA networks are closer to real social networks, we analyze the effect of the parametersp and q on these networks, respectively. Figure 8(a) illustrates how the spreader density changes as parameter q changes in BA networks. Given that other parameters are fixed, the bigger the value of q, the bigger the rumor’s influence. q is the maximal transformation probability that lurkers become spreaders; in reality, along with the increase of q, there are more lurkers becoming spreaders. As a result, the number of spreaders increases, which increases the influence of the rumor. Figure 8(b) describes how the density of spreader changes as p changes over time. Also, we can see in Figure 8(b) that the bigger the value of p, when other parameters are fixed, the bigger the rumor’s influence. p is the initial probability that an individual becomes a spreader. In reality, a large p means that the individual thinks that an event triggering rumors is very important in the spreading process, which can make more individuals change into spreaders and, therefore, increase the influence of the rumor. But p has little influence on the spread of rumor.Density of spreaders with differentp and q on BA networks.
(a)
(b)Figure9 displays the final spread size for regular networks and random networks as a function of the average degree k when c=0.1. In random networks, because k=P(N-1)≈PN, analyzing the effect of average degree k is equal to analyzing the effect of connected probability P on random network when network size is fixed. There are many similarities between Figures 9(a) and 9(b). It is obvious that R increases with k, but the final spreading size is larger in random networks at same average degree. Moreover, the curves show that R peaks suddenly with small k, which means that R is sensitive when k is small. Rumor spreading enhancement with the average degree could be more significant in random networks than that on regular networks when k is small because of the different network topologies, while for a large k, the changing scale for R is smaller, which means that the final size of R is not sensitive to the average degree when k is large. We have carefully checked the fact that, with an increase in k, R becomes increasingly larger and is finally equal to 1 when k is large enough. That is, when all individuals change to stiflers, the whole population is infected. In general, the final size of R is sensitive to k when k is small but not sensitive to k when k is large.Density of stiflers on regular networks and random networks.
(a)
(b)


### 3.2. Rumor Spreading in Social Networks

To study how rumors spread on social networks, we simulated the rumor spreading process on a snapshot of the China’s largest microblog site, Sina Microblog networks which had been crawled by Yuan and Liu [23]. This snapshot consisted of 118,517 nodes and 2728,213 edges, with the degree distribution obeying the power law p(k)~k-γ [24], γ=1.43. We chose online social networks because of the available network data and because we feel that their structure might be similar to that of other real world social networks.Figure10 corresponds to a case where P=0.05, q=0.65, and c=0.1, which shows the general trends for the four kinds of agents in our rumor spreading model. From the simulation, we find that as the rumor propagates at the beginning of the process, there is a sharp increase in the number of spreaders. With the further spreading of the rumor, the number of spreaders quickly reaches a peak at about t=30 and then decrease as time goes on. Finally, the number of spreaders is zero and the spreading of the rumor terminates. With an increase in time t, the number of unaware individuals S(t) always reduces while the number of stiflers R(t) always increases, finally reaching a balance at the end of the rumor spreading process. From this figure, it can be seen that almost all individuals become stiflers after about 70 time steps. The variation trend for the number of lurkers is similar to that of the spreaders, in that they increase at first and then decrease to zero. However, this increasing and decreasing trend for the lurkers is much more moderate than for the spreaders, and, as can be seen, the peak for the lurkers is much smaller than that for the spreaders.Figure 10
Densities of unaware individuals, spreaders, lurkers, and stiflers over time withc=0.1, p(1)=0.05, q=0.65, p2=0.15, and q=0.17.Figure11 shows how the densities of spreaders and stiflers change over time for variable and constant memory rate in the Sina Microbolg network. For the variable rate, it can be seen that the peak of the spreader is much higher, compared to the constant rate case. It can be seen from Figure 11(b) that final stifler density value R(t) is greater, which means that the population affected by the rumor is greater under a variable memory rate. At the same time, the rumor terminates faster when the memory rate is nonconstant, which shows that the addition of function p(t) has accelerated the pace of the entire rumor spreading process. Generally speaking, the rumor spreads faster and more broadly under variable rates on the Sina Microbolg. In real social networks, the number of people who spread rumors is small, but the rumors can spread widely as social networks have a high clustering coefficient which means the number of the lurkers who come into contact with the rumor is greater. Rumors spread quickly in the early stages through the friends of the people on social networks, and, some time later, the number of people who come into contact with the rumor reaches a maximum, after which people lose interest because local authorities have debunked the rumor or people have an increased awareness of the truth or otherwise of the rumor. Thus the speed of the transmission reduces. Finally, the spreaders gradually disappear, and the entire population becomes stiflers.Density of spreaders and stiflers on Sina Microblog network for the values of the memory rate is variable(p(t)=-0.6e-0.1(t-1)+0.65) and constant (p(t)=0.3), respectively.
(a)
(b)


## 3.1. Sensitivity Analysis in Artificial Networks

To compare the results, the proposed model was performed on two artificial networks with different values for parameterc. We chose regular networks and ER random networks with the same network size N=10000 and the same average degree k=16.(1)
Regular networks [19]: this is a nearest-neighbor coupled network, in which every node is of the same degree.(2)
Random networks [20]: two different nodes are connected with a probability of P=0.0016, so we get random networks with N nodes and PN(N-1)/2 edges. ER random networks degree distributions are an approximate Poisson distribution, k=P(N-1)≈PN. This distribution reaches its peak value at the average degree k [21].(3)
BA network [22]: the network exhibits a power-law degree distribution p(k)~k-γ with γ=3.Given that the other parameters are fixed, we compared the rumor spreading processes on the three networks with different values forc. The smaller the value of c, the slower the speed of increase for I(t). So here, we set c=0.05, c=0.1, and c=0.5. Figure 5 illustrates how the density of spreaders and stiflers changes over time for the different values of c in a regular network. From a macroscopic perspective, we found that the greater parameter c is, the higher the peak for I(t) is, and, as parameter c increases, the size of the stiflers increases, as larger c indicates a higher speed for p(t) to reach the maximum value q. Figure 5(a) describes how the density of spreaders changes with changes in parameter c. The green line indicates a scenario in which the memory speed is very fast. We have carefully checked that the results are not sensitive to c when c≥0.5. It can be seen that the higher parameter c is, the larger the spreader peak value is and the slower the rumor terminates. Figure 5(b) describes how the density of stiflers changes with c over time. The final value for the stifler density R(t) is greater, which indicates the number of people affected by a rumor when c is larger. Clearly, the larger the value of c when other parameters are fixed, the broader the rumor’s influence. A larger c indicates that the p(t) speed is higher and reaches the maximum value, less spreaders change into stiflers, and therefore the influence of the rumor increases.Density of spreaders and stiflers on regular networks.
(a)
(b)Figures6 and 7 illustrate how the density of spreaders and stiflers changes over time with different values for c in ER random networks and BA networks. Similar to the regular network, it can be seen that rumors spread more broadly when c is larger. However, there is a most significant impact from the importance of events on rumor spreading in random networks. That is to say, rumor spreading is sensitive to parameter c, which is related to the topology of random networks. From Figure 7(a), in the same situation, if c is larger, the rumors spread most quickly in BA networks than in the other two networks. And the number of I(t) follows the relation BA>random>regular. These results are consistent with the traditional understanding of epidemic and information spread in that the speed and range of the spread obey the relationship [8]. Figure 7(b) describes how the density of stiflers changes as c changes over time. When parameter c is fixed, we can see that the rumor spreads faster and the final size of the stiflers, R(t), is largest in BA networks. In other words, rumors spread faster and more broadly on BA networks than the other networks. For all the networks, the bigger the value of c, when other parameters are fixed, the bigger the rumors influence.Density of spreaders and stiflers on random networks.
(a)
(b)Density of spreaders and stiflers on BA networks.
(a)
(b)Because BA networks are closer to real social networks, we analyze the effect of the parametersp and q on these networks, respectively. Figure 8(a) illustrates how the spreader density changes as parameter q changes in BA networks. Given that other parameters are fixed, the bigger the value of q, the bigger the rumor’s influence. q is the maximal transformation probability that lurkers become spreaders; in reality, along with the increase of q, there are more lurkers becoming spreaders. As a result, the number of spreaders increases, which increases the influence of the rumor. Figure 8(b) describes how the density of spreader changes as p changes over time. Also, we can see in Figure 8(b) that the bigger the value of p, when other parameters are fixed, the bigger the rumor’s influence. p is the initial probability that an individual becomes a spreader. In reality, a large p means that the individual thinks that an event triggering rumors is very important in the spreading process, which can make more individuals change into spreaders and, therefore, increase the influence of the rumor. But p has little influence on the spread of rumor.Density of spreaders with differentp and q on BA networks.
(a)
(b)Figure9 displays the final spread size for regular networks and random networks as a function of the average degree k when c=0.1. In random networks, because k=P(N-1)≈PN, analyzing the effect of average degree k is equal to analyzing the effect of connected probability P on random network when network size is fixed. There are many similarities between Figures 9(a) and 9(b). It is obvious that R increases with k, but the final spreading size is larger in random networks at same average degree. Moreover, the curves show that R peaks suddenly with small k, which means that R is sensitive when k is small. Rumor spreading enhancement with the average degree could be more significant in random networks than that on regular networks when k is small because of the different network topologies, while for a large k, the changing scale for R is smaller, which means that the final size of R is not sensitive to the average degree when k is large. We have carefully checked the fact that, with an increase in k, R becomes increasingly larger and is finally equal to 1 when k is large enough. That is, when all individuals change to stiflers, the whole population is infected. In general, the final size of R is sensitive to k when k is small but not sensitive to k when k is large.Density of stiflers on regular networks and random networks.
(a)
(b)


## 3.2. Rumor Spreading in Social Networks

To study how rumors spread on social networks, we simulated the rumor spreading process on a snapshot of the China’s largest microblog site, Sina Microblog networks which had been crawled by Yuan and Liu [23]. This snapshot consisted of 118,517 nodes and 2728,213 edges, with the degree distribution obeying the power law p(k)~k-γ [24], γ=1.43. We chose online social networks because of the available network data and because we feel that their structure might be similar to that of other real world social networks.Figure10 corresponds to a case where P=0.05, q=0.65, and c=0.1, which shows the general trends for the four kinds of agents in our rumor spreading model. From the simulation, we find that as the rumor propagates at the beginning of the process, there is a sharp increase in the number of spreaders. With the further spreading of the rumor, the number of spreaders quickly reaches a peak at about t=30 and then decrease as time goes on. Finally, the number of spreaders is zero and the spreading of the rumor terminates. With an increase in time t, the number of unaware individuals S(t) always reduces while the number of stiflers R(t) always increases, finally reaching a balance at the end of the rumor spreading process. From this figure, it can be seen that almost all individuals become stiflers after about 70 time steps. The variation trend for the number of lurkers is similar to that of the spreaders, in that they increase at first and then decrease to zero. However, this increasing and decreasing trend for the lurkers is much more moderate than for the spreaders, and, as can be seen, the peak for the lurkers is much smaller than that for the spreaders.Figure 10
Densities of unaware individuals, spreaders, lurkers, and stiflers over time withc=0.1, p(1)=0.05, q=0.65, p2=0.15, and q=0.17.Figure11 shows how the densities of spreaders and stiflers change over time for variable and constant memory rate in the Sina Microbolg network. For the variable rate, it can be seen that the peak of the spreader is much higher, compared to the constant rate case. It can be seen from Figure 11(b) that final stifler density value R(t) is greater, which means that the population affected by the rumor is greater under a variable memory rate. At the same time, the rumor terminates faster when the memory rate is nonconstant, which shows that the addition of function p(t) has accelerated the pace of the entire rumor spreading process. Generally speaking, the rumor spreads faster and more broadly under variable rates on the Sina Microbolg. In real social networks, the number of people who spread rumors is small, but the rumors can spread widely as social networks have a high clustering coefficient which means the number of the lurkers who come into contact with the rumor is greater. Rumors spread quickly in the early stages through the friends of the people on social networks, and, some time later, the number of people who come into contact with the rumor reaches a maximum, after which people lose interest because local authorities have debunked the rumor or people have an increased awareness of the truth or otherwise of the rumor. Thus the speed of the transmission reduces. Finally, the spreaders gradually disappear, and the entire population becomes stiflers.Density of spreaders and stiflers on Sina Microblog network for the values of the memory rate is variable(p(t)=-0.6e-0.1(t-1)+0.65) and constant (p(t)=0.3), respectively.
(a)
(b)


## 4. Discussion

In our opinion, the cumulative effects of memory play an important role in the rumor spreading process. Although some variants of the SIR model do consider the memory effects, they do not consider the cumulative memory effect. That is to say, they reflect the memory effects using a constant parameter rather than a variable parameter which changes with time across the entire rumor spreading process. In order to better describe the rumor spreading process, we establish mean-field equations with variable memory effects rates in the network. As such, previous models which consider the memory effects rate as a constant become a special case of the new model. We compared our new rumor spreading model to the traditional model on an online social network. The simulation results showed that the functionp(t) induced a broader and faster rumor spreading process. Although rumor spreading simulates on a given snapshot of the Sina Microblog networks. The model does not take into account all rumor spreading features, such as the forgetting mechanism [25], or the nonredundancy of contact [16], but it largely reflects the cumulative memory effects in the variants of the SIR model. We also found that rumor spreading is more influenced by parameter c in random networks, which is caused by the specific topology of random networks. At the same time, it can be seen that the random networks are more effective in rumor spreading which supports Zanette’s experiment. In addition, our experiments show that a high average degree is very helpful in promoting the spreading of rumor. This is expected because as the average degree of networks increases, there are more friends for a user, which means there are more chances for spreaders to contact the unaware, making the spreading of the rumor much easier. However, the final number of people affected by the rumors is sensitive to the average degree k when k is small, but for a large k, nearly all individuals spread the rumor. In fact, the entire population knows the rumor quickly when k reaches a certain level, at which point the increased average degree has no effect.From the discussion of the cumulative effects of memory in this paper, governments should use this special property to control the spread of rumors. Because a sustained and rapid spread of information deepens people’s impression about this information and to a certain extent improves the credibility of the information, when an emergency arises, governments should release official information at a high frequency to ensure the sustained and rapid spreading of official information, which may help avoid the spread of panic and prevent the spread of unfounded rumors.


## 5. Conclusions

In this paper, we have proposed a rumor spreading model which considers changing memory effects over time. At the same time, the speed thatp(t) takes to reach a maximum value and the importance of an event were considered. The cumulative effect of memory reflects the repeatability of the spreading characteristics of the rumor, so the probability p(t) of a lurker becoming a spreader becomes increasingly larger and infinitely closer to a constant. Therefore, the probability is not a constant but a function of time in entire rumor spreading process. This conclusion was demonstrated in this study. We established a function p(t) to describe the cumulative effect of memory, and this function as a parameter was added to the rumor spreading model, and a dynamic analysis of the model was conducted.The simulations showed how the memory effects do affect rumor spreading when the parameters of the memory functionp(t) vary in regular networks, ER random networks, and BA networks. The simulation results indicated that a rumor spreads faster and more broadly in three network types when the initial memory rate p is larger; the larger the maximum value of the memory rate q, the faster the memory speed c, and the greater the peak value for the spreaders. These results suggest that rumor spreading is sensitive to parameter c in random networks, BA networks are most effective for rumor spreading, and the initial memory rate has little influence on the spread of rumor. We conducted a sensitivity analysis on the average degree in our model, the results of which demonstrated that the average degree plays a positive role in promoting rumor spreading. The final number of people affected by the rumors was found to be sensitive to the average degree k when k is small. Finally, the numerical solutions also demonstrated that the spreader peak value and final rumor spreading size on the Sina Microblog were much higher when the memory effects rate changes over time, compared to the case when the memory effects rate is considered a constant. Further study will be conducted which includes the other differences between rumor spreading and epidemic spreading. Further, investigation of rumor control strategies in different network topologies is also a future research direction.



---

*Source: 204395-2015-01-14.xml*