Steady-State Analysis of SECIR Rumor Spreading Model in Complex Networks

Yujiang Liu, Chunmei Zeng, Youquan Luo

Key Laboratory of Jiangxi Province for Numerical Simulation and Emulation Techniques, Gannan Normal University, Ganzhou, China.

DOI: 10.4236/am.2019.103007                

Abstract

In this paper, the SECIR rumor spreading model is formulated and analyzed, in which the social education level and the counterattack mechanism are taken into consideration. The results show that improving education level and increasing the ratio of counter are effective in reducing the risk of rumor propagation and enhancing the resistance to rumor propagation.

Keywords

SECIR Model, Rumor Spreading, Homogeneous Network, Counterattack Mechanism

1. Introduction

Rumor is a kind of social phenomenon that an unverified account or explanation of events spreads on a large-scale in a short time through people’s communication [1] [2] [3] . The spread of rumor can manipulate the public opinion in a locality, even can cause panic in some of the important public event [4] [5] [6] [7] [8] . Today, the increasing prevalence of social networking services, rumors spread by twitters, blogs, microblogs, WeChat and so on. In the Internet, the spreading of rumor is similar to epidemic spreading, but rumor’s spreading quantitative models have been rather limited in the complex network.

The standard model of rumor spreading is the Daley-Kendal (DK) model [9] [10] . The population in a local area is grouped ignorants, spreaders and stiflers. In this model, the rumor is propagated through the pair-wise contacts between spreaders and the others in the population. Spreader attempts to “infect” the other individuals with rumor. In the contact, ignorance becomes spreader, the others of individuals become stiflers. In the Maki-Thompson (MK) model [11] as the DK model’s variant, the spreader becomes stifler only who is the initiating spreader. The deficiencies of the DK model are considering homogeneous topology, and the simplified topology may not adequately describe the rumor’s spreading process in the Internet [12] - [17] . Zanette [18] [19] analyzed the MK model on a small-world network. His studies show that rumor “dies” in a small scale of its origin with varying network randomness. Morno et al. [20] [21] considered a rumor spreading model on a scale-free network. The results are the uniformity of the network which has a great impact on the rumor’s spreading process. Nekovee et al. [22] introduced the SIR model with forgetting mechanism and derived mean-field equation that describes the dynamics of rumor spreading process. His studies show that the SIR model is suitable for chain emails and large-scale information dissemination algorithms on the Internet. Zhao et al. [23] [24] refined the SIR rumor spreading model and took into account the remembering mechanism in addition to the forgetting mechanism. Zan et al. [25] considered the counterattack mechanism in the SIR model and introduced two models: Susceptible-infective-counterattack-Refractory (SICR) model and adjusted-SICR model. They derived mean-field equations to describe its dynamics in homogeneous network and involve steady-state analysis. Their studies show the self-resistance characteristic of networks to a rumor. Afassinou [26] extended the SIR model with the forgetting mechanism and population’s education rate and introduced SEIR model. He distinguishes two types of individuals in a population: educated individuals and non-educated individuals. His results show that improving the education rate of the population catalyzes the rumor spreading termination process. In social networks, when people with a higher degree of education heard a rumor which is in serious conflict with his/her belief, he/she is easier to counterattack the rumor, and even do the best to prevent the rumor propagation.

In this paper, inspiring of Zan et al. [25] , we consider two influential factors in rumor spreading process: the population’s education rate and the self-resistance feature of network. Motivatedly, we extend rumor spreading model―the SECIR model.

The remaining part of the paper is organized as follows. We formulate the propagation mechanism of the SECIR model in a social network, and derive a system of nonlinear ordinary differential equations that describe dynamics of rumor spreading process in Section 2. In Section 3, we analyze the steady-state of the SECIR model. We give some of the conclusions in the last section.

2. Model

We study the SECIR model in a closed homogeneously mixed population that we differentiate into five distinct classes: the rumor-mongers (spreader, S), those who are spreading the rumor, the non-educated ignorants individuals class (Ignorant, I), the people who never heard the rumor, the educated ignorants individuals class (Educatee, E), the people who never heard the rumor, but have more sophisticated behaviours with the non-educated ignorant individuals when they encountered the spreader, the counterattack class (Counter, C), those who do not agree but refute the rumor, and persuade others to agree with him (refute the rumor), and the stiflers class (Recovered, R), the ones who heard the rumor but have lost interest in disseminating it. For simplicity, we refer to the rumor spreading model as the SECIR model.

