Biological Modelling Project
b11902021 張嘉崴
Assignment 5
\lambda_{mx} = 0.951
Intro
Why didn't the Aztec Empire eradicate Europe?
Why didn't the Aztec Empire's
viruses eradicate Europe?
- population?
- difference in military forces?
- difference in tamable livestock.
Why didn't the Aztec Empire's
viruses eradicate Europe?
More livestocks
Higher probability of new viruses that can pass between humans emerge
The community with
stronger viruses wins
Really?
The community with
stronger viruses wins
Really?
001010
011000
The hamming distance should be 2
My vague blueprint
Make 2 communities A and B
Add tons of viruses/mutations with different susceptibilities, fatalities, transmission rate ... into A
Drop a fraction of people in A into B
How to mutate?
How to simulate the effect of mutation on the viruses' abilities to transmit, to kill...?
Model Design
- Individual Based (virus's perspective)
- SIR Model
- Evolutionary
Model Design
- Individual Based (virus's perspective)
- SIR Model
- Evolutionary
Model Design
for i in range(year):
human_reproduce()
human_move()
virus_simulation()
virus_infection()
human reproduction
n =n_s + n_i + n_r
\\
n_s \leftarrow n_s + kn_r + rn(1-\frac{n}{K})
what if \( n > K \)
(human move)
v_0
v_1
v_2
v_3
n_i \leftarrow n_i - deg(v_i) p + p\sum_{v_j \in N(v_i)} n_j
virus kill/stay/recover
Every virus has one real value \(g_a\)
Every node has one real value \(g_b\)
E = \frac{1}{1+e^{c(g_b-g_a)}}
E - \epsilon \,/\, 2\epsilon \,/\,(1-E) - \epsilon
deck
By willy62830
deck
- 157