|
|
| Line 1: |
Line 1: |
| | {{SZU-China/Nav}} | | {{SZU-China/Nav}} |
| | <html> | | <html> |
| − |
| |
| − | <style>
| |
| − | h1 {
| |
| − | color: #3880A0;
| |
| − | }
| |
| − |
| |
| − | h2 h3 h4 {
| |
| − | font-weight: bold;
| |
| − | }
| |
| − |
| |
| − | .row {
| |
| − | padding: 48px;
| |
| − | }
| |
| − |
| |
| − | body {
| |
| − | 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>
| |
| − |
| |
| − | <script src="https://2018.igem.org/common/MathJax-2.5-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
| |
| − |
| |
| − | <script type="text/x-mathjax-config">
| |
| − | MathJax.Hub.Config({
| |
| − | tex2jax: {inlineMath: [['$','$']]}
| |
| − | });
| |
| − | </script>
| |
| | | | |
| | | | |
| | + | <style> |
| | + | #icon{ |
| | + | transform:translateY(-49px); |
| | + | } |
| | | | |
| | + | .border-3 { |
| | + | border-width:3px !important; |
| | + | } |
| | | | |
| | + | </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="card">
| |
| − | <img class="card-img-top" src="https://static.igem.org/mediawiki/2018/e/ea/T--SZU-China--model1.jpg" />
| |
| − | </div>
| |
| − | </div>
| |
| | <div class="container-fluid"> | | <div class="container-fluid"> |
| − | <div class="row" style="background-image: url();"> | + | <div class="row"> |
| − | <div class="col-9 offset-1 shadow"> | + | <div class="col-9 offset-1 shadow "> |
| − | <div class="text-center">
| + | |
| − | <h1 id="header">Epidemic Model</h1>
| + | |
| − | </div>
| + | |
| − | | + | |
| | <div class="row"> | | <div class="row"> |
| − | | + | <div class="col-5 offset-1"> |
| − | <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,
| + | <div id=" card1" class=" card border- 3"> |
| − | <b>
| + | <div class=" view text- center"> |
| − | 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>
| + | <img id="icon" class="card-img-top" style="width: 96px;" src="https://static.igem.org/mediawiki/2018/0/05/T--SZU-China--Model1.png"/> |
| − | | + | </div> |
| − | </div>
| + | <div class="card-body text-center"> |
| − | | + | <a href="https://2018.igem.org/Team:SZU-China/Epidemic_Model"><h3>Epidemic Model </h3></a> |
| − | | + | <p>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. |
| − | <h2>Natural condition</h2>
| + | </p> |
| − | <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 $\dfrac {dN}{dt}=rN\left( 1-\dfrac {N}{K}\right)$ .After we research papers and government researches, we found that
| + | </div> |
| − | | + | |
| − | </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> | | </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="card">
| |
| − | <img class="card-img-top" src="https://static.igem.org/mediawiki/2018/d/d8/T--SZU-China--HP_Model_6.png" />
| |
| − | </div>
| |
| − | <p class="text-center">Figure 2: Diagram of epidemic model. Showing the relation between each group of individual.</p>
| |
| | | | |
| | </div> | | </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 $\dfrac {\beta IS}{S+I}$ .
| |
| − |
| |
| − | </li>
| |
| − | <li>7.The transition rate between infectious and dead is γ, its reciprocal ($\dfrac {1}{\gamma}$) 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>
| |
| − |
| |
| − | <hr>
| |
| − | <!--Parameters-->
| |
| − | <div class="row">
| |
| − | <h2 id="Parameters">Parameters</h2>
| |
| − | <table class="table table-bordered">
| |
| − | <thead>
| |
| − | <tr class="table-active text-center">
| |
| − | <th>Parameter</th>
| |
| − | <th>Value</th>
| |
| − | <th>Meaning</th>
| |
| − | <th>Source</th>
| |
| − | </tr>
| |
| − | </thead>
| |
| − | <tbody>
| |
| − | <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:
| |
| | | | |
| − | | + | <div class="col-5 offset-1"> |
| − | $$\dfrac {dS}{dt}=rN\left( 1-\dfrac {N}{K}\right) -\dfrac {\beta IS}{S+I}$$ $$\dfrac {dI}{dt}=\dfrac {\beta IS}{S+I}-\gamma I$$ $$\dfrac {dD}{dt}=\gamma I$$ $$N=S+I+D$$ | + | <div id="card2" class="card border-3"> |
| | + | <div class="view text-center"> |
| | + | <img id="icon" class="card-img-top" style="width: 96px;" src="https://static.igem.org/mediawiki/2018/5/53/T--SZU-China--Model2.png"/> |
| | + | </div> |
| | + | <div class="card-body text-center"> |
| | + | <a href="https://2018.igem.org/Team:SZU-China/Statistic_Model"><h3>Statistical Model</h3></a> |
| | + | <p>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> |
| | | | |
| − |
| |
| − | 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>
| |
| − | <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>
| |
| − |
| |
| − |
| |
| − |
| |
| − | <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>
| |
| − | <p class="text-center" style="color: #A29F9F;">
| |
| − | 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> | | </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>
| |
| − |
| |
| − | <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> |
| − | <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="row">
| + | |
| − | | + | </div> |
| − | <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.
| + | |
| − | </p>
| + | |
| | | | |
| − | </div>
| |
| − | <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">
| + | </body> |
| − | <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.
| + | <script> |
| − | </p>
| + | |
| | + | $("#card1").hover( |
| | + | function(){ |
| | + | $("#card1").addClass("shadow-lg rounded"); |
| | + | }, |
| | + | function(){ |
| | + | $("#card1").removeClass("shadow-lg rounded"); |
| | + | } |
| | + | ); |
| | | | |
| − | </div>
| + | $("#card2").hover( |
| − | <div class="indent">
| + | function(){ |
| − | <h2>References</h2>
| + | $("#card2").addClass("shadow-lg rounded"); |
| − | <div class="row">
| + | }, |
| − | <p>
| + | function(){ |
| − | [1] D. M. Müller-Graf, Jobet, E., Cloarec, A., Rivault, C., Baalen, M. V., & Morand, S. (2001). Population dynamics of host-parasite interactions in a cockroach-oxyuroid system. Oikos, 95(3), 431–440.
| + | $("#card2").removeClass("shadow-lg rounded"); |
| − | </p>
| + | } |
| − | <p>
| + | ); |
| − | [2] Campaña, Ana & Funderburk, Karen & Kaur, Amandeep & Puente, Patricia & 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>
| + | $("#icon").hover( |
| − | <p>
| + | function(){ |
| − | [4] Hu, Z., Liu, S., & 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.
| + | $("#icon").addClass("bounce "); |
| − | </p>
| + | }, |
| − | </div>
| + | function(){ |
| − | </div>
| + | $("#icon").removeClass("bounce"); |
| − |
| + | } |
| − |
| + | ); |
| − |
| + | </script> |
| − |
| + | |
| − |
| + | |
| − |
| + | |
| − |
| + | |
| − | <!--Model2-->
| + | |
| − |
| + | |
| − | <div class="indent">
| + | |
| − | <h1 style="color: #469789;">Statistic Model</h1>
| + | |
| − | </div>
| + | |
| − | | + | |
| − | <div class="row">
| + | |
| − | <h2>Introduction</h2>
| + | |
| − | <p>
| + | |
| − | Statistic analysis can give us a clear and scientific understanding of our complex data, that was the reason why we construct a statistical model for our experimental data.
| + | |
| − | | + | |
| − | </p>
| + | |
| − | <p>
| + | |
| − | The data from simulating room were subjected to statistic analysis to determine whether there are significantly different in migration rate, mortality and gnawing rate using
| + | |
| − | <i>One-way ANOVA</i>.
| + | |
| − | <font color="#469789">Levene’s test</font> to determine whether each group is homogeneity, then
| + | |
| − | <font color="#469789">T-test</font> for pooled data in pairwise gave us the value and finally to evaluate significant difference.
| + | |
| − | | + | |
| − | </p>
| + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <h2>Levene’s test</h2>
| + | |
| − | <p>Levene’s test is used to assess the equality of variances for a variable from two or more groups. In our model, we have four independent groups. Since in the following procedures for T-test, we assume that the variances of the populations from four samples are equal, so we need to test whether the following variances can be assume to be equal.</p>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <h4>Assumptions:</h4>
| + | |
| − | <p>
| + | |
| − | <font color="#469789">1. Due to the advantage of Levene’s test, we do not need to assume the normal data
| + | |
| − | </font>
| + | |
| − | </p>
| + | |
| − | <p>
| + | |
| − | <font color="#469789">2. Null hypothesis that all populations from samples have the same variances
| + | |
| − | </font>
| + | |
| − | $$ H_{0}:\sigma _{1}=\sigma _{2}=\cdots \sigma _{k}, \quad H_{A}:\sigma _{i}\neq \sigma _{j}, \quad \alpha =0.05 $$
| + | |
| − | </p>
| + | |
| − | <p>The Levene Statistic (W) of each groups is calculated by: $$ W{=}\dfrac {\left( N-k\right) }{\left( k-1\right) } \dfrac {\sum ^{k}_{i=1}N_{i}\left( Z_{i\cdot}-Z_{i\cdot\cdot}\right) ^{2}} {\sum ^{k}_{i=1}\sum ^{N_{i}}_{j=1}\left( Z_{ij}-Z_{i\cdot}\right) ^{2}} \sim F\left( \alpha ,k-1,N-k\right) $$
| + | |
| − | </p>
| + | |
| − | <p>Where</p>
| + | |
| − | | + | |
| − | <ul>
| + | |
| − | <li> $k$ is the number of groups which need to be tested, equal to 4 in our data
| + | |
| − | </li>
| + | |
| − | <li>$N_{i}$ is the number of cases in the ith group</li>
| + | |
| − | <li>$N$ is the total number of cases in all groups, equal to 12 in our data</li>
| + | |
| − | <li>$Y_{ij}$ is the value of the measured variable for the $i$th case from the $j$th group</li>
| + | |
| − | <li>$Z_{ij}=\left| Y_{ij}-\overline {Y_{i}}\right|$, $\overline {Y_{i}}$ is the mean of the $i$th group</li>
| + | |
| − | <li>$\overline {Z_{i\cdot}}$ is the mean of the $Z_{ij}$ for group $i$</li>
| + | |
| − | <li>$\overline {Z_{\cdot\cdot}}$ is the mean of all data</li>
| + | |
| − | </ul>
| + | |
| − | <div class="indent">
| + | |
| − | <p>After calculate Levene Statistic of each group with SPSS, the W value are shown below:</p>
| + | |
| − | <table class="table">
| + | |
| − | <caption>Test of Homogeneity of Variances</caption>
| + | |
| − | <thead>
| + | |
| − | <tr>
| + | |
| − | <th scope="col"></th>
| + | |
| − | <th scope="col">Levene Statistic</th>
| + | |
| − | <th scope="col">df1</th>
| + | |
| − | <th scope="col">df2</th>
| + | |
| − | <th scope="col">Sig.</th>
| + | |
| − | </tr>
| + | |
| − | </thead>
| + | |
| − | <tbody>
| + | |
| − | <tr>
| + | |
| − | <th scope="">Migration Rate</th>
| + | |
| − | <td>2.070</td>
| + | |
| − | <td>3</td>
| + | |
| − | <td>8</td>
| + | |
| − | <td>.183</td>
| + | |
| − | </tr>
| + | |
| − | <tr>
| + | |
| − | <th scope="">Mortality</th>
| + | |
| − | <td>3.122</td>
| + | |
| − | <td>3</td>
| + | |
| − | <td>8</td>
| + | |
| − | <td>.088</td>
| + | |
| − | </tr>
| + | |
| − | <tr>
| + | |
| − | <th scope="">Gnawing Rate</th>
| + | |
| − | <td>.000</td>
| + | |
| − | <td>3</td>
| + | |
| − | <td>3</td>
| + | |
| − | <td>1.000</td>
| + | |
| − | </tr>
| + | |
| − | </tbody>
| + | |
| − | </table>
| + | |
| − | </div>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | | + | |
| − | <div class="indent">
| + | |
| − | <p>
| + | |
| − | Fortunately, all
| + | |
| − | <font color="#469789">$ W < F\left( 0.05,3,8\right) =4.066 $ or Sig.>0.05</font>, indicating the null hypothesis is accepted, meaning
| + | |
| − | <font color="#469789">$ \sigma _{1},\sigma _{2},\sigma _{3},\sigma _{4} $,</font>
| + | |
| − | do not have statistically different, so we can say they have the same variance σ2, and we can continues our procedure.
| + | |
| − | | + | |
| − | </p>
| + | |
| − | </div>
| + | |
| − | | + | |
| − | <div class="indent">
| + | |
| − | <h2>T-test for pooled data</h2>
| + | |
| − | <p>The independent T-test for pooled data is used to determine whether the mean difference between two groups is statistically significantly different
| + | |
| − | </p>
| + | |
| − | <h5>
| + | |
| − | Assumptions:
| + | |
| − | </h5>
| + | |
| − | <p style="color: #469789;">1.The sample groups are pairwise independent</p>
| + | |
| − | <p style="color: #469789;">2.All populations from samples have the same variances</p>
| + | |
| − | <p>The T-value (t) of eachpair is calculated by:$$t_{n_{1}+n_{2}-2}=\dfrac { \overline {x_{1}}-\overline {x_{2}}}
| + | |
| − | {\sqrt {\dfrac {\left( n_{1}-1\right) s^{2}_{1}+\left( n_{2}-1\right) s^{2}_{2}}
| + | |
| − | {\left( n_{1}-1\right) + \left( n_{2}-1\right) }\left( \dfrac {1}{n_{1}}+\dfrac {1}{n_{2}}\right) }}$$ </p>
| + | |
| − |
| + | |
| − | | + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <p>When $ n_{1}=n_{2}=n=3 $ in this case, it can be simplify as:
| + | |
| − | $$ t_{2\left( n-1\right) }=\dfrac {\overline {x_{1}}-\overline {x_{2}}}{\sqrt {\dfrac {s^{2}_{1}+s^{2}_{2}}{n}}} $$
| + | |
| − | </p>
| + | |
| − |
| + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <p>By using two-tailed test, if p-value is lower than 0.05 or $ t > t_{4,0.05\left( two-tailed \right) }=2.776 $, we can say two mean have statistically significantly different.
| + | |
| − | </p>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | | + | |
| − | <div class="indent">
| + | |
| − | <h3>Results</h3>
| + | |
| − | <p>Levene’s test, gave us results that each group is homogeneity at P-value of 0.05, indicating we can assume they have equal variances.T-test for pooled data gave us following results:
| + | |
| − | </p>
| + | |
| − | <p>The graphs in figure 1 divided in matrix cell Zij show the visualized results of above procedures, t value lager than 2.776 showing red means Zi statistically significantly higher than Zj, while lower than -2.776 showing green means Zi statistically significantly lower than Zj. More red or green means different more significantly. Where, Z1, Z2, Z3, Z4 represent cockroach killing chalk, BAYER-Premise, M.anisopliae emulsifiable powder and blank group respectively.
| + | |
| − | </p>
| + | |
| − | <div class="row">
| + | |
| − | <div class="col-5 offset-1">
| + | |
| − | <div class="card h-100">
| + | |
| − | <img class="card-img-top" src="https://static.igem.org/mediawiki/2018/0/03/T--SZU-China--S_1.png" />
| + | |
| − | </div>
| + | |
| − | </div>
| + | |
| − | <div class="col-5 offset-1">
| + | |
| − | <div class="card h-100">
| + | |
| − | <img class="card-img-top" src="https://static.igem.org/mediawiki/2018/1/17/T--SZU-China--S_2.png" />
| + | |
| − | </div>
| + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <div class="col-6 offset-3">
| + | |
| − | <div class="card">
| + | |
| − | <img class="card-img-top" src="https://static.igem.org/mediawiki/2018/e/e6/T--SZU-China--S_3.png" />
| + | |
| − | </div>
| + | |
| − | <p class="text-center" style="color: #A29F9F;">Figure 1: Results of independent T-test for migration rate, mortality and gnawing rate. Matrix cell Zij means Zi compared to Zj. More red or green means different more significantly.
| + | |
| − | </p>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | </div>
| + | |
| − | </div>
| + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <p>We can obviously see that cockroaches will behave strong migration tendency after people using cockroach killing chalk or BAYER-Premise. Migration rate infected by M.anisopliae emulsifiable powder was significantly lower than that by cockroach killing chalk or BAYER-Premise (P
| + | |
| − | <0.05), even lower than blank group. While, there was no significant difference in migration rate between cockroach killing chalk and BAYER-Premise (P>0.05).
| + | |
| − | </p>
| + | |
| − | <p>
| + | |
| − | Similarly, gnawing rate in M.anisopliae group was significantly higher than that in cockroach killing chalk or BAYER-Premise group (P
| + | |
| − | <0.05), while between the latter two groups, there was no significant difference(P>0.05). Thus, M.anisopliae will spread easier in cockroaches population causing an epidemic-like disease. As for mortality within three days, there was no significant difference between each two groups (P>0.05).
| + | |
| − | | + | |
| − | </p>
| + | |
| − | <p>
| + | |
| − | The results above show that M.anisopliae emulsifiable powder has higher efficiency than other traditional cockroach-killing methods. It has significantly lower migration rate, higher gnawing rate to spread out naturally, and has similar mortality after using it.
| + | |
| − | </p>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | <div class="indent">
| + | |
| − | <h3>References</h3>
| + | |
| − | <p>[1] Levene, Howard (1960). “Robust tests for equality of variances”. In Ingram Olkin; Harold Hotelling; et al. Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling. Stanford University Press. pp. 278–292.</p>
| + | |
| − | <p>[2]
| + | |
| − | <a href="https://en.wikipedia.org/wiki/Levene%27s_test">https://en.wikipedia.org/wiki/Levene%27s_test</a>
| + | |
| − | </p>
| + | |
| − | </div>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | | + | |
| − | </div>
| + | |
| − | <!--
| + | |
| − | <div class="col-2">
| + | |
| − | | + | |
| − | <nav id="myScrollspy" class="navbar navbar-light bg-light position-fixed">
| + | |
| − | <nav class="nav nav-pills flex-column">
| + | |
| − | | + | |
| − | <a class="nav-link" href="#Assumptions">Assumptions</a>
| + | |
| − | <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>
| + | |
| − |
| + | |
| − | -->
| + | |
| − | </div>
| + | |
| − | <script>
| + | |
| − | $('body').scrollspy({
| + | |
| − | target: '#myScrollspy',
| + | |
| − | offset: 10
| + | |
| − | });
| + | |
| − | | + | |
| − | var scroll = new SmoothScroll('a[href*="#"]');
| + | |
| − | </script>
| + | |
| − | | + | |
| − | <style>
| + | |
| − | h1 {
| + | |
| − | color: #469789;
| + | |
| − | font-weight: 500;
| + | |
| − | }
| + | |
| − |
| + | |
| − | .indent {
| + | |
| − | padding: 33px;
| + | |
| − | }
| + | |
| − | </style>
| + | |
| − | </div>
| + | |
| − | | + | |
| − | </body> | + | |
| | | | |
| | | | |
| | </html> | | </html> |