According to the MK model, we assume the rumor spreads by directed contact of the spreads with others in the population, and the contacts between rumor-mongers and the rest of the population are governed by the following dynamics (As shown in Figure 1):

・ Whenever a non-educated ignorant contact a spreader, the ignorant will be a relatively large probability ( ${\beta }_1$ ) into the spreader, but also with a smaller probability ( ${\alpha }_1$ ) of transition to a stifler;

・ Whenever an educated ignorant contact a spreader, the ignorant will be one of the three class: spreader or stifler analog the non-educated ignorant, and he can be evolved into counterattack, and the change probability is ${\beta }_2$ , ${\alpha }_2$ and $\theta$ , respectively;

・ Whenever a spreader contact a spreader, a stifler or a counter, the spreader will be into stifler with the probability g, g and $\eta$ , respectively;

Let $I\left(t\right)$ , $E\left(t\right)$ , $C\left(t\right)$ , $C\left(t\right)$ , and $R\left(t\right)$ respectively represent the density of the corresponding compartment in the total population. Namely that we have $I\left(t\right)+E\left(t\right)+C\left(t\right)+S\left(t\right)+R\left(t\right)=1$ . Note that the resistance to the rumor of the educated ignorants more than the non-educated ignorants, we assume that

${\beta }_1>{\beta }_2,{\alpha }_1>{\alpha }_2\text{and}\eta >g.$ (1)

In accordance with the above rules, the mean-field equations of the SECIR model can be described as follows:

$\begin{array}{l}\frac{\text{d}I\left(t\right)}{\text{d}t}=-\left({\beta }_1+{\alpha }_1\right)\bar{k}I\left(t\right)S\left(t\right)\\ \frac{\text{d}E\left(t\right)}{\text{d}t}=-\left({\beta }_2+{\alpha }_2+\theta \right)\bar{k}E\left(t\right)S\left(t\right)\\ \frac{\text{d}C\left(t\right)}{\text{d}t}=\theta \bar{k}E\left(t\right)S\left(t\right)\\ \frac{\text{d}S\left(t\right)}{\text{d}t}={\beta }_1\bar{k}I\left(t\right)S\left(t\right)+{\beta }_2\bar{k}E\left(t\right)S\left(t\right)-g\bar{k}S\left(t\right)\left(S\left(t\right)+R\left(t\right)\right)-\eta \bar{k}S\left(t\right)C\left(t\right)\\ \frac{\text{d}R\left(t\right)}{\text{d}t}={\alpha }_1\bar{k}I\left(t\right)S\left(t\right)+{\alpha }_2\bar{k}E\left(t\right)S\left(t\right)+g\bar{k}S\left(t\right)\left(S\left(t\right)+R\left(t\right)\right)+\eta \bar{k}S\left(t\right)C\left(t\right)\end{array}$ (2)

Figure 1. SECIR rumor spreading model.

We assume that all in the population are I or E but only one spreader at the beginning of the rumor spreading, and the ratio of E to the sum of I and E is $ϵ$ when $t=0$ . Namely, when $t=0$ , the initial condition of rumor spreading is given as follows:

$I\left(0\right)=\frac{N-1}{N}\left(1-ϵ\right),E\left(0\right)=\frac{N-1}{N}ϵ,S\left(0\right)=\frac{1}{N},C\left(0\right)=0,R\left(0\right)=0.$ (3)

Note that for an ignorant, he/she can be a spreader or stifler, so we have ${\beta }_1+{\alpha }_1<1$ , and $1-{\beta }_1-{\alpha }_1$ is the probability that no one tells him/her the rumor. For the same reason, ${\beta }_2+{\alpha }_2+\theta <1$ is considered.

3. Model Analysis

We then postulate that the number of individuals in the I class that have heard about the rumor is the same as the number of individuals in the E class that have heard about the rumor. This simply translates to

$\begin{array}{c}{\beta }_1+{\alpha }_1={\beta }_2+{\alpha }_2+\theta \\ =M\end{array}$ (4)

From the first and the second equation of (2), we have

$\frac{\text{d}E\left(t\right)}{\text{d}I\left(t\right)}=\frac{E\left(t\right)}{I\left(t\right)}$ (5)

Solve the above differential equation with the initial conditions (3), we obtain

