Dynamics of a New Rumor Propagation Model with the Spread of Truth

Abstract

A mathematical model described the propagation of information including rumor and truth presented and its properties investigated. We explored exists of the equilibria, local stability and global asymptotical stability, and obtained the propagation threshold of rumor spreading. Numerical simulation is shown to demonstrate our results. Uncertainty and sensitivity analysis shows the importance of the parameters in our model.

Keywords

Share and Cite:

Liu, Y. , Zeng, C. and Luo, Y. (2018) Dynamics of a New Rumor Propagation Model with the Spread of Truth. Applied Mathematics, 9, 536-549. doi: 10.4236/am.2018.95038.

1. Introduction

Rumors, as daily events have been occurred. With the development of new media such as Internet, mobile phone has broadened the public access to information, optimized the circulation of information, but also exacerbated the spread of rumors. In emergency communication, financial markets fluctuation and all kinds of contagions, rumors play an important role [1] [2] [3] . Therefore, understanding the propagation of rumors and how to control it effectively is a very meaningful topic.

The remainder of this paper is organized as follows. In Section 2, we derived the model considering the influence of rumors, after which a detailed steady-state analysis is carried out in Section 3, and the threshold of the rumor propagation is obtained. In Section 4, numerical simulation, uncertainty and sensitivity analysis of the model are presented. Finally, Section 5 summarizes this work.

2. Model

We assume that there are two kinds of spreaders in the whole population which called spreader with rumor and spreader with truth, respectively. As shown in Figure 1, the total population is divided into four different compartments name ignorants, spreaders with rumor, spreaders with truth and recovers, noted with $I\left(t\right)$ , ${S}_{+}\left(t\right)$ , ${S}_{-}\left(t\right)$ , $R\left(t\right)$ respectively. When an ignorant contacts a spreader with truth, the ignorant becomes the other side with rate ${\beta }_{+}$ , otherwise, if he contacts a spread with rumor, the ignorant becomes the other party with rate ${\beta }_{-}$ . When a spreader with truth contacts with spreader with rumor, they will become a

Figure 1. The flow diagram of the rumor propagation model.

removed with probability ${\lambda }_{+}$ and ${\lambda }_{-}$ . The spreaders ( ${S}_{+}\{S}_{-}$ ) spontaneously become the removed at a rate ( ${\gamma }_{+}\{\gamma }_{-}$ ) for losing interesting.

According to the dynamic interact as mentioned above, we established the S2IR rumor spreading model based on the above assumptions. The model is described as follows:

