Difference between revisions of "Team:SZU-China/Model"

 
(36 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
{{SZU-China/Nav}}
+
{{SZU-China/header}}
 
<html>
 
<html>
  
−
<style>
+
 
−
h1 {
+
<style>
−
color: #3880A0;
+
#icon{
−
}
+
transform:translateY(-49px);
−
+
}
−
h2 h3 h4 {
+
−
font-weight: bold;
+
.border-3 {
−
}
+
    border-width:3px !important;
−
+
   
−
.row {
+
}
−
padding: 48px;
+
.card[class*="border"] {
−
}
+
    border: 1px solid #469789;
−
+
    -webkit-box-shadow: none;
−
body {
+
    box-shadow: none;
−
position: relative;
+
−
font-family: 'Segoe UI', sans-serif !important;
+
−
font-weight: 400;
+
−
}
+
−
+
−
@font-face {
+
−
font-family: 'Segoe UI';
+
−
src: url('https://static.igem.org/mediawiki/2018/b/b8/T--SZU-China--segoeui.ttf');
+
−
}
+
−
+
−
#HQ_page p {
+
−
font-size: 20px;
+
−
color: #393939;
+
−
}
+
−
hr {
+
−
width: auto;
+
−
border-top: 1px solid rgba(0, 0, 0, 0.125);
+
−
}
+
−
+
−
.container-full{
+
−
padding: 0px;
+
−
}
+
−
body{
+
−
padding-top: 40px;
+
 
}
 
}
−
</style>
 
  
 +
.card-text{
 +
font-size: 20px !important;
 +
}
 +
#HQ_page p {
 +
   
 +
    text-align: center;
 +
}
 +
#card1 > div.card-body.text-center{
 +
padding-top: 0px;
 +
}
 +
#card2 > div.card-body.text-center{
 +
padding-top: 0px;
 +
}
 +
#card1 > div.view.text-center{
 +
height: 60px;
 +
}
 +
#card2 > div.view.text-center{
 +
height: 60px;
 +
}
 +
</style>
 +
 +
 
<body data-spy="scroll" data-target="#myScrollspy" data-offset="10">
 
<body data-spy="scroll" data-target="#myScrollspy" data-offset="10">
−
<div class="container-full">
+
<div class="container-full">
−
+
 
<div class="card">
 
<div class="card">
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/8/81/T--SZU-China--pre_design.jpg" />
+
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/f/fc/T--SZU-China--Model-home.jpg" />
 
</div>
 
</div>
 
</div>
 
</div>
 +
 
<div class="container-fluid">
 
<div class="container-fluid">
−
<div class="row" style="background-image: url();">
+
<div class="col-10 offset-1 shadow">
−
<div class="col-9 offset-1 shadow ">
+
<div class="row"></div>
−
<div class="text-center">
+
<div class="text-center">
−
<h1 id="header">Epidemic Model</h1>
+
<h1 class="h1" style="color: #469789; text-align: center;">Model</h1>
−
</div>
+
</div>
−
 
+
<div class="row">
−
<div class="row">
+
−
 
+
−
<h2>Introduction</h2>
+
−
<p>We developed a mathematical epidemic model to predict the population dynamics of cockroaches infected by our Metarhizium anisopliae. We then performed numerical simulations on the model and sensitivity analysis on some key parameters. By doing so,
+
−
<b>
+
−
we expected to estimate the efficiency with lethal time after using our ‘cockroach terminator’ indoors, and analyzed the impact of each relative factors that influence our product so as to modify our design</b>.</p>
+
−
 
+
−
</div>
+
−
 
+
−
<div class="row">
+
−
<h2>Natural condition</h2>
+
−
<p>In natural condition indoors, due to environmental resistance like food, water and space, the population of cockroaches is more likely to follow a S-shaped growth curve (sigmoid growth curve), which can be formalized mathematically by logistic function dN/dt=rN(1-N/K).After we research papers and government researches, we found that
+
−
 
+
−
</p>
+
−
<div class="card border-success" style="border-color: #469789;">
+
−
The average number of cockroach in Shenzhen family house is 38 The grow rate (r) is 0.61, And by assumption, we set the carrying capacity (K) to be 40.
+
−
 
+
−
</div>
+
−
<p>After simulating the numerical solution (Figure 1), we can see that the cockroach population growing rapidly and reach to stationary state of 40 less than fifteen days. It’s necessary for as to control this worrisome infection.
+
−
 