$\begin{array}{c}\frac{E\left(t\right)}{I\left(t\right)}=\frac{E\left(0\right)}{I\left(0\right)}\\ =\frac{ϵ}{1-ϵ}\end{array}$ (6)

From the first and the third Equation (2), we obtain

$\frac{\text{d}C\left(t\right)}{\text{d}t}=-\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\frac{\text{d}E\left(t\right)}{\text{d}t}$

With the initial conditions (3), we can derive the relational expression between $C\left(t\right)$ and $E\left(t\right)$ by separation of variable,

$C\left(t\right)=\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\left(ϵ-E\left(t\right)\right)$ (7)

From (6) and (7), we have

$\begin{array}{c}S\left(t\right)+R\left(t\right)=1-\left(I\left(t\right)+C\left(t\right)+E\left(t\right)\right)\\ =1-\frac{1-ϵ}{ϵ}E\left(t\right)-\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\left(ϵ-E\left(t\right)\right)-E\left(t\right)\\ =\left(1-\frac{\theta ϵ}{{\alpha }_2+{\beta }_2+\theta }\right)-\left(\frac{1}{ϵ}-\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\right)E\left(t\right)\end{array}$ (8)

Note that $\tau =\frac{g}{ϵ}+\frac{\left(\eta -g\right)\theta }{{\alpha }_2+{\beta }_2+\theta }$ and $\beta =\frac{1-ϵ}{ϵ}{\beta }_1+{\beta }_2$ , then we have

$\begin{array}{c}\frac{\text{d}S\left(t\right)}{\text{d}t}={\beta }_1\bar{k}I\left(t\right)S\left(t\right)+{\beta }_2\bar{k}E\left(t\right)S\left(t\right)-g\bar{k}S\left(t\right)\left(S\left(t\right)+R\left(t\right)\right)-\eta \bar{k}S\left(t\right)C\left(t\right)\\ ={\beta }_1\bar{k}\frac{1-ϵ}{ϵ}E\left(t\right)S\left(t\right)+{\beta }_2\bar{k}E\left(t\right)S\left(t\right)-\eta \bar{k}S\left(t\right)\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\left(ϵ-E\left(t\right)\right)\\ -g\bar{k}S\left(t\right)\left(\left(1-\frac{\theta ϵ}{{\alpha }_2+{\beta }_2+\theta }\right)-\left(\frac{1}{ϵ}-\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\right)E\left(t\right)\right)\\ =\left(\frac{1-ϵ}{ϵ}{\beta }_1+{\beta }_2+\frac{g}{ϵ}+\frac{\left(\eta -g\right)\theta }{{\alpha }_2+{\beta }_2+\theta }\right)\bar{k}E\left(t\right)S\left(t\right)-\left(g+\frac{\theta \left(\eta -g\right)ϵ}{{\alpha }_2+{\beta }_2+\theta }\right)\bar{k}S\left(t\right)\\ =\left(\beta +\tau \right)\bar{k}E\left(t\right)S\left(t\right)-ϵ\tau \bar{k}S\left(t\right)\end{array}$ (9)

From the second equation of (2), therefore (9) becomes

$\frac{\text{d}S\left(t\right)}{\text{d}t}=\frac{\frac{ϵ\tau }{E\left(t\right)}-\left(\beta +\tau \right)}{M}\frac{\text{d}E\left(t\right)}{\text{d}t}$

Solving the differential equations above by the method of separation of variables, we have

$S\left(t\right)=\frac{ϵ\tau }{M}\left(\ln E\left(t\right)-\ln ϵ\right)-\frac{\beta +\tau }{M}\left(E\left(t\right)-ϵ\right)$ (10)

By the second equation of (2), it is easy to see that $\frac{\text{d}E\left(t\right)}{\text{d}t}<0$ and $E\left(t\right)$ are monotonically decreasing and continuous function. Let $\frac{\text{d}S\left(t\right)}{\text{d}t}=0$ , we obtain $E\left(t\right)=\frac{ϵ\tau }{\beta +\tau }$ . It is easy to see that, the peak value of spreader is

$S_{\max }=\frac{ϵ\beta }{M}+\frac{ϵ\tau }{M}\ln \left(\frac{\tau }{\beta +\tau }\right)$ (11)

