Two Species Model of Epidemic Spreading
description
Transcript of Two Species Model of Epidemic Spreading
![Page 1: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/1.jpg)
Two Species Model of Epidemic Spreading
Ahn, Yong-Yeol. 2005.4.14
![Page 2: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/2.jpg)
CSSPL, http://stat.kaist.ac.kr
Epidemics and Contact Process
• Epidemics is about self-replicators which can die.
• Contact process : A simple model consists of branching and spontaneous annihilation.
![Page 3: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/3.jpg)
CSSPL, http://stat.kaist.ac.kr
Contact Process
1
λ
v
λ =ν/ δ
![Page 4: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/4.jpg)
CSSPL, http://stat.kaist.ac.kr
Transition of C. P. (1d)
![Page 5: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/5.jpg)
CSSPL, http://stat.kaist.ac.kr
Contact Process
iiidt
di)1(
Approx. Malthus–Verhulst eq.
For λ < 0 infection vanishes.For λ > 0 steady state ( i = 1 – 1/ λ)
This equation is “wrong”, but correctly describes the qualitative dynamics of contact process, which is that…
Contact process exhibits an absorbing phase transition.
We can NOT solve exactly this terribly simple process!
![Page 6: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/6.jpg)
CSSPL, http://stat.kaist.ac.kr
Absorbing Phase Transition
• Order parameter
• Correlation length
• Correlation time
ρ
p
activeinactive
cpp :
x
cx pp :
t
ct pp :
![Page 7: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/7.jpg)
CSSPL, http://stat.kaist.ac.kr
Directed Percolation
• Contact process is same as the directed percolation process.
• In the lattice with randomly occupied edges,
• Follow only occupied edges in only designated direction.
• This process can continue ad infinitum?
Directed Percolation problem
![Page 8: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/8.jpg)
CSSPL, http://stat.kaist.ac.kr
Surface Roughening Model and DP
XNo spontaneous infection no digging condition
- Let’s see only one layer
![Page 9: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/9.jpg)
CSSPL, http://stat.kaist.ac.kr
S. R. and Coupled DP
A A A A A
B
D D
C
B B B B B B B B
C C C C C C C CCC
D D D D D D D D DDD
C
(RSOS : Restricted solid-on-solid interface)
active
inactive
B
![Page 10: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/10.jpg)
CSSPL, http://stat.kaist.ac.kr
Asymmetrically Coupled DP
A B C D ..
Positive coupling
Simulation results in 1D
A : DP critical behavior B, C, … : dressed negative coupling is irrelevant
Negative coupling
![Page 11: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/11.jpg)
CSSPL, http://stat.kaist.ac.kr
ACDP : Diagram
BA
AntigenContagion
AntibodyWorm killing worm
Detection and initiation
Healing
![Page 12: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/12.jpg)
CSSPL, http://stat.kaist.ac.kr
ACDP on Networks
• Antigen – antibody interaction• Worm – worm killing worm
• Motto:
“Never flood the network”
![Page 13: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/13.jpg)
CSSPL, http://stat.kaist.ac.kr
Phase Diagram
• If there is no virus, vaccine spreading has a epidemic threshold.
• If there is no vaccine, virus has no epidemic threshold.
Ap
0
0A
B
Bp
0
0A
B
0
0A
B
Ma
c
b
There is the happy region where there is no virus and no vaccine.
Vaccine Annihilation probability
Virus Annihilation probability
![Page 14: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/14.jpg)
CSSPL, http://stat.kaist.ac.kr
Mean Field Theory
ABBBB
BAAAAA
a
a
~
,~
2
2
Aa
Ba
0A B
0
/A
B B Ba b
0
0A
B
M
DP
DP
DP
2, 1/ 2, 1
1, DP1/ 2, 1A xA tA
B xB tB
At ‘M’
Mean field eq.
![Page 15: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/15.jpg)
CSSPL, http://stat.kaist.ac.kr
On a Scale-Free Network
• We should consider each degree seperately.
![Page 16: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/16.jpg)
CSSPL, http://stat.kaist.ac.kr
SIS Model on a Scale-free Network
))(())(1()()(
titiktit
tikk
k
Because a scale-free network is very inhomogenous, we must consider the density of infected nodes seperately for each degree.
))(( ti is the probability that any given link points to an infected node.
(We don’t need to consider the excess degree because we’re dealing with SIS model. )
![Page 17: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/17.jpg)
CSSPL, http://stat.kaist.ac.kr
SIS Model on a Scale-free Network
))(())(1()()(
titiktit
tikk
k
))(( ti
The probability that the selected node is not infected
![Page 18: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/18.jpg)
CSSPL, http://stat.kaist.ac.kr
SIS Model on a Scale-free Network
At Steady state,
Where Average degree
Edge selection
With these two equations and continum approximation,
We get
ρk = ik
![Page 19: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/19.jpg)
CSSPL, http://stat.kaist.ac.kr
SIS Model on a Scale-free Network
• By solving previous equation, we get, “no epidemic threshold
on the scale-free network”
Homogeneous network
Scale-freeVanishing
Epidemic threshold
![Page 20: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/20.jpg)
CSSPL, http://stat.kaist.ac.kr
CDP on a Scale-free Network
spontaneous annihilationwith prob. Ap
A
(1 )
infecting or activating vaccine at “all neighbors” with prob. 1 Ap
B
branching “one offspring”with prob. 1 Bp
spontaneous annihilationwith prob. Bp
From J.D. Noh’s ppt
predation
![Page 21: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/21.jpg)
CSSPL, http://stat.kaist.ac.kr
Different Branching Rule
• Because we do not want the vaccine occupying the whole network and vaccine is more active in the condition that virus is abundant, we assign a different rule for the vaccine.
• Both rules are possible in our CDP framework
Previous SIS rule Modified rule
![Page 22: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/22.jpg)
CSSPL, http://stat.kaist.ac.kr
On a Scale-Free Network
kk
Bk
kk
Ak
BkkB
AkkAkBk
AkA
BkkBkAk
k
kkPykkPx
kypkypypy
kpxpxpx
)|( ,)|(
where,
)1)(1()1()1(
)1)(1()1(
Mean field equation
![Page 23: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/23.jpg)
CSSPL, http://stat.kaist.ac.kr
When There is no Virus
kk
Bk
kk
Ak
BkkB
AkkAkBk
AkA
BkkBkAk
k
kkPykkPx
kypkypypy
kpxpxpx
)|( ,)|(
where,
)1)(1()1()1(
)1)(1()1(
' '''
)'(
'
)|'(
k kkk
Bk k
kPy
k
kkPy
Without degree-degree correlation,
Mean field eq.
![Page 24: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/24.jpg)
CSSPL, http://stat.kaist.ac.kr
Mean Field eq.
'
' )ky/)'(
(where )1)(1(k
kBkkBkBk k
kPyy
k
kypypy
By assuming a steady state, multiplying p(k), and summing over k, we get a equation
kp
pky
B
Bk 1
21
Note that the critical point pB=½.
![Page 25: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/25.jpg)
CSSPL, http://stat.kaist.ac.kr
MF eq.
• Calculating with 2<γ<3,
) ,1
(where
,1
)( 1
2
kuk
y
p
pu
uuPduky
B
B
k
kk
m
And by direct calculation,
1
)21(~ Bpy pc
As γ decreases, β increases.
![Page 26: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/26.jpg)
CSSPL, http://stat.kaist.ac.kr
MF eq.
• Calculating with γ=3 β = 1 with logarithmic correction
• Calculating with γ>3, β = 1
γ2 3
β = 1/(γ-2) β = 1
β = 1 w/ log. correction
![Page 27: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/27.jpg)
CSSPL, http://stat.kaist.ac.kr
Diffusion – Annihilation Process
'' )
)'( (where )1)(1(
kk
Bk
BkkBkBk k
kPykypypy
Out CDP eq.
by multiflying P(k) and summing over k,
2/1Bp
k
kyypypypy kBBB )1()1(
When ,
k
kyyy k
2
1
![Page 28: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/28.jpg)
CSSPL, http://stat.kaist.ac.kr
Diffusion – Annihilation Process
),)](21[)()(
ttk
kt
dt
tdkk
k
))(2)(
ttdt
td
And by summing over P(k),
Meanwhile, diffusion-annihilation process can be described as,
k
kyyy k
2
1
Same eq.
![Page 29: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/29.jpg)
CSSPL, http://stat.kaist.ac.kr
DA Process and Time Evolution
• At critical point, DA process and our model is described by the same mf eq. And we apply the DA result to our model.
γ2 3
tt
1~)(
ttt
ln
1~)(
1
~)( tt
![Page 30: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/30.jpg)
CSSPL, http://stat.kaist.ac.kr
Simulation Result
Critical point = 0.425 (mean field : 0.5)
δ = 1.3 (mf: 1 with log. correction )
γ =3
δ=1.3
If we draw with logarithmic term, δ ~ 1.1
![Page 31: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/31.jpg)
CSSPL, http://stat.kaist.ac.kr
Finite Size Scaling
p
p
pppp
c
c
~
)( ~
)(
)( ~),(
LL
LpLt
Assuming a scaling form,
LL
LpLpfLLp
:
: )(),(
( θ=β/ ν )
![Page 32: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/32.jpg)
CSSPL, http://stat.kaist.ac.kr
Finite Size Scaling
β =1.3
γ =3
θ~β/ ν
Mf: β = 1 with log. correction
![Page 33: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/33.jpg)
CSSPL, http://stat.kaist.ac.kr
Another γ s
• γ =2.5 Critical point = 0.425 (mean field : 0.5)
δ = 2.0 (mf: 1/(γ-2) = 2 )
![Page 34: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/34.jpg)
CSSPL, http://stat.kaist.ac.kr
Another γ s
• γ =2.5
β = 1.75 (mf: 2)
![Page 35: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/35.jpg)
CSSPL, http://stat.kaist.ac.kr
Another γ s
• γ =2.75, 4
δ ~ 1.5
pc = 0.425 pc = 0.426
δ ~ 1.3
![Page 36: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/36.jpg)
CSSPL, http://stat.kaist.ac.kr
Summary : δ
γ2 3 1
cor.) log.(w/ 1
1
~1.3 (mf: 2)
~1.3 (1.1 w l.c.) (mf:1 w lc)
2 (mf: 2)
2.5 2.75 4
~1.5(mf:1.33)
![Page 37: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/37.jpg)
CSSPL, http://stat.kaist.ac.kr
Summary : β
γ2 3β = 1/(γ-2) β = 1
β = 1 w/ log. correction
~1.3
1.75(mf:2)
2.5 2.75
Not yet Not yet
4
![Page 38: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/38.jpg)
CSSPL, http://stat.kaist.ac.kr
Summary
• We put the CDP model on a scale-free network for simulating the behavior of an epidemic and reacting anti-pathogen.
• The simulation results are consistent to the mean field calculation.
![Page 39: Two Species Model of Epidemic Spreading](https://reader035.fdocument.org/reader035/viewer/2022062500/56815013550346895dbdf899/html5/thumbnails/39.jpg)
CSSPL, http://stat.kaist.ac.kr
References
• U. C. Täuber et al., Phys. Rev. Lett. 80, 2165 (1998); Y. Y. Goldschmidt, ibid. 81, 2178 (1998).
• J. D. Noh and H. Park, CondMat/051542 (to be published in Phys. Rev. Lett.).
• J. Marro and R. Dickman, “Nonequilibrium Phase Transition in Lattice Models”.
• Papers in http://janice.kaist.ac.kr/bbs/viewtopic.php?t=30 and related papers.
• Haye Hinrichsen, “Nonequilibrium Critical Phenomena and Phase Transitions into Absorbing States”, arXiv:cond-mat/0001070.
• Arno Karlen, “Man and Microbes”.• J. D. Murray, “Mathematical Biology”.• M. E. J. Newman, SIAM Review, 45, 167 (2003).• Malcolm Gladwell, “The Tipping Point”.• About percolation : 손승우 , http://stat.kaist.ac.kr/~sonswoo/prl.html .