Rumor spread in social network

Rumor is an important form of social communications, and spread of rumors plays a significant role in a variety of human affairs. There are two rumor models that are widely used, i.e. DK model and MK model. Particularly, we can view rumor spread as a stochastic process in social networks.

Model

Spread of Rumors

An standard model of rumor spreading was introduced by Daley and Kendall, which is called DK model. Assume there are N people in total. And those people in the network are categorized into three groups: ignorants, spreaders and stiflers, which are denoted as S, I, and R respectively hereinafter:

• S: people who are ignorant of the rumor;
• I: people who actively spread the rumor;
• R: people who have heard the rumor, but no longer are interested in spreading it.

The rumor is propagated through the population by pair-wise contacts between spreaders and others in the population. Any spreader involved in a pair-wise meeting attempts to “infect” the other individual with the rumor. In the case this other individual is an ignorant, he or she becomes a spreader. In the other two cases, either one or both of those involved in the meeting learn that the rumor is known and decided not to tell the rumor anymore, thereby turning into stiflers.

One famous variant is Maki-Thompson(MK) model. In this model, rumor is spread by directed contacts of the spreaders with others in the population. Furthermore, when a spreader contacts another spreader only the initiating spreader becomes a stifler. Therefore, three types of interactions can happen with certain rates.

1. Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): S+I\xrightarrow{\alpha}2I

1. which says when a spreader meet an ignorant, the ignorant will become a spreader.

2. Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): I+I\xrightarrow{\beta}I+R

1. which says when two spreaders meet with each other, one of them will become a stifler.

2. Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): I+R\xrightarrow{\beta}2R

1. which says when a spreader meet a stifler, the spreader will lose the interest in spreading the rumor, so become a stifler.

Of course we always have conservation of individuals:
$N=I+S+R$

The change in each class in a small time interval is:
$\Delta S \approx - \Delta t \alpha IS/N$
$\Delta I \approx \Delta t [{\alpha IS \over N} - {\beta I^2 \over N} - {\beta IR \over N}]$
$\Delta R \approx \Delta t [{\beta I^2 \over N}+{\beta IR \over N}]$

Since we know $S$, $I$ and $R$ sum up to $N$, we can reduce one equation from the above, which leads to a set of differential equations using relative variable $x=I/N$ and $y=S/N$ as follows
${dx \over dt} = x \alpha y - \beta x^2 - \beta x(1-x-y)$
${dy \over dt} = - \alpha xy$
which we can write ${dx \over dt} = (\alpha + \beta)xy - \beta x$
${dy \over dt} = - \alpha xy$

Compared with the ordinary SIR model, we see that the only difference to the ordinary SIR model is that we have a factor $\alpha + \beta$ in the first equation instead of just $\alpha$. We immediately see that the ignorants can only decrease since $x,y\ge 0$ and ${dy \over dt}\le 0$. Also, if
$R_0={\alpha +\beta \over \beta} >1$
which means
${\alpha \over \beta}>0$
the rumour model exhibits an “epidemic” even for arbitrarily small rate parameters.

Rumor Spread in Social Network

We model the process introduced above on a network in discrete time, that is, we can model it as a DTMC. Say we have a network with N nodes, then we can define $X_i(t)$ to be the state of node i at time t. Then $X(t)$ is a stochastic process on $S=\{S,I,R\}^N$. At a single moment, some node i and node j interact with each other, and then one of them will change its state. Thus we define the function $f$ so that for $x$ in $S$,$f(x,i,j)$ is when the state of network is $x$, node i and node j interact with each other, and one of them will change its state. The transition matrix depends on the number of ties of node i and node j, as well as the state of node i and node j. For any $y=f(x,i,j)$, we try to find $P(x,y)$. If node i is in state I and node j is in state S, then $P(x,y)=\alpha A_{ji}/k_i$; if node i is in state I and node j is in state I, then $P(x,y)=\beta A_{ji}/k_i$; if node i is in state I and node j is in state R, then $P(x,y)=\beta A_{ji}/k_i$. For all other $y$, $P(x,y)=0$.
The procedure on a network is as follows:

1. We initial rumor to a single node $i$;
2. We pick one of its neighbors as given by the adjacency matrix, so the probability we will pick node $j$ is
   

   $p_j={A_{ji} \over k_i}$
   where $A_{ji}$ is from the adjacency matrix and $A_{ji}=1$ if there is a tie from $i$ to $j$, and Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): k_i= \textstyle \sum_{j=1}^N A_{ij}

   is the degree for node $i$;
3. Then have the choice:
   1. If node $j$ is an ignorant, it becomes a spreader at a rate $\alpha$;
   2. If node $j$ is a spreader or stifler, then node $i$ becomes a stifler at a rate $\beta$.
4. We pick another node who is a spreader at random, and repeat the process.

We would expect that this process spreads the rumor throughout a considerable fraction of the network. Note however that if we have a strong local clustering around a node, what can happen is that many nodes become spreaders and have neighbors who are spreaders. Then, every time we pick one of those, they will recover and can extinguish the rumor spread. On the other hand, if we have a network that is small world, that is, a network in which the shortest path between two randomly chosen nodes is much smaller than that one would expect, we can expect the rumor spread far away.

Also we can compute the final number of people who once spread the news, this is given by
$r_\infty=1-e^{-({\alpha +\beta \over \beta})r_\infty}$
In networks the process that does not have a threshold in a well mixed population, exhibits a clear cut phase-transition in small worlds. The following graph illustrates the asymptotic value of $r_\infty$ as a function of the rewiring probability $p$.