$\begin{array}{c}\frac{\text{d}R\left(t\right)}{\text{d}t}={\alpha }_1\bar{k}I\left(t\right)S\left(t\right)+{\alpha }_2\bar{k}E\left(t\right)S\left(t\right)+g\bar{k}S\left(t\right)\left(S\left(t\right)+R\left(t\right)\right)+\eta \bar{k}S\left(t\right)C\left(t\right)\\ ={\alpha }_1\bar{k}\frac{1-ϵ}{ϵ}E\left(t\right)S\left(t\right)+{\alpha }_2\bar{k}E\left(t\right)S\left(t\right)+\eta \bar{k}S\left(t\right)\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\left(ϵ-E\left(t\right)\right)\\ +g\bar{k}S\left(t\right)\left(\left(1-\frac{\theta ϵ}{{\alpha }_2+{\beta }_2+\theta }\right)-\left(\frac{1}{ϵ}-\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\right)E\left(t\right)\right)\\ =\left(\frac{1-ϵ}{ϵ}{\alpha }_1+{\alpha }_2-\frac{g}{ϵ}-\frac{\left(\eta -g\right)\theta }{{\alpha }_2+{\beta }_2+\theta }\right)\bar{k}E\left(t\right)S\left(t\right)+\left(g+\frac{\theta \left(\eta -g\right)ϵ}{{\alpha }_2+{\beta }_2+\theta }\right)\bar{k}S\left(t\right)\\ =\left(\alpha -\tau \right)\bar{k}E\left(t\right)S\left(t\right)+ϵ\tau \bar{k}S\left(t\right)\end{array}$ (12)

where $\alpha =\frac{1-ϵ}{ϵ}{\alpha }_1+{\alpha }_2$ . From the second equation of (2), therefore (12) becomes

$\frac{\text{d}R\left(t\right)}{\text{d}t}=\left(\frac{\tau -\alpha }{{\alpha }_2+{\beta }_2+\theta }-\frac{ϵ\tau }{{\alpha }_2+{\beta }_2+\theta }\frac{1}{E\left(t\right)}\right)\frac{\text{d}E\left(t\right)}{\text{d}t}$ (13)

Solving the above differential equations, we get

$R\left(t\right)=\frac{\alpha -\tau }{{\alpha }_2+{\beta }_2+\theta }\left(ϵ-E\left(t\right)\right)+\frac{ϵ\tau }{{\alpha }_2+{\beta }_2+\theta }\left(\ln ϵ-\ln E\left(t\right)\right)$ (14)

It is easy know that $S_{\infty }=0$ , let $t\to \infty$ , then (8) becomes

$R_{\infty }=\left(\frac{1}{ϵ}-\frac{\theta }{{\alpha }_2+{\beta }_2+\theta }\right)\left(ϵ-E_{\infty }\right)$ (15)

Let $t\to \infty$ , Substituting (15) into (14), it becomes

$\left(\frac{1}{ϵ}-\frac{\alpha -\tau +\theta }{{\alpha }_2+{\beta }_2+\theta }\right)\left(ϵ-E_{\infty }\right)=\frac{ϵ\tau }{{\alpha }_2+{\beta }_2+\theta }\left(\ln ϵ-\ln E_{\infty }\right)$ (16)

$\begin{array}{c}\alpha +\beta =\frac{1-ϵ}{ϵ}\left({\beta }_1+{\alpha }_1\right)+\left({\beta }_2+{\alpha }_2\right)\\ =\frac{1-ϵ}{ϵ}\left({\beta }_2+{\alpha }_2\right)+\left({\beta }_2+{\alpha }_2\right)\\ =\frac{1}{ϵ}\left({\beta }_2+{\alpha }_2+\theta \right)-\theta \end{array}$ (17)

Solve from (17), we get $ϵ=\frac{{\beta }_2+{\alpha }_2+\theta }{\beta +\alpha +\theta }$ , substitute into (16), we have

$\left(ϵ-E_{\infty }\right)=\frac{ϵ\tau }{\beta +\tau }\left(\ln ϵ-\ln E_{\infty }\right)$ (18)

Let $A\left(t\right)=I\left(t\right)+E\left(t\right)=\frac{1-ϵ}{ϵ}E\left(t\right)+E\left(t\right)=E\left(t\right)/ϵ$ , from (18) we get the final size

$A_{\infty }=\frac{\tau }{\beta +\tau }\ln A_{\infty }+1$ (19)

Theorem 1. For $0<\sigma <1$ , the equation $x=\sigma \ln x+1$ has two solutions, $x=1$ and a nontrivial solution $x_1$ , where $0<x_1<\sigma$ .

