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

Line 1: Line 1:
<!DOCTYPE html>
 
 
<html>
 
<html>
<head>
 
<link rel="stylesheet" type="text/css" href="https://2018.igem.org/wiki/index.php?title=Template:SZU-China/CSS&action=raw&ctype=text/css" />
 
  
<meta charset="utf-8" />
+
<head>
<title></title>
+
</head>
+
<body>
+
+
+
         
+
         
+
         
+
+
          <!--<ul class="nav navbar-nav ml-auto">
+
            <li class="nav-item">
+
              <a class="nav-link" href="http://builtwithbootstrap.com/" target="_blank">Built With Bootstrap</a>
+
            </li>
+
            <li class="nav-item">
+
              <a class="nav-link" href="https://wrapbootstrap.com/?ref=bsw" target="_blank">WrapBootstrap</a>
+
            </li>
+
          </ul>-->
+
  
 +
    <style>
 +
        #sideMenu, #top_title {
 +
            display: none;
 +
        }
 +
 +
        #content {
 +
            width: 100%;
 +
            margin: 0px;
 +
            padding: 0px;
 +
            background: #FFF;
 +
        }
 +
 +
        .def li{
 +
            font-style:italic;
 +
            text-decoration:none;
 +
        }
 +
 
 +
 +
   
 +
    </style>
 +
 +
    <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
 +
    <link rel="stylesheet" type="text/css" href="https://2017.igem.org/Template:SZU-China/sub/CSS?action=raw&amp;ctype=text/css">
 +
    <link href="https://2017.igem.org/Template:SZU-China/HP/CSS?action=raw&amp;ctype=text/css" rel="stylesheet" type="text/css">
 +
    <link href='http://fonts.font.im/css?family=Open+Sans' rel='stylesheet' type='text/css'>
 +
    <script type="text/javascript" src="https://2017.igem.org/Template:SZU-China/jquery/Javascript?action=raw&amp;ctype=text/javascript"></script>
 +
 +
    <script type="text/javascript">
 +
        $(document).ready(function () {
 +
            $("#nav-list-ul li").hover(function () {
 +
                $(this).find('ul.sub-nav').show();
 +
            }, function () {
 +
                $(this).find('ul.sub-nav').hide();
 +
            });
 +
 +
            $("div.navbar").hover(function () {
 +
                $(this).css("box-shadow", "1px 1px 7px #000")
 +
            }, function () {
 +
                $(this).css("box-shadow", "0px 0px 0px #000")
 +
            });
 +
 +
            $("#top").click(function () {
 +
                $('body,html').animate({ scrollTop: 0 }, 200);
 +
                return false;
 +
            });
 +
        });
 +
    </script>
 +
 +
</head>
 +
 +
<body style="background: #FEFEFE; font-family: "open sans", "Helvetica Neue" , helvetica, arial, sans-serif;">
 +
 +
    <div class="navbar" style="top:5px;position:fixed;z-index:5">
 +
        <div class="navbar-nav">
 +
            <ul id="nav-list-ul" class="nav-list">
 +
 +
                <li>
 +
                    <a href="#">TEAM</a>
 +
                    <span class="caret"></span>
 +
                    <ul class="sub-nav" style="top:-5px">
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Team">MEMBERS</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Collaborations">COLLABORATIONS</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Attributions">ATTRIBUTION</a></li>
 +
                    </ul>
 +
                </li>
 +
 +
 +
                <li>
 +
                    <a href="#">PRACTICE</a>
 +
                    <span class="caret"></span>
 +
                    <ul class="sub-nav" style="top:-5px">
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/HP/Gold_Integrated">INTEGRATED HP</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/HP/Silver">SILVER HP</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Engagement">ENGAGEMENT</a></li>
 +
                    </ul>
 +
                </li>
 +
 +
 +
                <li>
 +
                    <a href="#">EXPERIMENT</a>
 +
                    <span class="caret"></span>
 +
                    <ul class="sub-nav" style="top:-5px">
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Procedure">PROCEDURE</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Results">RESULTS</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Demonstrate">DEMONSTRATE</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Protocol">PROTOCOL</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Notebook">NOTEBOOK</a></li>
 +
 +
                    </ul>
 +
                </li>
 +
                <li>
 +
                    <a href="#">PROJECT</a>
 +
                    <span class="caret"></span>
 +
                    <ul class="sub-nav" style="top:-5px">
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Description">DESCRIPTION</a></li>
 +
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Design">DESIGN</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Model">MODEL</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Safety">SAFTY</a></li>
 +
                    </ul>
 +
                </li>
 +
 +
 +
                <li>
 +
                    <a href="#">ACHIEVEMENT</a>
 +
                    <span class="caret"></span>
 +
                    <ul class="sub-nav" style="top:-5px">
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Achievement">MEDAL</a></li>
 +
                        <li><a href="https://2017.igem.org/Team:SZU-China/Parts">PARTS</a></li>
 +
                    </ul>
 +
                </li>
 +
            </ul>
 +
        </div>
 +
        <div>
 +
            <a href="https://2017.igem.org/Team:SZU-China"><img src="https://static.igem.org/mediawiki/2017/1/17/T--SZU-China--team-logo.png" style="height: 36px;width: auto;padding-top: 12px;padding-left:10px;"></a>
 
         </div>
 
         </div>
      </div>
 
 
     </div>
 
     </div>
 