$\begin{array}{l}\stackrel{˙}{I}=\Lambda -{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I-{\mu }_{1}I\\ {\stackrel{˙}{S}}_{+}={\beta }_{+}{S}_{+}I-{\lambda }_{+}{S}_{-}{S}_{+}-{\gamma }_{+}{S}_{+}-{\mu }_{2}{S}_{+}\\ {\stackrel{˙}{S}}_{-}={\beta }_{-}{S}_{-}I-{\lambda }_{-}{S}_{-}{S}_{+}-{\gamma }_{-}{S}_{-}-{\mu }_{3}{S}_{-}\\ \stackrel{˙}{R}={\lambda }_{+}{S}_{-}{S}_{+}+{\lambda }_{-}{S}_{-}{S}_{+}-{\gamma }_{-}+{\gamma }_{+}{S}_{+}+{\gamma }_{-}{S}_{-}-{\mu }_{4}R\end{array}$ (1)

The parameters of the model descript as follows:

$\Lambda$ , the constant recruitment rate of the population;

${\mu }_{1}$ , ${\mu }_{2}$ , ${\mu }_{3}$ , the natural fade away rate to recover of the ignorants, spreaders with rumor and spreaders with truth respectively;

${\beta }_{+}$ , ${\beta }_{-}$ , the force of infection from ignorants to spreaders of rumor, spreaders of truth respectively;

${\lambda }_{+}$ , the fade away rate to recover when a spreader of rumor disseminate the rumor to a spreader of truth;

${\lambda }_{-}$ , the fade away rate to recover when a spreader of truth disseminate the truth to a spreader of rumor;

${\gamma }_{+}$ , ${\gamma }_{-}$ , the forgetting rate of spreader of rumor and spreader of truth, respectively.

The initial condition of (1) is given as $I\left(0\right)>0,{S}_{+}\left(0\right)\ge 0,{S}_{-}\left(0\right)\ge 0,R\left(0\right)\ge 0$ , for the fourth equation of system of (1) is not comprised in the first three equation. Then we can consider the simplified model as follows:

$\begin{array}{l}\stackrel{˙}{I}=\Lambda -{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I-{\mu }_{1}I\\ {\stackrel{˙}{S}}_{+}={\beta }_{+}{S}_{+}I-{\lambda }_{+}{S}_{-}{S}_{+}-{\gamma }_{+}{S}_{+}-{\mu }_{2}{S}_{+}\\ {\stackrel{˙}{S}}_{-}={\beta }_{-}{S}_{-}I-{\lambda }_{-}{S}_{-}{S}_{+}-{\gamma }_{-}{S}_{-}-{\mu }_{3}{S}_{-}\end{array}$ (2)

The initial condition of (2) is given as $I\left(0\right)>0,{S}_{+}\left(0\right)\ge 0,{S}_{-}\left(0\right)\ge 0$ . Note $\mu =max\left\{{\mu }_{1},{\mu }_{2},{\mu }_{3}\right\}$ Then we can easily verify that the system of (2) is mathematically well posed in the positive invariant region $G=\left\{\left(I,{S}_{+},{S}_{-}\right)|I>0,{S}_{+}\ge 0,{S}_{-}\ge 0,I+{S}_{+}+{S}_{-}\le \Lambda /\mu \right\}$ . And all the solutions with $I\left(0\right)>0,{S}_{+}\left(0\right)\ge 0,{S}_{-}\left(0\right)\ge 0$ approaches or enter the region G, so it is sufficient to consider solutions in G.

3. Equilibria and Stability

Noting that

${R}_{0}=\mathrm{max}\left\{{R}_{01},{R}_{02}\right\}$ (3)

where ${R}_{01}=\frac{\Lambda {\beta }_{+}}{{\mu }_{1}\left({\gamma }_{+}+{\mu }_{2}\right)}$ , ${R}_{02}=\frac{\Lambda {\beta }_{-}}{{\mu }_{1}\left({\gamma }_{-}+{\mu }_{3}\right)}$ , and

${R}^{*}=\mathrm{min}\left\{{R}_{1}^{*},{R}_{2}^{*}\right\}$ (4)

where ${R}_{1}^{*}=\frac{\Lambda {\beta }_{+}{\lambda }_{-}+{\beta }_{+}\left({\mu }_{2}+{\gamma }_{+}\right)\left({\mu }_{3}+{\gamma }_{-}\right)}{{\beta }_{-}{\left({\mu }_{2}+{\gamma }_{+}\right)}^{2}+{\lambda }_{-}{\mu }_{1}\left({\mu }_{2}+{\gamma }_{+}\right)}$ , ${R}_{2}^{*}=\frac{\Lambda {\beta }_{-}{\lambda }_{+}+{\beta }_{-}\left({\mu }_{3}+{\gamma }_{-}\right)\left({\mu }_{2}+{\gamma }_{+}\right)}{{\beta }_{+}{\left({\mu }_{3}+{\gamma }_{-}\right)}^{2}+{\lambda }_{+}{\mu }_{1}\left({\mu }_{3}+{\gamma }_{-}\right)}$ .

Letting the left-hand side of the differential equations of model (2) equal to zero yields the following equations

$\begin{array}{l}0=\Lambda -{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I-{\mu }_{1}I\\ 0={\beta }_{+}{S}_{+}I-{\lambda }_{+}{S}_{-}{S}_{+}-{\gamma }_{+}{S}_{+}-{\mu }_{2}{S}_{+}\\ 0={\beta }_{-}{S}_{-}I-{\lambda }_{-}{S}_{-}{S}_{+}-{\gamma }_{-}{S}_{-}-{\mu }_{3}{S}_{-}\end{array}$ (5)

Obviously, the system exists a trivial equilibrium ${E}_{0}\left({I}^{0},0,0\right)$ where ${I}^{0}=\Lambda /{\mu }_{1}$ . If ${R}_{01}>1$ , let ${S}_{-}=0$ , we easily get the rumor-free equilibrium ${E}_{1}\left({I}^{1},{S}_{+}^{1},0\right)$
where ${I}^{1}=\left({\gamma }_{+}+{\mu }_{2}\right)/{\beta }_{+}$ , and

${S}_{+}^{1}=\Lambda /\left({\gamma }_{+}+{\mu }_{2}\right)-{\mu }_{1}/{\beta }_{+}=\frac{{\mu }_{1}}{{\beta }_{+}}\left({R}_{01}-1\right)$ . Similarly, If ${R}_{02}>1$ , we get the counter-free equilibrium ${E}_{2}\left({I}^{2},0,{S}_{-}^{2}\right)$ where ${I}^{2}=\left({\gamma }_{-}+{\mu }_{3}\right)/{\beta }_{-}$ , and ${S}_{-}^{2}=\Lambda /\left({\gamma }_{-}+{\mu }_{3}\right)-{\mu }_{1}/{\beta }_{-}=\frac{{\mu }_{1}}{{\beta }_{-}}\left({R}_{02}-1\right)$ .

In the following, we shall study the existence of the positive equilibrium ${E}^{*}\left({I}^{*},{S}_{+}^{*},{S}_{-}^{*}\right)$ of system (1). From the second and the third equation of (5), we have ${S}_{+}^{*}=\left({\beta }_{-}{I}^{*}-\left({\mu }_{3}+{\gamma }_{-}\right)\right)/{\lambda }_{-}$ , ${S}_{-}^{*}=\left({\beta }_{+}{I}^{*}-\left({\mu }_{2}+{\gamma }_{+}\right)\right)/{\lambda }_{+}$ .

Substituting them in the first equation of (5), we have

$F\left(I\right)=A{\left({I}^{*}\right)}^{2}+B{I}^{*}-C$ (6)

where $A={\beta }_{+}{\beta }_{-}\left({\lambda }_{+}+{\lambda }_{-}\right)$ , $B=\left({\lambda }_{+}+{\lambda }_{-}\right){\mu }_{1}-{\beta }_{+}{\lambda }_{+}\left({\mu }_{3}+{\gamma }_{-}\right)-{\beta }_{-}{\lambda }_{-}\left({\mu }_{2}+{\gamma }_{+}\right)$ , and $C=\Lambda {\lambda }_{+}{\lambda }_{-}$ . Obviously, $A>0,C>0$ then exist a unique positive solution $I={I}^{*}$ , yield $F\left({I}^{*}\right)=0$ . For the ${S}_{+}^{*}>0$ and ${S}_{-}^{*}>0$ , we have $F\left(\left({\mu }_{2}+{\gamma }_{+}\right)/{\beta }_{+}\right)<0$ and $F\left(\left({\mu }_{3}+{\gamma }_{-}\right)/{\beta }_{-}\right)<0$

From (6) we get that ${R}_{1}^{*}>1$ , ${R}_{2}^{*}>1$ , namely ${R}^{*}>0$ . Thus, we obtain the following theorem.

Theorem 1. System (2) has the following equilibrium:

1) if ${R}_{0}<1$ , system (2) has only the trivial equilibrium ${E}_{0}$ ;

2) if ${R}_{01}>1$ , ${R}_{02}<1$ , system (2) has a trivial equilibrium ${E}_{0}$ and a truth-free equilibrium ${E}_{1}$ ;

3) if ${R}_{01}<1$ , ${R}_{02}>1$ , system (2) has a trivial equilibrium ${E}_{0}$ and a rumor-free equilibrium ${E}_{2}$ ;