Proof. Obviously $x=1$ is a solution of $x=\sigma \ln x+1$ .

Let $f\left(x\right)=x-\sigma \ln x-1$ , and take the derivative of $f\left(x\right)$ with respect to x: $f^{\prime }\left(x\right)=1-\sigma /x$ , $f^{″}\left(x\right)=\sigma /x^2>0$ .

Let $f^{\prime }\left(x\right)=0$ , we obtain the unique minimum point $x=\sigma$ , and the function $f\left(x\right)$ is a convex function, we have $f\left(\sigma \right)=\sigma -\sigma \ln \sigma -1<\sigma +\sigma \left(1/\sigma -1\right)-1=0$ , and $f\left(0^{+}\right)=\infty$ . According to the Mean Value Theorem, $f\left(x\right)$ have a nontrivial solution $x_1$ , where $0<x_1<\sigma$ . ,

Theorem 2. If the parameters are satisfied (1), (3) and (4), we have

(1) $S_{\max }$ is decrease with $ϵ$ , if the other parameters keep constant.

(2) $S_{\max }$ is increase with ${\beta }_1\mathrm{,}{\beta }_2$ but decrease with ${\alpha }_1\mathrm{,}{\alpha }_2$ , if the other parameters keep constant.

(3) $S_{\max }$ is decrease with $g\mathrm{,}\eta$ , if the other parameters keep constant.

Proof. (1) From (11), differential with $ϵ$ , we have

$\begin{array}{c}\frac{\text{d}{S}_{\max }}{\text{d}ϵ}=-\frac{{\beta }_1-{\beta }_2}{M}+\frac{ϵ\tau }{M}{\left(\ln \frac{\tau }{\beta +\tau }\right)}_{ϵ}+\frac{{\left(ϵ\tau \right)}^{\prime }{}_{ϵ}}{M}\ln \frac{\tau }{\beta +\tau }\\ =-\frac{{\beta }_1-{\beta }_2}{M}+\frac{ϵ\tau }{M}\left(\frac{-\beta g+\tau {\beta }_1}{Mϵ\left(\beta +\tau \right)}\right)+\frac{\left(\eta -g\right)\theta }{M^2}\ln \frac{\tau }{\beta +\tau }\end{array}$

$\begin{array}{c}\le -\frac{{\beta }_1-{\beta }_2}{M}+\frac{ϵ\tau }{M}\left(\frac{-\beta g+\tau {\beta }_1}{Mϵ\left(\beta +\tau \right)}\right)+\frac{\left(\eta -g\right)\theta }{M^2}\left(\frac{\tau }{\beta +\tau }-1\right)\\ =-\frac{\left(\left(1-ϵ\right){\beta }_1+{\beta }_2ϵ\right)\left({\beta }_1-{\beta }_2\right)}{M\left(\beta +\tau \right)ϵ}\end{array}$ (20)

Therefore, $\frac{\text{d}{S}_{\max }}{\text{d}ϵ}<0$ , which implies that $S_{\max }$ decreases as $ϵ$ increases.

(2) Since that $\tau =\frac{g}{ϵ}+\frac{\left(\eta -g\right)\theta }{{\alpha }_2+{\beta }_2+\theta }$ and $\beta =\frac{1-ϵ}{ϵ}{\beta }_1+{\beta }_2$ , we have ${{\tau }^{\prime }}_{{\beta }_i}=0$ , ${{\tau }^{\prime }}_{{\alpha }_i}=0$ , $i=1,2$ , ${{\beta }^{\prime }}_{{\beta }_i}>0$ , ${{\beta }^{\prime }}_{{\alpha }_i}<0$ . Taking the derivative of (11), with respect to ${\beta }_i$ , we have $\frac{\text{d}{S}_{\max }}{\text{d}{\beta }_i}=\frac{ϵ}{M}{{\beta }^{\prime }}_{{\beta }_i}-\frac{ϵ\tau }{M}\frac{{{\beta }^{\prime }}_{{\beta }_i}}{\beta +\tau }=\frac{ϵ}{M}\frac{\beta }{\beta +ϵ}{{\beta }^{\prime }}_{{\beta }_i}>0$ . It is easy to see that $S_{\max }$ increases as ${\beta }_i\left(i=1,2\right)$ increases.

By the same way, we can proof that $S_{\max }$ is decrease with ${\alpha }_1\mathrm{,}{\alpha }_2$ .