+
 
 
+
    <div style="position: relative;top: 80px;">
 
+
        <section id="overview" style="padding:96px 0 40px 0;background-color:white;">
 
+
            <div class="container">
 
+
                <div class="row">
  <div class="row">
+
                    <div class="col-sm-12">
<div class="col-9 offset-1">
+
                        <h3 class="uppercase color-primary mb40 " style="margin-bottom: 40px;font-size:50px"><center>MODEL</center> </h3>
<div class="bs-docs-section">
+
 
  <div data-spy="scroll" data-target="#list-example" data-offset="0" class="scrollspy-example">
+
 
<div>
+
                    </div>
<div class="bs-component">
+
                </div>
<div class="jumbotron">
+
            </div>
<div class="container"></div>
+
        </section>
<div class="center-block">
+
 
<h1> </h1>
+
 
<p> </p>
+
        <p></p><br /><br />
<h1 id="header" >Model</h1>
+
 
<p>We set up a mathematical model to predict the population dynamics of cockroaches before and after using our product. By doing so, we can estimate the lethal time of our cockroaches terminator, analyse the relationships among each relative factors so as to modify our product.</p>
+
        <section style="background-color: rgba(245,245,245,0.45); padding: 96px 0; ">
</div>
+
            <div style="position: relative; margin: 30px 15%;">
</div>
+
                <center><p class="font1" style="font-weight: 500; text-align: center; font-size: 25px; color: rgb(75, 151, 165); ">Introduction</p></center><br /><br />
</div>
+
                <p style="font-size:16px">This year, our team creates a mathematical representation of our concrete self-healing system. This representation, or model, constructs a judging scale with which we can utilize to regulate the four main environmental factors affecting our final concrete healing efficiency (reflected by mineralization activity). </p><br /><br />
+
                <p style="font-size:16px">
<div>
+
                    The four factors are: <br /><img src="https://static.igem.org/mediawiki/2017/c/c6/T--SZU-China--equation.jpg " width="25%" /></br>
<h2 id="Assumptions">Assumptions</h2>
+
                </p>
<ul>
+
                <p style="font-size:16px">Based on this model, we can also design the best ‘package'– the microcapsule shells with adequate nutrition combination.</p>
<li>1. The number of cockroach has reached the highest value in stable stage</li>
+
                <br />
<li>2. Ignore natural birth and death rates in our system</li>
+
               
<li>3. Infectious individuals can not recover</li>
+
                <p style="font-size:16px">This modelling process, presented below, can be seen as a feedback from the wet lab (experiment result) and the dry lab(modelling analysis).</p><br /><br />
<li>4. Other factors that may affect the experiment are ignored</li>
+
 
</ul>
+
                <p style="font-size:16px">The following page shows how we conducted modelling approaches to achieve our goals. To begin with, we make some definition.</p><br /><br />
<h2 id="Natural condition">Natural condition</h2>
+
                <br /><br />