4) if ${R}_{01}<1$ , ${R}_{02}>1$ , system (2) has a trivial equilibrium ${E}_{0}$ , a truth-free equilibrium ${E}_{1}$ and a rumor-free equilibrium ${E}_{2}$ ;

5) if ${R}_{01}>1$ , ${R}_{02}>1$ , ${R}^{*}>1$ , system (2) has a trivial equilibrium ${E}_{0}$ , a truth-free equilibrium ${E}_{1}$ , a rumor-free equilibrium ${E}_{2}$ and coexist equilibrium ${E}^{*}$ .

The general Jacobian of (2) is given by

${J}_{\left(2.2\right)}=|\begin{array}{ccc}-{\beta }_{+}{S}_{+}-{\beta }_{-}{S}_{-}-{\mu }_{1}& -{\beta }_{+}I& -{\beta }_{-}I\\ {\beta }_{+}{S}_{+}& {\beta }_{+}I-{\lambda }_{+}{S}_{-}-{\gamma }_{+}-{\mu }_{2}& -{\lambda }_{+}{S}_{+}\\ {\beta }_{-}{S}_{-}& -{\lambda }_{-}{S}_{-}& {\beta }_{-}I-{\lambda }_{-}{S}_{+}-{\gamma }_{-}-{\mu }_{3}\end{array}|$ (7)