(3) Similarly (2), we have ${{\beta }^{\prime }}_{\eta }=0$ , ${{\beta }^{\prime }}_g=0$ , ${{\tau }^{\prime }}_x>0,\left(x=\eta \text{or}g\right)$ , then

$\begin{array}{c}\frac{\text{d}{S}_{\max }}{\text{d}x}=\frac{ϵ}{M}\ln \frac{\tau }{\beta +\tau }{{\tau }^{\prime }}_x+\frac{ϵ\tau }{M}\left(\frac{{{\tau }^{\prime }}_x}{\tau }-\frac{{{\tau }^{\prime }}_x}{\beta +\tau }\right)\\ =\left(\frac{ϵ}{M}\ln \frac{\tau }{\beta +\tau }+\frac{ϵ\beta }{M\left(\beta +\tau \right)}\right){{\tau }^{\prime }}_x\\ <\left(\frac{ϵ}{M}\left(\frac{\tau }{\beta +\tau }-1\right)+\frac{ϵ\beta }{M\left(\beta +\tau \right)}\right){{\tau }^{\prime }}_x\\ =0\end{array}$ (21)

which means that $S_{\max }$ is decrease with $g\mathrm{,}\eta$ The proof is complete. ,

Theorem 3. If the parameters are satisfied (1), (3) and (4), we have

(1) the other parameters keep constant, the final state $A_{\infty }$ is increased with $ϵ$ .

(2) the other parameters keep constant, the final state $A_{\infty }$ is decreased with ${\beta }_1\mathrm{,}{\beta }_2$ but decrease with ${\alpha }_1\mathrm{,}{\alpha }_2$ .

(3) the other parameters keep constant, the final state $A_{\infty }$ is increased with $g\mathrm{,}\eta$ .

Proof. Let $\sigma =\frac{\tau }{\beta +\tau }$ , from (19) we have

$\frac{\text{d}{A}_{\infty }}{\text{d}\chi }=\frac{\text{d}\sigma }{\text{d}\chi }\ln A_{\infty }+\sigma \frac{1}{A_{\infty }}\frac{\text{d}{A}_{\infty }}{\text{d}\chi }=\frac{A_{\infty }\ln A_{\infty }}{A_{\infty }-\sigma }\frac{\text{d}\sigma }{\text{d}\chi }$ (22)

From Theorem 1, $0<A_{\infty }<\sigma$ , it is obviously that $\frac{A_{\infty }\ln A_{\infty }}{A_{\infty }-\sigma }>0$ , so, if we know that the plus-minus sign of $\frac{\text{d}\sigma }{\text{d}\chi }$ , then the plus-minus sign of $\frac{\text{d}{A}_{\infty }}{\text{d}\chi }$ can be determined, and furthermore, the monotonicity between $A_{\infty }$ and parameters $\chi$ can be proved.

(1)

$\frac{\text{d}\sigma }{\text{d}ϵ}=\frac{{{\tau }^{\prime }}_{ϵ}\beta -\tau {{\beta }^{\prime }}_{ϵ}}{{\left(\beta +\tau \right)}^2}=\frac{{\beta }_1\tau -g\beta }{{ϵ}^2{\left(\beta +\tau \right)}^2}=\frac{\left({\beta }_1-{\beta }_2\right)g+\frac{\left(\eta -g\right)\theta }{M}{\beta }_1}{{ϵ}^2{\left(\beta +\tau \right)}^2}$ (23)

so $\frac{\text{d}\sigma }{\text{d}ϵ}>0$ , that is, $\sigma$ increases as $ϵ$ increases, and also implies that $A_{\infty }$ increases as $ϵ$ increases.

(2)

$\frac{\text{d}\sigma }{\text{d}{\beta }_i}=-\frac{\tau {{\beta }^{\prime }}_{{\beta }_i}}{{\left(\beta +\tau \right)}^2}<0,\left(i=1,2\right)$ (24)

similarly, we can get $\frac{\text{d}\sigma }{\text{d}{\alpha }_i}>0$ . and so $\frac{A_{\infty }}{\text{d}{\beta }_i}<0$ , $\frac{A_{\infty }}{\text{d}{\alpha }_i}<0$ , where $i=1,2$ , which proves (2).

(3)

$\frac{\text{d}\sigma }{\text{d}x}=\frac{{{\tau }^{\prime }}_x\beta }{{\left(\beta +\tau \right)}^2},x\in \left\{\eta ,g\right\}$ (25)