+
−
</p>
+
−
<div class="col-6 offset-3">
+
−
<div class="card">
+
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/a/a0/T--SZU-China--Model_E_1.png" />
+
−
</div>
+
−
<p class="text-center" style="color: #A29F9F;">Figure 1: The population growth curve of cockroaches in house from one to large amount.</p>
+
−
</div>
+
−
</div>
+
−
<div class="row">
+
−
<h2>After control</h2>
+
−
<p>After we using M.anisopliae emulsifiable powder to control the cockroach population indoors, M.anisopliae will behave like a disease and spread out inside the population and finally kill the cockroaches. Our SID epidemic model with ordinary differential equation (ODE) was constructed based on SIR (Susceptible, Infectious, Recovered) epidemic model, following are some basic properties and assumptions of our model.
+
−
</p>
+
−
 
+
−
</div>
+
−
 
+
−
<!--Assumptions-->
+
−
<div class="indent">
+
−
<h2 id="Assumptions">Assumptions</h2>
+
−
<ul>
+
−
<li>1. The number of cockroach has reached the highest value in stable stage</li>
+
−
<li>2. Ignore natural birth and death rates in our system</li>
+
−
<li>3. Infectious individuals can not recover</li>
+
−
<li>4. Other factors that may affect the experiment are ignored</li>
+
−
</ul>
+
−
 
+
−
</div>
+
−
 
+
−
<hr>
+
−
<!--With infection-->
+
−
<div class="indent">
+
−
<h2 id="With infection">With infection</h2></div>
+
−
<div class="row">
+
−
 
+
 
<div class="col-8 offset-2">
 
<div class="col-8 offset-2">
−
<div class="card">
+
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/d/d8/T--SZU-China--HP_Model_6.png" />
+
<p style="font-size: 18px; font-weight: bold;">
−
</div>
+
On this page, two mathematical models were constructed to analysis the efficiency of our producet based on simulation and experiment. <br />Click the frame  to see the details.
−
<p class="text-center">Figure 2: Diagram of epidemic model. Showing the relation between each group of individual.</p>
+
</p>
−
 
+
−
</div>
+
−
 
+
−
<p>Our model was constructed based on SIR epidemic model (Susceptible, Infectious, Recovered) , following are some basic properties:</p>
+
−
<ul>
+
−
<li>1.The population growth of cockroach follows sigmoid growth curve.</li>
+
−
<li>2.The population of cockroach has reached close to the stationary stage.</li>
+
−
<li>3.Ignore natural death rate in our system since stubborn vitality of cockroach.</li>
+
−
<li>4.Infectious individuals can not recover, and ignore they give birth.
+
−
</li>
+
−
<li>5.Naturally all cockroaches are susceptible individuals, they can infect by M.anisopliae becoming infectious individuals.
+
−
</li>
+
−
<li>
+
−
6.The number of individual being infected in a contact between a susceptible and an infectious subject is determined by standard incidence rate βIS/(S+I).
+
−
 
+
−
</li>
+
−
<li>7.The transition rate between infectious and dead is γ, its reciprocal (1/γ) determines the average infectious period, which is estimate by experiment data.
+
−
</li>
+
−
<li>
+
−
8.Other factors that may affect the experiment are negligible.
+
−
</li>
+
−
</ul>
+
−
 
+
 
</div>
 
</div>
−
 
+
</div>
−
<hr>
+
<div class="row">
−
<!--Parameters-->
+
 
<div class="row">
 