Theorem 2. The boundary equilibrium points of (1) have the following local stability properties:

1) if ${R}_{0}<1$ , ${E}_{0}\left({I}_{0},0,0\right)$ is a stable node, and if ${R}_{0}>1$ , ${E}_{0}\left({I}_{0},0,0\right)$ is a saddle node.

2) ${E}_{1}\left({I}_{1},{S}_{+}^{1},0\right)$ is a stable if and only if ${R}_{1}^{*}>1$ .

3) ${E}_{2}\left({I}_{1},0,{S}_{-}^{2}\right)$ is a stable if and only if ${R}_{2}^{*}>1$ .

4) ${E}^{*}\left({I}^{*},{S}_{+}^{*},{S}_{-}^{*}\right)$ is always unstable if it exists.

Proof. 1) At ${E}_{0}$ , we have

${J}_{\left(2.2\right)}|{E}_{0}=|\begin{array}{ccc}-{\mu }_{1}& -\frac{\Lambda {\beta }_{+}}{{\mu }_{1}}& -\frac{\Lambda {\beta }_{-}}{{\mu }_{1}}\\ 0& \frac{\Lambda {\beta }_{+}}{{\mu }_{1}}-\left({\gamma }_{+}+{\mu }_{2}\right)& 0\\ 0& 0& \frac{\Lambda {\beta }_{-}}{{\mu }_{1}}-\left({\gamma }_{-}+{\mu }_{3}\right)\end{array}|$ (8)

The eigenvalues of the Equation (8) are ${\lambda }_{1}=-{\mu }_{1}$ , ${\lambda }_{2}=\left(\Lambda {\beta }_{+}-{\mu }_{1}\left({\gamma }_{+}+{\mu }_{2}\right)\right)/{\mu }_{1}=\left({\gamma }_{+}+{\mu }_{2}\right)/\left({R}_{01}-1\right)$ , ${\lambda }_{3}=\left(\Lambda {\beta }_{-}-{\mu }_{1}\left({\gamma }_{-}+{\mu }_{3}\right)\right)/{\mu }_{1}=\left({\gamma }_{-}+{\mu }_{3}\right)/\left({R}_{02}-1\right)$ , and therefore the ${E}_{0}$ is an unstable node if ${R}_{0}<1$ .

2) At ${E}_{1}$ , the Jacobian of (2) is

${J}_{\left(2.2\right)}|{E}_{1}=|\begin{array}{ccc}-\frac{\Lambda {\beta }_{+}}{{\gamma }_{+}+{\mu }_{2}}& -\left({\gamma }_{+}+{\mu }_{2}\right)& -\frac{{\beta }_{-}\left({\gamma }_{+}+{\mu }_{2}\right)}{{\beta }_{+}}\\ \frac{\Lambda {\beta }_{+}}{{\gamma }_{+}+{\mu }_{2}}-{\mu }_{1}& 0& -\frac{\left(\Lambda {\beta }_{1}-{\mu }_{1}\left({\lambda }_{+}+{\mu }_{2}\right)\right){\gamma }_{+}}{{\beta }_{+}\left({\lambda }_{+}+{\mu }_{2}\right)}\\ 0& 0& \frac{\left({R}_{1}^{*}-1\right)\left({\beta }_{-}\left({\lambda }_{+}+{\mu }_{2}\right)+{\gamma }_{-}{\mu }_{1}\right)}{{\beta }_{+}}\end{array}|$ (9)