<p>In natural condition indoors, due to environmental resistance like food, water, 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.</p>
+
                <center><span class="font1" style="font-weight: 500; text-align: center; font-size: 25px; color: rgb(75, 151, 165); ">Definition</span></center><br /><br />
</div>
+
               
<div class="text-center">
+
                <br/><br/>
<img class="rounded" style="width: 420px;" src="../img/Model_1.png"/>
+
                <ul style="font-size:16px;" class="def">Variables and nomenclature:<li>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;concentration of spores - c[Spore]</li> <li>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;concentration of carbon source(C3H5O3Na) - c[C3H5O3Na]</li><li> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;concentration of nitrogen source(NaNO3) - c[NaNO3]</li><li>&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;pH of the media - pH</li></ul><br />
</div>
+
                <ul style="font-size: 16px;" class="def">Goal of model:<li>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Use preliminary data to guide future experiments.</li></ul>
+
                <br /><br />
<div>
+
                <center><span class="font1" style="font-weight: 500; text-align: center; font-size: 25px; color: rgb(75, 151, 165); ">Procedure</span></center><br /><br />
<h2 id="With infection">With infection</h2>
+
                <p style="font-size:16px">
<p>Our model was constructed based on SIR epidemic model (Susceptible, Infectious, Recovered) , following are some basic properties:</p>
+
                    1.&nbsp;&nbsp;Standardizing variables: transforming variables into the same scale. Here we utilize the z-score standardizing method, such that
<ul>
+
                </p>
<li>1. Naturally all cockroaches are susceptible individuals, they can infect by M.anisopliae becoming infectious individuals.</li>
+
                <div style="text-align:center"><img src="https://static.igem.org/mediawiki/2017/0/07/T--SZU-China--standardization.jpg" width="30%" ><br /></div><br />
<li>2. The number of individual being infected in a contact between a susceptible and an infectious subject is simulate by standard incidence      .</li>
+
                <p style="font-size:16px">
<li>3. The transition rate between Infectious and dead is γ, its reciprocal (1/γ) determines the average infectious period, which is estimate by experiment data.</li>
+
                    where the symbols in the formula are variable, standardized variable, sample mean, andstandard deviation respectively.
+
                </p>
</ul>
+
                <p style="font-size:16px">2.&nbsp;&nbsp;Fitting functions of each variable with polynomial function.</p>
<div class="text-center">
+
                <p style="font-size:16px">3.&nbsp;&nbsp;Getting the overall relationship using linear least square method.</p><br /><br /><br />
<img width="560px" src="../img/Model_2.png"/>
+
                <center> <span class="font1" style="font-weight: 500; text-align: center; font-size: 25px; color:rgb(75, 151, 165); ">Results</span></center><br /><br />
</div>
+
 
+
 
+
                <div style="text-align: center;margin: 24px 0;z-index:-1;">
</div>
+
                    <img src="https://static.igem.org/mediawiki/2017/2/24/T--SZU-China--model1.png " width="90%" style="box-shadow: 0px 0px 2px #1E1E1E;">
<div>
+
                </div>
<h2 id="Parameters">Parameters</h2>
+
 
<p>Model was simulate during 30 days, with total number of 60.</p>
+
                <p style="font-size:16px">The graph above depicts the polynomial regression of 4 factors. In general, for each of those four factors, the mineralization activity shows the similar tendency of going up first and down later with the increase of each factor.</p><br /><br />
<table class="table table-bordered">
+
 
<thead>
+
                <p style="font-size:16px">The fitted functions are showed as follows:</p>
<tr class="table-active">
+
                <div style="text-align:center"><img src="https://static.igem.org/mediawiki/2017/1/1d/T--SZU-China--fitting.jpg" width="60%" ><br /></div>
<th>Parameter</th>
+
             
<th>Value</th>
+
                <p style="font-size:16px">We then fitted x with y linearly to get the overall regression equation which can describe the weight of four variables respectively.</p><br /><br />
<th>Meaning</th>
+
 
</tr>
+
                <center><span class="font1" style="font-weight:bold;font-size:16px">Overall regression equation:</span></center><br /><br />
</thead>
+
                <div style="text-align:center"><img src="https://static.igem.org/mediawiki/2017/a/ad/T--SZU-China--regression.jpg" width="60%" ><br /></div>