<div class="row">
−
<h2 id="Parameters">Parameters</h2>
+
<div class="col-4 offset-1">
−
<table class="table table-bordered">
+
−
<thead>
+
<div id="card1" class="card border-3 h-100">
−
<tr class="table-active text-center">
+
<div class="view text-center">
−
<th>Parameter</th>
+
<img id="icon" class="card-img-top" style="width: 96px;" src="https://static.igem.org/mediawiki/2018/4/48/T--SZU-China--Model_home1.png"/>
−
<th>Value</th>
+
</div>
−
<th>Meaning</th>
+
<div class="card-body text-center">
−
<th>Source</th>
+
<a style="color: #469789;" href="https://2018.igem.org/Team:SZU-China/Epidemic_Model"><h3 >Epidemic Model </h3></a>
−
</tr>
+
<p class="card-text" >We developed a epidemic model with ordinary differential equation to predict the population dynamics of cockroaches infected by Metarhizium anisopliae. We then performed numerical simulations on the model and sensitivity analysis on some key parameters to find they impacts.
−
</thead>
+
</p>
−
<tbody>
+
</div>
−
<tr>
+
−
<td>β</td>
+
−
<td>0.775</td>
+
−
<td>transmission rate, which is the probability of getting the infection in a contact between susceptible and an infectious
+
−
</td>
+
−
<td>[2] [3]</td>
+
−
</tr>
+
−
<tr>
+
−
<td>γ</td>
+
−
<td>1/7</td>
+
−
<td>mortality, which is the the transition rate between I and D, its reciprocal (1/γ) determines the average infectious period
+
−
</td>
+
−
<td>B</td>
+
−
</tr>
+
−
<tr>
+
−
<td>r</td>
+
−
<td>0.61 d<sup>-1</sup></td>
+
−
<td>growth rate</td>
+
−
<td>[1]</td>
+
−
</tr>
+
−
<tr>
+
−
<td>K</td>
+
−
<td>40</td>
+
−
<td>carrying capacity</td>
+
−
<td>A</td>
+
−
</tr>
+
−
<tr>
+
−
<td>S(t)</td>
+
−
<td>variable</td>
+
−
<td>the number of susceptible individuals over time</td>
+
−
<td>-</td>
+
−
</tr>
+
−
<tr>
+
−
<td>I(t)</td>
+
−
<td>variable</td>
+
−
<td>the number of infectious individuals over time</td>
+
−
<td>-</td>
+
−
</tr>
+
−
<tr>
+
−
<td>D(t)</td>
+
−
<td>variable</td>
+
−
<td>the number of dead individuals over time</td>
+
−
<td>-</td>
+
−
</tr>
+
−
<tr>
+
−
<td>N(t)</td>
+
−
<td>S(t)+I(t)+D(t)</td>
+
−
<td>population size</td>
+
−
<td>-</td>
+
−
</tr>
+
−
<tr>
+
−
<td>S(0)</td>
+
−
<td>35</td>
+
−
<td>the initial number of susceptible individuals</td>
+
−
<td>A</td>
+
−
</tr>
+
−
<tr>
+
−
<td>I(0)</td>
+
−
<td>3</td>
+
−
<td>the initial number of infectious individuals</td>
+
−
<td>-</td>
+
−
</tr>
+
−
<tr>
+
−
<td>D(0)</td>
+
−
<td>0</td>
+
−
<td>the initial number of dead individuals</td>
+
−
<td>-</td>
+
−
</tr>
+
−
<tr>
+
−
<td colspan="4" class="text-left">
+
−
<sup>[]</sup>See in the reference code at the bottom<br>
+
−
<sup>A</sup>Government researches from Shenzhen Center for Disease Control and prevention<br>
+
−
<sup>B</sup>Estimate by experimental data
+
−
</td>
+
−
</tr>
+
−
</tbody>
+
−
</table>
+
−
<div class="indent">
+
−
<p class="text-center">Table 1: Parameters and variables for the model, all value are from research papers or experimental date
+
−
</p>
+
−
</div>
+
−
 
+
−
</div>
+
−
<div class="row">
+
−
<p>By the above relations between each group of number, the system based on assumption described above can be expressed by the following set of ordinary differential equations:
+
−
</p>
+
−
<div class="row"></div>
+
−
<p>This system is non-linear, and the analytic solution does not exist, but we can compute the numerical solution.
+
−
</p>
+
−
<p>We used improved Euler method to estimate the discrete solution in one step size to simulate the change of number with time in days (see the results). And Runge-Kutta methods solution for variable-step by MATLAB to constructed the sensitivity analysis.
+
−
</p>
+
−
 
+
−
</div>
+
−
<hr>
+
−
<!--Results-->
+
−
+
−
+
−
<div class="indent">
+
−
<h2 id="Results">Results</h2></div>
+
−
<p>The following curves (Figure 3) show dynamics of number change of each kinds of individuals.
+
−
</p>
+
−
<p>
+
−
We see that the infectious individuals grow fast before first 7 day, and then began to drop. The total number of cockroaches continuously going down.
+
−
 
+
−
</p>
+
−
<p>
+
−
We specify the median lethal time (LT50), which means 50% of cockroaches dead, where in this condition is about 8 day, and lethal time 80% (LT80) is about 15 day. This result is reasonable with our expectation, it proved that our product indeed control the total number of cockroaches efficiently in doors.
+
−
 
+
−
</p>
+
−
<div class="col-6 offset-3">
+
−
<div class="card">
+
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/c/cd/T--SZU-China--Model_E_3.png" />
+
−
 
+
 
</div>
 
</div>
−
<p class="text-center">
 
−
Figure 3: The epidemic model simulating change of each group of individual within 40 days.
 
  
−
</p>
 
−
</div>
 
−
 
−
 
−
<hr>
 
−
 
−
<div class="row">
 
−
<h2 id="Sensitivity Analysis">Sensitivity Analysis</h2>
 