then, $\frac{\text{d}\sigma }{\text{d}x}>0$ , it means that $A_{\infty }$ increases as $\eta \left(g\right)$ increases. The proof is complete.

4. Numerical Simulation

We assume $N={10}^5$ , the average degree of network $\bar{k}=10$ , and the initial condition of the model follows equation of (2).

Figure 2 shows the general trends of the five kinds of agents in the SEICR rumor spreading model. We see the density of spreaders begin to expand rapidly from the initial rumor spread. As the rumor spread further, the density of spreaders reaches a peak and thereafter declines. Finally the density of spreaders is to zero and this leads to the termination of rumor spreading. And over the course of the rumor spreading, the density of ignorants and educatees always decreases, and finally evolves to zero. but the density of counter and the recovered

Figure 2. Density of spreaders, educatees, ignorants, counters and recovereds over time, with ${\alpha }_1={\beta }_1=0.5,{\alpha }_2={\beta }_2=0.4,g=0.2,\eta =0.1,ϵ=0.5$ .

always decreases, and eventually evolve to a stable value.

Figures 3-5 show how the densities of spreaders change over time for different system parameters include $ϵ\mathrm{,}{\alpha }_1\mathrm{,}{\beta }_1\mathrm{,}{\alpha }_2\mathrm{,}{\beta }_2\mathrm{,}g$ and $\eta$ , and the change is consistent with theorem 2. It is interest that the higher parameter ${\beta }_1\mathrm{,}{\beta }_2$ , the earlier the outbreak, the larger the peak of the outbreak, but the shorter the outbreak period (Figure 4). However, if the parameter X is lowered, the outbreak period will not come earlier, and the outbreak period will be longer (in Figure 5).

5. Conclusion

In this paper, considering the social education level and the counterattack mechanism,

Figure 3. Density of spreaders over time under different population’s education rate $ϵ$ , with ${\alpha }_1={\beta }_1=0.5,{\alpha }_2={\beta }_2=0.4,g=0.2,\eta =0.1$ .

Figure 4. Density of spreaders over time under different ${\alpha }_1\mathrm{,}{\beta }_1\mathrm{,}{\alpha }_2\mathrm{,}{\beta }_2$ , with $g=0.2,\eta =0.1,\epsilon =0.5$ .

Figure 5. Density of spreaders over time under different $g\mathrm{,}\eta$ , with ${\alpha }_1={\beta }_1=0.5,{\alpha }_2={\beta }_2=0.4,g=0.2,\eta =0.1,ϵ=0.5$ .

we analyze the dynamics of rumor propagation. The results of simulations show that improving education level and increasing the ratio of counter are effective in reducing the risk of rumor propagation and enhancing the resistance to rumor propagation.

Acknowledgements