The eigenvalues are ${\lambda }_{1}=\left(1-{R}_{1}^{*}\right)\left({\beta }_{+}\left({\gamma }_{-}+{\mu }_{3}\right)\right)/{\beta }_{2}$ , ${\lambda }_{2,3}=-\frac{{\mu }_{1}\left({R}_{01}±\sqrt{{R}_{01}^{2}-\left({R}_{01}-1\right)\left({\lambda }_{+}+{\mu }_{2}\right)/{\mu }_{1}}}{2}$ , obviously they all have negative real parts and ${E}_{1}$ is local stability if ${R}_{1}^{*}>1$ and ${R}_{01}>1$ .

In a similar fashion, we can conclude that if ${R}_{2}^{*}>1$ and ${R}_{02}>1$ . ${E}_{2}$ is local stability.

For the positive equilibrium ${E}^{*}$ ,

$|\begin{array}{ccc}{\frac{\Lambda }{I}}^{*}& -{\beta }_{+}{I}^{*}& -{\beta }_{-}{I}^{*}\\ \frac{{\beta }_{-}{I}^{*}-\left({\lambda }_{-}+{\mu }_{3}\right){\beta }_{+}}{{\beta }_{-}}& 0& -\frac{\left({\beta }_{-}{I}^{*}-\left({\lambda }_{-}+{\mu }_{3}\right)\right){\gamma }_{+}}{{\beta }_{-}}\\ \frac{\left({\beta }_{+}{I}^{*}-\left({\lambda }_{+}+{\mu }_{2}\right)\right){\beta }_{-}}{{\beta }_{+}}& -\frac{\left({\beta }_{+}{I}^{*}-\left({\lambda }_{+}+{\mu }_{2}\right)\right){\gamma }_{-}}{{\beta }_{+}}& 0\end{array}|$ (10)

The characteristic polynomial of (10) can be calculated as follows is

${\lambda }^{3}+\frac{\Lambda }{{I}^{*}}{\lambda }^{2}+\left({\beta }_{+}^{2}{S}_{+}^{*}{I}^{*}+{\beta }_{-}^{2}{S}_{-}^{*}{I}^{*}-{\lambda }_{+}{\lambda }_{-}{S}_{+}^{*}{S}_{-}^{*}\right)\lambda -f=0$ (11)

which $f=\left({\beta }_{+}{\beta }_{-}\left({\lambda }_{+}+{\lambda }_{-}\right){S}_{+}^{*}{S}_{-}^{*}\right)+\frac{\Lambda {S}_{+}^{*}{S}_{-}^{*}{\lambda }_{+}{\lambda }_{-}}{{I}^{*}}$ , the eigenvalues satisfied

${\lambda }_{1}{\lambda }_{2}{\lambda }_{3}=f>0$ , It can be obtained by the Routh-Hurwitz criterion: It exists at least one eigenvalue of the polynomial equation has positive real part, then ${E}^{*}$ is always unstable.,

Next, we will consider the global stability of quiet equilibrium ${E}_{0}$ of system (2).

Theorem 3. The quiet equilibrium ${E}_{0}=\left({I}_{0},0,0\right)$ is globally asymptotically stable provided that ${R}_{0}<1$ .

Proof. Define the Lyapunov function

$V\left(I,{S}_{+},{S}_{-}\right)={V}_{1}\left(I,{S}_{+},{S}_{-}\right)+{V}_{2}\left(I,{S}_{+},{S}_{-}\right)$ (12)

where ${V}_{1}\left(I,{S}_{+},{S}_{-}\right)=I-{I}_{0}-{I}_{0}ln\frac{I}{{I}_{0}}$ , ${V}_{2}\left(I,{S}_{+},{S}_{-}\right)={S}_{+}+{S}_{-}$ , then differentiating ${V}_{1}$ with respect to t along solutions of system (2) gives

$\begin{array}{c}\frac{\text{d}{V}_{1}}{\text{d}t}=\left(1-\frac{{I}_{0}}{I}\right)\left(\Lambda -{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I-{\mu }_{1}I\right)\\ =\left(1-\frac{{I}_{0}}{I}\right)\left(-{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I\right)+{\mu }_{1}{I}_{0}\left(1-\frac{{I}_{0}}{I}\right)\left(1-\frac{I}{{I}_{0}}\right)\\ \le \left(1-\frac{{I}_{0}}{I}\right)\left(-{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I\right)\\ =-{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I+{\beta }_{-}{S}_{-}{I}_{0}+{\beta }_{+}{S}_{+}{I}_{0}\end{array}$ (13)

$\frac{\text{d}{V}_{2}}{\text{d}t}={\beta }_{+}{S}_{+}I-{\lambda }_{+}{S}_{-}{S}_{+}-{\gamma }_{+}{S}_{+}-{\mu }_{2}{S}_{+}+{\beta }_{-}{S}_{-}I-{\lambda }_{-}{S}_{-}{S}_{+}-{\gamma }_{-}{S}_{-}-{\mu }_{3}{S}_{-}$ (14)

Then, we have

$\begin{array}{c}\frac{\text{d}V}{\text{d}t}=\frac{\text{d}{V}_{1}}{\text{d}t}+\frac{\text{d}{V}_{2}}{\text{d}t}\\ \le {\beta }_{+}{S}_{+}{I}_{0}-{\lambda }_{+}{S}_{-}{S}_{+}-\left({\gamma }_{+}+{\mu }_{2}\right){S}_{+}+{\beta }_{-}{S}_{-}{I}_{0}-{\lambda }_{-}{S}_{-}{S}_{+}-\left({\gamma }_{-}-{\mu }_{3}\right){S}_{-}\\ =\left({\gamma }_{+}+{\mu }_{2}\right)\left({R}_{01}-1\right){S}_{+}+\left({\gamma }_{-}+{\mu }_{3}\right)\left({R}_{02}-1\right){S}_{-}-\left({\lambda }_{+}+{\lambda }_{-}\right){S}_{-}{S}_{+}\\ \le \left({\gamma }_{+}+{\mu }_{2}\right)\left({R}_{0}-1\right){S}_{+}+\left({\gamma }_{-}+{\mu }_{3}\right)\left({R}_{0}-1\right){S}_{-}-\left({\lambda }_{+}+{\lambda }_{-}\right){S}_{-}{S}_{+}\end{array}$ (15)

Therefore, when ${R}_{0}<1$ , $\frac{\text{d}V}{\text{d}t}\le 0$ , and the equality holds only for $I={I}_{0}$ ,

${S}_{+}=0$ and ${S}_{-}=0$ , According to the LaSalle invariant principle, we have
$li{m}_{t\to \infty }I\left(t\right)={I}_{0}$ , $li{m}_{t\to \infty }{S}_{+}\left(t\right)=0$ and $li{m}_{t\to \infty }{S}_{-}\left(t\right)=0$ . Thus, the proof is completed.,

4. Numerical Studies

4.1. Numerical Simulation

We have formulated a S2IR rumor-spread model. We have established some threshold conditions for the steady state. To illustrate the analytical results, we do some numerical simulations. From Table 1, we can set ${\mu }_{1}=0.3$ , ${\mu }_{2}=0.3$ , ${\mu }_{3}=0.3$ , ${\gamma }_{+}=0.4$ , ${\gamma }_{-}=0.4$ , ${\lambda }_{+}=0.38$ , ${\lambda }_{-}=0.1$ , then system (2) becomes

$\begin{array}{l}\stackrel{˙}{I}=\Lambda -{\beta }_{-}{S}_{-}I-{\beta }_{+}{S}_{+}I-0.3I\\ {\stackrel{˙}{S}}_{+}={\beta }_{+}{S}_{+}I-0.38{S}_{-}{S}_{+}-0.4{S}_{+}-0.3{S}_{+}\\ {\stackrel{˙}{S}}_{-}={\beta }_{-}{S}_{-}I-0.1{S}_{-}{S}_{+}-0.4{S}_{-}-0.3{S}_{-}\end{array}$ (16)

Let $\Lambda =0.2$ , ${\beta }_{+}=0.705$ , ${\beta }_{-}=0.725$ , then ${R}_{0}=0.6905<1$ . By theorem 1.(a), theorem 2.(a) and theorem 3, system (16) has only the trivial equilibrium ${E}_{0}\left(0.6667,0,0\right)$ and it is globally asymptotically stable (see Figure 2), it means that the message(rumor or truth) will eliminate.

Let $\Lambda =0.312$ , ${\beta }_{+}=0.856$ , ${\beta }_{-}=0.656$ . We have ${R}_{01}=1.272>1$ , ${R}_{02}=0.974<1$ , ${R}_{1}^{*}=1.303>1$ , ${R}_{2}^{*}=0.7996<1$ , It presents the system (16) has have two equilibria, ${E}_{0}\left(1.04,0,0\right)$ and ${E}_{1}\left(0.8178,0.095,0\right)$ where ${E}_{0}$ is unstable but ${E}_{1}$ is stable according to theorem 1.(b) and theorem 2.(a)-2.(b) (see Figure 3(a)). Let $\Lambda =0.323$ , ${\beta }_{+}=0.6215$ , ${\beta }_{-}=0.8545$ , then system ${R}_{01}=0.956<1$ , ${R}_{02}=1.314>1$ , ${R}_{1}^{*}=0.738<1$ , ${R}_{2}^{*}=1.362>1$ . It shows that the system has

Figure 2. The trajectories, equilibria and stability of system (4.1) with $\Lambda =0.2,\text{\hspace{0.17em}}{\beta }_{+}=0.705,\text{\hspace{0.17em}}{\beta }_{-}=0.725$ .

(a)(b)

Figure 3. the trajectories, equilibria and stability of system (16) with (a) $\Lambda =0.323,\text{\hspace{0.17em}}{\beta }_{+}=0.6215,\text{\hspace{0.17em}}{\beta }_{-}=0.8545$ , (b) $\Lambda =0.323,\text{\hspace{0.17em}}{\beta }_{+}=0.6215,\text{\hspace{0.17em}}{\beta }_{-}=0.8545$ .

two equilibria, ${E}_{0}\left(1.077,0,0\right)$ and ${E}_{2}\left(0.819,0,0.11\right)$ where ${E}_{0}$ is unstable but ${E}_{2}$ is stable according to theorem 1.c and theorem 2.a, 2.c (see Figure 3(b)).

Let $\Lambda =2.1$ , ${\beta }_{+}=0.6$ , ${\beta }_{-}=0.9$ , then ${R}_{0}=9>0$ , ${R}_{1}^{*}=0.9091<1$ , ${R}_{2}^{*}=3.101>1$ . By theorem 1.(d), system (16) has three equilibria ${E}_{0}\left(7,0,0\right)$ , ${E}_{1}\left(1.167,2.5,0\right)$ and ${E}_{2}\left(0.778,0,2.667\right)$ and ${E}_{0}$ , ${E}_{1}$ is unstable, but ${E}_{2}$ is stable (see Figure 4(a)). Let $\Lambda =2.8$ , ${\beta }_{+}=0.76$ , ${\beta }_{-}=0.72$ , then

${R}_{0}=10.13>0$ , ${R}^{*}={R}_{1}^{*}=1.566>1$ , ${R}_{2}^{*}=2.474>1$ . By theorem 1.(e), system (16) has three equilibria ${E}_{0}\left(9.33,0,0\right)$ , ${E}_{1}\left(0.921,3.605,0\right)$ , ${E}_{2}\left(0.972,0,3.583\right)$ and ${E}^{*}\left(1.244,1.956,0.646\right)$ , which ${E}_{0}$ and ${E}^{*}$ are not stable but ${E}_{1}$ and

(a)(b)

Figure 4. The trajectories, equilibria and stability of system (16) with (a) $\Lambda =2.1,\text{\hspace{0.17em}}{\beta }_{+}=0.6,\text{\hspace{0.17em}}{\beta }_{-}=0.9$ (b) $\Lambda =2.8,\text{\hspace{0.17em}}{\beta }_{+}=0.76,\text{\hspace{0.17em}}{\beta }_{-}=0.72$ .

${E}_{2}$ are bistable. The phase portrait is represented in Figure 4(b).

4.2. Uncertainty and Sensitivity Analysis

In this part, uncertainty and sensitivity analysis, based on Latin hypercube sampling (LHS) and partial rank correlation coefficients (PRCC) scheme, will explore the dependence of ${R}_{0}$ , ${R}_{1}^{*}$ and ${R}_{2}^{*}$ due to the change of input parameters in the estimation of the uncertainty. Each input parameters is sampling 1000 times. A uniform distribution function was used and the variation ranges of all parameters are given in Table 1.

Table 1. Input parameter sample values for simulation.

Table 2. Input parameter sample values for simulation.

The PRCC results which illustrate the dependence on each parameters of ${R}_{0}$ , ${R}_{1}$ , ${R}_{2}$ respectively. We considered $|\text{PRCC}|\ge 0.4$ as indicating that the high correlation between input parameters and output variables, $0.2\le |\text{PRCC}|<0.4$ as moderate correlations and $|\text{PRCC}|<0.2$ as no relations.

The results of the simulation of PRCCs are shown in Table 2 and Figure 5. For ${R}_{0}$ , the absolute values of PRCCs of each parameters fall below 0.2. It suggests that the impact of ${R}_{0}$ is the result of the interaction of each parameters and it is very difficult to decrease ${R}_{0}$ by regulating and controlling a few parameters. The parameters with the most positive impact on ${R}_{1}$ are the fade away rate of truth-spreaders ${\lambda }_{-}$ and the force of infection of rumor-spreaders ${\beta }_{+}$ , while the most negative impact on ${R}_{1}$ are the natural fade away of truth-spreaders ${\mu }_{2}$ , the forgetting rate of rumor-spreader ${\gamma }_{+}$ and the force of infection of truth-spreader ${\beta }_{-}$ . It means that the five factors play a critical role on the process of rumor-spreading. The parameters that have the moderate influence with ${R}_{1}$ are the recruitment $\Lambda$ , the natural fade away of truth-spreader ${\mu }_{3}$ and the forgetting rate of truth-spreader ${\gamma }_{-}$ . Compare with ${R}_{2}$ , we find the parameters ${\beta }_{+}$ , ${\beta }_{-}$ , ${\mu }_{2}$ , ${\mu }_{3}$ , ${\gamma }_{+}$ , ${\gamma }_{-}$ have opposite influences on ${R}_{1}$ and ${R}_{2}$ . This exhibitions the propagation have the nature of competitive exclusion between rumor with truth. The parameters of the fade away rate ${\lambda }_{+}$ , ${\lambda }_{-}$ have no relations with ${R}_{1}$ , ${R}_{2}$ respectively, but have high positive relations with ${R}_{2}$ , ${R}_{1}$ respectively. If both ${\lambda }_{+}$ and ${\lambda }_{-}$ are large, it is possible that both ${R}_{1}$ and ${R}_{2}$ are greater than 1 and this will lead to the occurrence of bistable on ${E}_{1}$ and ${E}_{2}$ .

Figures 6(a)-(c) shows the LHSs frequency distributions. The mean value for ${R}_{0}$ , ${R}_{1}$ and ${R}_{2}$ are ${M}_{{R}_{0}}=2.4214$ , ${M}_{{R}_{1}}=1.4797$ , ${M}_{{R}_{2}}=1.4793$ respectively. The standard deviation with ${R}_{0}$ is 0.6642 and greater than with ${R}_{1}$ and ${R}_{2}$ (0.3120, 0.3037 respectively) which indicates that the derived frequency distribution for ${R}_{0}$ is dispersed than the distributions for ${R}_{1}$ and ${R}_{2}$ .

5. Conclusion

Rumor propagation has been investigated through different types of mathematical models. In our study, we consider a rumor propagation model with truth-spreading and determine the threshold which governs the dynamics of the system. On one hand, the propagation of rumor and truth is mutually

Figure 5. PRCCs illustrating the dependence of ${R}_{0},\text{\hspace{0.17em}}{R}_{1}^{*}$ and ${R}_{2}^{*}$ for the rumor spread model on each parameter. Symbol * shows the PRCC value is not zero significantly ( $P<0.05$ ).

Figure 6. Uncertainty results based on Latin hypercube sampling. The top row consists of the frequency diagrams for (a) ${R}_{0}$ , (b) ${R}_{1}^{*}$ and (c) ${R}_{2}^{*}$ . The following row shows estimates of CDFs for output variables ((d) ${R}_{0}$ , (e) ${R}_{1}^{*}$ and (f) ${R}_{2}^{*}$ ).

exclusive. Therefore, the appropriate increase in the spread of truth is conducive to the elimination of the spread of rumors. On the other hand, the spread of rumor and truth can coexist for a long time under certain conditions. At this point, the improvement of the spread of the truth is also conducive to the spread of rumors in a certain range. This requires the government and related organizations to further develop the corresponding work.

Acknowledgements

The research has been supported by The Social Science Planning Projects of Jiangxi Province (14XW08) and The 12th Five-year Education Scientific Planning Project of Jiangxi Province (15ZD3LYB031).

Conflicts of Interest

The authors declare no conflicts of interest.