−
<p>We use sensitivity analysis to analyze the impacts of some important parameter values (β, γ) on our model outcomes lethal time 50% and 80% (LT50 and LT80), in order to figure out what we should need to improve in future, and which parameters we need to be more certain.
 
−
</p>
 
−
<div class="col-6">
 
−
<p>We ran our model with transmission rate (β) and mortality (γ) ranging from 0 to 1 and plotted these values against the time in LT50 and LT80. What’s more, we introduced sensitivity coefficient (△t*p/t *△p ratio of the relative change of the lethal time to the relative change of the parameter) into consideration, increasing our parameters 10 percent, and find the relative change of the lethal time. The figures below show the changing tendency of lethal time with respect to each parameter and the sensitivity coefficient in ten percent up.
 
−
</p>
 
−
</div>
 
−
<div class="col-6">
 
−
<div class="card">
 
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/c/c1/T--SZU-China--Model_E_4.png" />
 
−
</div>
 
−
<p class="text-center" style="color: #A29F9F;">Figure 4: Change of the lethal time when transmission rate (β) ranging from 0 to 1. Blake dot represent our parameter used in Table 1.
 
−
</p>
 
−
</div>
 
−
 
−
</div>
 
−
 
−
<div class="row">
 
−
<div class="col-6">
 
−
<div class="card">
 
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/8/8e/T--SZU-China--Model_E_5.png" />
 
−
 
−
</div>
 
−
<p class="text-center" style="color: #A29F9F;">Figure 5: Change of the lethal time when mortality (γ) ranging from 0 to 1. Upper x-axis is 1/γ, indicating the average mortality day. Blake dot represent our parameter used in Table 1, and dash line showing when lager than 0.8 because ODE system can not find solution in our interval.
 
−
</p>
 
−
</div>
 
−
<div class="col-6">
 
−
<p>We can conclude the intuitive results after scanning either parameter from 0 to 1. For transmission rate (Figure 4), as beta increasing, both LTD50 and LT80 going down. However, the change of lethal time become more smoothly when beta closer to one. The black dot in curves indicates our simulating parameter value 0.775 in epidemic model, which can consider as a good approximation and high efficiency that we do not need to increase so much.
 
−
</p>
 
 
</div>
 
</div>
 
 
−
</div>
+
<div class="col-4 offset-2 ">
−
 
+
<div id="card2" class="card border-3 h-100">
−
<div class="row">
+
<div class="view text-center">
−
<div class="col-6">
+
<img id="icon" class="card-img-top" style="width: 96px;" src="https://static.igem.org/mediawiki/2018/4/4b/T--SZU-China--Model_home2.png"/>
−
<p>As for mortality (Figure 5), lethal time will increase dramatically when gamma gets closer to zero or one (when gamma closer to one, the ODE system can not find real solution in our interval), but change smoothly in a low region in the middle. That’s sound great for as since the black dot, indicating our simulating parameter value 1/7 in epidemic model, lines near the middle. So, if we can control to increase mortality a little bite, it can make a difference. This is just what we do in our project.
+
</div>
−
</p>
+
<div class="card-body text-center">
−
 
+
<a style="color: #469789;" href="https://2018.igem.org/Team:SZU-China/Statistic_Model"><h3  >Statistic Model</h3></a>
−
</div>
+
<p class="card-text">We constructed a statistical model for our experimental data from simulating room. It can give us a clear and scientific understanding that whether there had significant difference in migration rate, mortality and gnawing rate between each products.</p>
−
<div class="col-6">
+
−
<div class="card">
+
−
<img class="card-img-top" src="https://static.igem.org/mediawiki/2018/3/38/T--SZU-China--Model_E_6.png" />
+
  
 
</div>
 
</div>
−
<p class="text-center" style="color: #A29F9F;">Figure 6: Sensitivity coefficient of transmission rate (β) and mortality (γ) to LT50 and LT80.
 
−
</p>
 
−
 
 
</div>
 
</div>
−
<p>Moreover, the sensitivity coefficient in Figure 6 shows the coefficient is less than -0.8, which indicate a low sensitivity of beta and gamma, while gamma has higher sensitivity compared with beta. In conclusion, the parameters we estimated previously did not greatly affect the lethal time of our model, we can say our parameters are well approximate.
 
−
</p>
 
−
</div>
 
−
 
−
<div class="indent">
 
−
<h2>Discussion</h2>
 