<tbody class="text-center">
+
 
<tr>
+
                <p style="font-size:16px">with</p><br /><br />
<td>S(t)</td>
+
 
<td></td>
+
                <div style="text-align:center"><img src="https://static.igem.org/mediawiki/2017/b/be/T--SZU-China--regressiontransform.jpg" width="30%" ><br /></div>
<td>the number of susceptible individuals over time</td>
+
 
</tr>
+
                <p style="font-size:16px">where a<sub>i1</sub>represents the coefficient of the highest power term of curve-fitting equation, such that the coefficient of the highest power term are of the same scale.</p><br /><br />
<tr>
+
                <br/>
<td>I(t)</td>
+
 
<td></td>
+
                <p style="font-size:16px">From the equation above, we can see that nitrogen source has the maximum weight, while pH has the minimum weight, which means nitrogen source is the most essential nutrition for B.subtilis spore. Additionally, the low weight of pH shows that spores are not sensitive to the change of pH in a relatively apt range, although there is a sharp decline in activity when pH reaches 11.</p><br />
<td>the number of infectious individuals over time</td>
+
 
</tr>
+
                <br /><p style="font-size:16px">In summary, as long as we keep the environment below the boundary high-pH point, it makes no much difference how much we have improved the alkaline resistance of <i>B.subtilis</i> spore. In this way, this model instructs us on how to modify our formulation, such that the ingredients can be made full use of.</p><br />
<tr>
+
                <span></span><br />
<td>D(t)</td>
+
                <span></span><br />
<td></td>
+
                <span></span><br />
<td>the number of dead individuals over time</td>
+
                <br /><br /><br /><br /><br /><br />
</tr>
+
            </div>
<tr>
+
        </section>
<td>β</td>
+
    </div>
<td>0.75</td>
+
 
<td>transmission rate, which is the probability of getting the infection in a contact between susceptible and an infectious</td>
+
 
</tr>
+
    <div class="footer" style="position:absolute;margin-bottom:0px;padding-bottom:0">
<tr>
+
        <div class="sponsor">
<td>γ</td>
+
            <img src="https://static.igem.org/mediawiki/2017/0/08/2017SZUChina-sponsor1.jpg" class="sponsor-img">
<td>1/8</td>
+
            <img src="https://static.igem.org/mediawiki/2017/d/d4/2017SZUChina-sponsor2.jpg" class="sponsor-img">
<td>mortality, which is the the transition rate between I and D, its reciprocal (1/γ) determines the average infectious period</td>
+
            <img src="https://static.igem.org/mediawiki/2017/d/d5/2017SZUChina-sponsor3.jpg" class="sponsor-img">
</tr>
+
            <img src="https://static.igem.org/mediawiki/2017/1/1e/2017SZUChina-sponsor4.jpg" class="sponsor-img">
<tr>
+
            <img src="https://static.igem.org/mediawiki/2017/a/a3/2017SZUChina-sponsor5.jpg" class="sponsor-img">
<td>S(0)</td>
+
        </div>
<td>55</td>
+
        <div class="footer-info">
<td>the initial number of susceptible individuals</td>
+
            <address style="color:#888">3688 Nanhai Ave, ShenZhen, GuangDong Province, China</address>
</tr>
+
            <span class="copyright">Copyright © 2017.<a href="https://2017.igem.org/Team:SZU-China" class="mainpage">SZU-China</a> iGEM</span>
<tr>
+
        </div>
<td>I(0)</td>
+
    </div>
<td>5</td>
+
 
<td>the initial number of infectious individuals</td>
+
    <div class="arrow-div">
</tr>
+
        <a id="top" href="#" title="Back to top">
<tr>
+
            <img src="https://static.igem.org/mediawiki/2017/f/f6/T--SZU-China--arrow.png" class="arrow-img">
<td>r</td>
+
        </a>
<td>0.3</td>
+
    </div>
<td>growth rate</td>
+
 
</tr>
+
</body>
<tr>
+
 
<td>N=S+I</td>
+
 
<td></td>
+
 
<td>population size</td>
+
 