The research has been supported by The National Natural Science Foundation of China (11561004) and The 12th Five-year Education Scientific Planning Project of Jiangxi Province (15ZD3LYB031).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Peterson, W.A. and Gist, N.P. (1951) Rumor and Public Opinion. American Journal of Sociology, 57, 159-167. https://doi.org/10.1086/220916
[2]	Kesten, H. and Sidoravicius, V. (2005) The Spread of a Rumor or Infection in a Moving Population. Annals of Probability, 33, 2402-2462. https://doi.org/10.1214/009117905000000413
[3]	Kostka, J., Oswald, Y.A. and Wattenhofer, R. (2008) Word of Mouth: Rumor Dissemination in Social Networks. International Colloquium on Structural Information & Communication Complexity, Villars-sur-Ollon, 17-20 June 2008, 185-196. https://doi.org/10.1007/978-3-540-69355-0_16
[4]	Ganesh, A.J., Kermarrec, A.M. and Massoulie, L. (2002) Hiscamp: Self-Organizing Hierarchical Membership Protocol. Workshop on ACM Sigops European Workshop. https://doi.org/10.1145/1133373.1133398
[5]	Thomas, S.A. (2007) Lies, Damn Lies, and Rumors: An Analysis of Collective Efficacy, Rumors, and Fear in the Wake of Katrina. Sociological Spectrum, 27, 679-703. https://doi.org/10.1080/02732170701534200
[6]	Bhavnani, R. and Kuklinski, J.H. (2009) Rumor Dynamics in Ethnic Violence. Journal of Politics, 71, 876. https://doi.org/10.1017/S002238160909077X
[7]	Galam, S. (2002) Modelling Rumors: The No Plane Pentagon French Hoax Case. Physica A: Statistical Mechanics & Its Applications, 320, 571-580.
[8]	Kimmel, A.J. (2004) Rumors and Rumor Control: A Manager’s Guide to Understanding and Combatting Rumors. Lawrence Erlbaum Associates.
[9]	Daley, D.J. and Kendall, D.G. (1965) Stochastic Rumours. IMA Journal of Applied Mathematics, 1, 42-55. https://doi.org/10.1093/imamat/1.1.42
[10]	Daley, D.J. and Gani, J. (1999) Epidemic Modelling: An Introduction. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511608834
[11]	Maki, D.P. and Thompson, M. (1973) Mathematical Models and Applications: With Emphasis on the Social, Life, and Management Sciences. Prentice Hall, Upper Saddle River.
[12]	Pittel, B. (1990) On a Daley-Kendall Model of Random Rumours. Journal of Applied Probability, 27, 14-27. https://doi.org/10.2307/3214592
[13]	Newman, M.E., Forrest, S. and Balthrop, J. (2002) Email Networks and the Spread of Computer Viruses. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 66, Article ID: 035101. https://doi.org/10.1103/PhysRevE.66.035101
[14]	Sudbury, A. (1985) The Proportion of the Population Never Hearing a Rumour. Journal of Applied Probability, 22, 443-446. https://doi.org/10.2307/3213787
[15]	Holger, E., Lutz-Ingo, M. and Stefan, B. (2002) Scale-Free Topology of E-Mail Networks. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 66, Article ID: 035103.
[16]	Smith, R.D. (2002) Instant Messaging as a Scale-Free Network.
[17]	Fang, W., Yamir, M. and Yaoru, S. (2006) Structure of Peer-to-Peer Social Networks. Physical Review E Statistical Nonlinear & Soft Matter Physics, 73, Article ID: 036123.
[18]	Zanette, D.H. (2001) Critical Behavior of Propagation on Small-World Networks. Physical Review E Statistical Nonlinear & Soft Matter Physics, 64, Article ID: 050901. https://doi.org/10.1103/PhysRevE.64.050901
[19]	Zanette, D.H. (2002) Dynamics of Rumor Propagation on Small-World Networks. Physical Review E Statistical Nonlinear & Soft Matter Physics, 65, Article ID: 041908. https://doi.org/10.1103/PhysRevE.65.041908
[20]	Yamir, M., Maziar, N. and Alessandro, V. (2004) Efficiency and Reliability of Epidemic Data Dissemination in Complex Networks. Physical Review E Statistical Nonlinear & Soft Matter Physics, 69, Article ID: 055101.
[21]	Moreno, Y., Pastorsatorras, R. and Vespignani, R. (2001) Epidemic Outbreaks in Complex Heterogeneous Networks. The European Physical Journal B Condensed Matter and Complex Systems, 26, 521-529. https://doi.org/10.1140/epjb/e20020122
[22]	Nekovee, M., Moreno, Y., Bianconi, G. and Marsili, M. (2008) Theory of Rumour Spreading in Complex Social Networks. Physica A Statistical Mechanics & Its Applications, 374, 457-470. https://doi.org/10.1016/j.physa.2006.07.017
[23]	Zhao, L., Cui, H., Qiu, X. and Wang, X. (2013) Sir Rumor Spreading Model in the New Media Age. Physica A Statistical Mechanics & Its Applications, 392, 995-1003.
[24]	Zhao, L., Wang, J., Chen, Y., Qin, W., Cheng, J. and Cui, H. (2012) Sihr Rumor Spreading Model in Social Networks. Physica A Statistical Mechanics & Its Applications, 391, 2444-2453. https://doi.org/10.1016/j.physa.2011.12.008
[25]	Zan, Y., Wu, J., Ping, L. and Yu, Q. (2014) Sicr Rumor Spreading Model in Complex Networks: Counterattack and Self-Resistance. Physica A Statistical Mechanics & Its Applications, 405, 159-170. https://doi.org/10.1016/j.physa.2014.03.021
[26]	Afassinou, K. (2014) Analysis of the Impact of Education Rate on the Rumor Spreading Mechanism. Physica A Statistical Mechanics & Its Applications, 414, 43-52. https://doi.org/10.1016/j.physa.2014.07.041