−
<p>By our improved epidemic model, we successfully predict the population dynamics of cockroaches infected by our <i>Metarhizium anisopliae</i> indoors. Estimate the key parameters beta to be 0.775, gamma to be 1/7, can be considered as a good approximation not only we estimate them from research papers and experimental data, but also we constructed sensitivity analysis to these parameters and showed a low sensitivity. These results give us a clear relation about each factor and to direct our project design. Numerical simulations on the model showed our product kill cockroaches indoors with a lethal time 50% (LT50) about 8 day, and lethal time 80% (LT80) about 15 day. Thanks to the assumption, we can simplify our model, however in future, we will consider more factor that may influence our prediction to make our model more real.
 
−
</p>
 
−
 
−
</div>
 
−
<div class="indent">
 
−
<h2>References</h2>
 
−
<div class="row">
 
−
<p>
 
−
[1] D. M. Müller-Graf, Jobet, E., Cloarec, A., Rivault, C., Baalen, M. V., &#38; Morand, S. (2001). Population dynamics of host-parasite interactions in a cockroach-oxyuroid system. Oikos, 95(3), 431–440.
 
−
</p>
 
−
<p>
 
−
[2] Campaña, Ana &#38; Funderburk, Karen &#38; Kaur, Amandeep &#38; Puente, Patricia &#38; Ríos-Soto, Karen. (2017). A Household Model of Cockroach Infestation and Its Effects on Atopic Asthma Symptoms.
 
−
 
−
</p>
 
−
<p>
 
−
[3] Ni, Yu-Ting, et al. (2014). Metarhizium anisopliae as a Potential Microbial Agent for Managing the Brown-banded Cockroach, Supella longipalpa (F.).
 
−
 
−
</p>
 
−
<p>
 
−
[4] Hu, Z., Liu, S., &#38; Wang, H. (2008). Backward bifurcation of an epidemic model with standard incidence rate and treatment rate. Nonlinear Analysis Real World Applications, 9(5), 2302-2312.
 
−
</p>
 
 
</div>
 
</div>
 
</div>
 
</div>
 +
</div>
  
−
</div>
+
 +
  
−
<div class="col-2">
 
  
−
<nav id="myScrollspy" class="navbar navbar-light bg-light position-fixed">
+
</div>
−
<nav class="nav nav-pills flex-column">
+
  
−
<a class="nav-link" href="#Assumptions">Assumptions</a>
+
</div>
−
<a class="nav-link" href="#Natural condition">Natural condition</a>
+
−
<a class="nav-link" href="#With infection">With infection</a>
+
−
<a class="nav-link" href="#Parameters">Parameters</a>
+
−
<a class="nav-link" href="#Results">Results</a>
+
−
<a class="nav-link" href="#Sensitivity Analysis">Sensitivity Analysis</a>
+
−
</nav>
+
−
</nav>
+
  
−
</div>
+
</body>
−
</div>
+
−
<script>
+
−
$('body').scrollspy({
+
−
target: '#myScrollspy',
+
−
offset: 10
+
−
});
+
 +
<script>
 +
 +
$("#card1").hover(
 +
    function(){
 +
        $("#card1").addClass("shadow-lg rounded");
 +
    },
 +
    function(){
 +
        $("#card1").removeClass("shadow-lg rounded");
 +
    }
 +
);
  
−
var scroll = new SmoothScroll('a[href*="#"]');
+
$("#card2").hover(
−
</script>
+
    function(){
 +
        $("#card2").addClass("shadow-lg rounded");
 +
    },
 +
    function(){
 +
        $("#card2").removeClass("shadow-lg rounded");
 +
    }
 +
);
  
−
<style>
 
−
h1 {
 
−
color: #469789;
 
−
font-weight: 500;
 
−
}
 
−
 
−
.indent {
 
−
padding: 33px;
 
−
}
 
−
</style>
 
−
</div>
 
−
 
−
</body>
 
  
 +
$("#icon").hover(
 +
    function(){
 +
        $("#icon").addClass("bounce ");
 +
    },
 +
    function(){
 +
        $("#icon").removeClass("bounce");
 +
    }
 +
);
 +
</script>
 +
  
  
 
</html>
 
</html>

Latest revision as of 18:33, 17 October 2018

Model

On this page, two mathematical models were constructed to analysis the efficiency of our producet based on simulation and experiment.
Click the frame to see the details.

Epidemic Model

We developed a epidemic model with ordinary differential equation to predict the population dynamics of cockroaches infected by Metarhizium anisopliae. We then performed numerical simulations on the model and sensitivity analysis on some key parameters to find they impacts.

Statistic Model

We constructed a statistical model for our experimental data from simulating room. It can give us a clear and scientific understanding that whether there had significant difference in migration rate, mortality and gnawing rate between each products.