</tr>
+
<tr>
+
<td>K</td>
+
<td>70</td>
+
<td>carring capacity</td>
+
</tr>
+
</tbody>
+
</table>
+
<P>The system without so-called vital dynamics (birth and death) described above can be expressed by the following set of ordinary differential equations:</P>
+
<p>This system is non-linear, and the analytic solution does not exist, but we can compute the numerical solution by MATLAB. (see results)</p>
+
</div>
+
<div>
+
<h2 id="Results">Results</h2>
+
<p>The following curves show dynamics number change of each kinds of individuals. We see that the infectious individuals grow fast before first 6 day, and then began to drop. The total number of cockroaches continuously going down. We specify the median lethal time (LT50), which in this condition is 11.1 days.</p>
+
<div class="text-center">
+
<img class="rounded" style="width: 420px;" src="../img/Model_3.png"/>
+
</div>
+
</div>
+
<div>
+
<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 (LT50). The figures below show the tendency of dead number with respect to each parameter change. </p>
+
+
+
</div>
+
<div class="card-group">
+
<div class="card">
+
<div class="card-header">
+
<h3>1. change gamma</h3>
+
</div>
+
+
<img class="card-img-top" style="width: 512px;" src="../img/Model_4.png" />
+
<div class="card card-body">
+
<table >
+
<thead class="table table-active">
+
<tr>
+
<th>change α</th>
+
<th>α</th>
+
<th>β</th>
+
<th>LT50</th>
+
<th>ΔLT50</th>
+
<th>Ratio</th>
+
</tr>
+
</thead>
+
<tbody>
+
<tr>
+
<td>+20%</td>
+
<td>0.750</td>
+
<td>0.150</td>
+
<td>10.100</td>
+
<td>0.090</td>
+
<td>0.450</td>
+
</tr>
+
<tr>
+
<td></td>
+
<td>0.075</td>
+
<td>0.125</td>
+
<td>11.100</td>
+
<td></td>
+
<td></td>
+
</tr>
+
<tr>
+
<td>-20%</td>
+
<td>0.750</td>
+
<td>0.100</td>
+
<td>12.700</td>
+
<td>0.138</td>
+
<td>0.690</td>
+
</tr>
+
+
</tbody>
+
</table>
+
</div>
+
</div>
+
<div class="card">
+
<div class="card-header"><h3>2. change beta</h3></div>
+
<img class="card-img-top" style="width: 512px;" src="../img/Model_5.png" />
+
<div class="card card-body">
+
<table >
+
<thead class="table table-active">
+
<tr>
+
<th>change β</th>
+
<th>α</th>
+
<th>β</th>
+
<th>LT50</th>
+
<th>ΔLT50</th>
+
<th>Ratio</th>
+
</tr>
+
</thead>
+
<tbody>
+
<tr>
+
<td>+20%</td>
+
<td>0.900</td>
+
<td>0.125</td>
+
<td>10.300</td>
+
<td>0.072</td>
+
<td>0.360</td>
+
</tr>
+
<tr>
+
<td></td>
+
<td>0.750</td>
+
<td>0.125</td>
+
<td>11.100</td>
+
<td></td>
+
<td></td>
+
</tr>
+
<tr>
+
<td>-20%</td>
+
<td>0.600</td>
+
<td>0.125</td>
+
<td>12.400</td>
+
<td>0.117</td>
+
<td>0.585</td>
+
</tr>
+
+
</tbody>
+
</table>
+
</div>
+
</div>
+
</div>
+
+
<div>
+
+
</div>
+
<div>
+
+
</div>
+
<p>The last term Ratio is the normalized sensitivities-the ratio of the relative change of the output to the relative change of the parameter.</p>
+
</div>
+
</div>
+
</div>
+
</div>
+
+
<div class="col-2">
+
<!--
+
<div class="position-fixed text-muted">
+
<a class="text-muted" href="#">Modle</a>
+
    <a class="text-muted" href="#Assumptions">Assumptions</a>
+
    <a class="" href="#Natural condition">Natural condition</a>
+
    <a class="" href="#With infection">With infection</a>
+
    <a class="" href="#Parameters">Parameters</a>
+
    <a class="" href="#Results">Results</a>
+
    <a class="" href="#Sensitivity Analysis">Sensitivity Analysis</a>
+
</div>
+
-->
+
+
+
<div>
+
<ul class=" position-fixed border-left">
+
<p><a class="text-muted" href="#">Modle</a></p>
+
<p><a class="text-muted" href="#Assumptions">Assumptions</a></p>
+
<p><a class="text-muted" href="#Natural condition">Natural condition</a></p>
+
+
<p><a class="text-muted" href="#With infection">With infection</a></p>
+
<p><a class="text-muted" href="#Parameters">Parameters</a></p>
+
<p><a class="text-muted" href="#Results">Results</a></p>
+
<p><a class="text-muted" href="#Sensitivity Analysis">Sensitivity Analysis</a></p>
+
</ul>
+
</div>
+
+
<!--
+
<nav id="navbar" class="sidebar flex-column position-fixed border-left ">
+
 
+
  <nav class="nav nav-pills flex-column text-muted">
+
  <a class="nav-link text-muted" href="#">Modle</a>
+
    <a class="nav-link text-muted" 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 class="col-2">
+
<div id="list-example" class="list-group list-group-flush position-fixed">
+
<a class="list-group-item list-group-item-action active" href="#">Model</a>
+
<a class="list-group-item list-group-item-action" href="#Assumptions">Assumptions</a>
+
            <a class="list-group-item list-group-item-action" href="#Natural condition">Natural condition</a>
+
            <a class="list-group-item list-group-item-action" href="#With infection">With infection</a>
+
          <a class="list-group-item list-group-item-action" href="#Parameters">Parameters</a>
+
          <a class="list-group-item list-group-item-action" href="#Results">Results</a>
+
            <a class="list-group-item list-group-item-action" href="#Sensitivity Analysis">Sensitivity Analysis</a>
+
         
+
</div>
+
</div>
+
-->
+
</div>
+
+
+
+
<div></div>
+
</div>
+
+
 
+
 
+
</body>
+
 
</html>
 
</html>

Revision as of 14:34, 28 September 2018

MODEL



Introduction



This year, our team creates a mathematical representation of our concrete self-healing system. This representation, or model, constructs a judging scale with which we can utilize to regulate the four main environmental factors affecting our final concrete healing efficiency (reflected by mineralization activity).



The four factors are:

Based on this model, we can also design the best ‘package'– the microcapsule shells with adequate nutrition combination.


This modelling process, presented below, can be seen as a feedback from the wet lab (experiment result) and the dry lab(modelling analysis).



The following page shows how we conducted modelling approaches to achieve our goals. To begin with, we make some definition.





Definition




    Variables and nomenclature:
  •         concentration of spores - c[Spore]
  •         concentration of carbon source(C3H5O3Na) - c[C3H5O3Na]
  •         concentration of nitrogen source(NaNO3) - c[NaNO3]
  •          pH of the media - pH

    Goal of model:
  •         Use preliminary data to guide future experiments.


Procedure


1.  Standardizing variables: transforming variables into the same scale. Here we utilize the z-score standardizing method, such that



where the symbols in the formula are variable, standardized variable, sample mean, andstandard deviation respectively.

2.  Fitting functions of each variable with polynomial function.

3.  Getting the overall relationship using linear least square method.




Results


The graph above depicts the polynomial regression of 4 factors. In general, for each of those four factors, the mineralization activity shows the similar tendency of going up first and down later with the increase of each factor.



The fitted functions are showed as follows:


We then fitted x with y linearly to get the overall regression equation which can describe the weight of four variables respectively.



Overall regression equation:



with




where ai1represents the coefficient of the highest power term of curve-fitting equation, such that the coefficient of the highest power term are of the same scale.




From the equation above, we can see that nitrogen source has the maximum weight, while pH has the minimum weight, which means nitrogen source is the most essential nutrition for B.subtilis spore. Additionally, the low weight of pH shows that spores are not sensitive to the change of pH in a relatively apt range, although there is a sharp decline in activity when pH reaches 11.



In summary, as long as we keep the environment below the boundary high-pH point, it makes no much difference how much we have improved the alkaline resistance of B.subtilis spore. In this way, this model instructs us on how to modify our formulation, such that the ingredients can be made full use of.