Difference between revisions of "Team:NCTU Formosa/Disease Occurrence Model"

 
(18 intermediate revisions by 4 users not shown)
Line 17: Line 17:
 
     margin: 0;
 
     margin: 0;
 
     padding: 0;
 
     padding: 0;
    height: 100vh;
 
 
     overflow-x: hidden;
 
     overflow-x: hidden;
 +
    height: 100vh;
 
}
 
}
  
Line 24: Line 24:
 
overflow-x: hidden !important;
 
overflow-x: hidden !important;
 
-ms-overflow-x: hidden !important;
 
-ms-overflow-x: hidden !important;
 +
}
 +
#fut{
 +
background: white;
 +
width: 100vw;
 +
height: 50px;
 +
}
 +
 +
.line{
 +
position: relative;
 +
left: 10vw;
 +
top: 10vmax;
 +
}
 +
 +
.latex{
 +
    font-size: 20px;
 +
    max-width: 60%;
 +
}
 +
 +
.latex2{
 +
    font-size: 15px;
 +
    max-width: 60%;
 
}
 
}
  
p{
+
p {
 
     text-align: justify;
 
     text-align: justify;
 
     font-size: 1.3em !important;
 
     font-size: 1.3em !important;
 
     font-family: Arial, sans-serif;
 
     font-family: Arial, sans-serif;
 +
    margin: 10px 0 0;
 
}
 
}
  
.subtitle{
+
h1{
     height: 35px;
+
    border: 0;
     max-width: 90%;
+
}
 +
 
 +
h4{
 +
    text-align: center;
 +
    margin-top: 0px;
 +
    margin-bottom: 50px !important;
 +
    font-size: 15px;
 +
}
 +
 
 +
h3{
 +
    text-align: center;
 +
    margin-top: 50px;
 +
    margin-bottom: 0px !important;
 +
    font-size: 15px;
 +
}
 +
 
 +
@font-face {
 +
  font-family: neo_latina;
 +
  src: url(https://static.igem.org/mediawiki/2017/4/46/Neo-latina-demo-FFP.ttf);
 +
}
 +
 
 +
.subtitle>h6{
 +
     font-family: neo_latina;
 +
     font-size: 3.9em;
 
     position: relative;
 
     position: relative;
     left: 10vw;
+
     left: 25px;
 
     margin-top: 50px;
 
     margin-top: 50px;
 
     margin-bottom: 10px;
 
     margin-bottom: 10px;
 +
    line-height: 1em;
 
}
 
}
 +
 +
.subtitle>h5{
 +
    font-family: neo_latina;
 +
    font-size: 1.95em;
 +
    position: relative;
 +
    left: 50px;
 +
    margin-top: -10px;
 +
    margin-bottom: 10px;
 +
    line-height: 1em;
 +
}
 +
 +
 +
 +
.dsocm_content{
 +
    margin: 10vh 10vw 0vh;
 +
}
 +
 +
.sub_text{
 +
    width:65vw;
 +
    position: relative;
 +
    left: 5vw;
 +
}
 +
 +
.sub2_text{
 +
    width:60vw;
 +
    position: relative;
 +
    left: 10vw;
 +
    line-height: 15px;
 +
    margin-bottom: 10px;
 +
}
 +
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------*/
  
Line 168: Line 245:
 
.hide>img{
 
.hide>img{
 
     margin-bottom: 50px;
 
     margin-bottom: 50px;
 +
}
 +
/*----------------------------------------------------------------------------*/
 +
span>a:link {
 +
    color: #85CFEA;
 
}
 
}
  
 +
span>a:visited {
 +
    color: #85CFEA;
 +
}
 +
 +
span>a:hover {
 +
    color: blue;
 +
}
 +
 +
span>a:active {
 +
    color: blue;
 +
}
 +
 +
span>a{
 +
    text-decoration: none;
 +
}
  
  
Line 176: Line 272:
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------*/
 
#first:hover{
 
    width: 23vw;
 
    left: 8.5vw;
 
}
 
 
#third:hover{
 
    width: 23vw;
 
    left: 68.5vw;
 
}
 
  
 
</style>
 
</style>
 
<!----------------------------------------------------------------------------->
 
<!----------------------------------------------------------------------------->
 
<script>
 
<script>
 +
$(document).ready(function(){
 +
    $(".show_pic").click(function(){
 +
        $(".hide").show(500);
 +
        $(".show").hide();
 +
    });
 +
});
 +
 +
$(document).ready(function(){
 +
    $(".hide_pic").click(function(){
 +
        $(".hide").hide();
 +
        $(".show").show(500);
 +
    });
 +
});
 +
 +
$(document).ready(function(){
 +
    $(".show2_pic").click(function(){
 +
        $(".hide2").show(500);
 +
        $(".show2").hide();
 +
    });
 +
});
 +
 +
$(document).ready(function(){
 +
    $(".hide2_pic").click(function(){
 +
        $(".hide2").hide();
 +
        $(".show2").show(500);
 +
    });
 +
});
  
 
</script>
 
</script>
 +
 +
    <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js" type="text/javascript">
 +
        MathJax.Hub.Config({
 +
            extensions: ["tex2jax.js", "TeX/AMSmath.js", "TeX/AMSsymbols.js"],
 +
            jax: ["input/TeX", "output/HTML-CSS"],
 +
            tex2jax: {
 +
                inlineMath: [
 +
                    ['$', '$'],
 +
                    ["\\(", "\\)"]
 +
                ],
 +
                displayMath: [
 +
                    ['$$', '$$'],
 +
                    ["\\[", "\\]"]
 +
                ],
 +
            },
 +
            "HTML-CSS": {
 +
                availableFonts: ["TeX"]
 +
            }
 +
        });
 +
    </script>
 
<!----------------------------------------------------------------------------->
 
<!----------------------------------------------------------------------------->
  
Line 215: Line 348:
  
 
         <!----------------------------------------------------------------------------->
 
         <!----------------------------------------------------------------------------->
         <img id="dsocm_pic" class="subtitle" src="https://static.igem.org/mediawiki/2017/0/0e/HpAg_study.png">
+
         <div class="subtitle">
 +
            <h6>Prediction Modeling</h6></div>
  
 
         <div class="dsocm_dsocm">
 
         <div class="dsocm_dsocm">
 
             <div class="show">
 
             <div class="show">
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In order to cure the fungal diseases in reality, we used the method called Convolutional Neural Network(CNN) to catch the weather patterns that was hard to be recognized by humans. With this disease occurrence model, we can
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In order to cure the fungal diseases in reality, we used the method called Convolutional Neural Network(CNN) to catch the weather patterns that was hard to be recognized by humans. With this disease occurrence model, we can get the daily possibility of disease occurrences. Also, we set up a warning and auto-spraying system by combining the model with the IoTtalk to apply peptides into farmlands.
                    get the daily possibility of disease occurrences. Also, we set up a warning and auto-spraying system by combining the model with the IoTtalk to apply peptides into farmlands.
+
 
                 </p>
 
                 </p>
  
Line 241: Line 374:
 
                 <p>
 
                 <p>
 
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;For the disease occurrence prediction modeling, NCTU_Formosa built a model predicting the diseases relating to everyday weather information that was based on the neural network. There were two kinds of data in our predict
 
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;For the disease occurrence prediction modeling, NCTU_Formosa built a model predicting the diseases relating to everyday weather information that was based on the neural network. There were two kinds of data in our predict
                     system, which were the recorded diseases happening data that were collected from the government agency and the weather data from the <a href="http://e-service.cwb.gov.tw/HistoryDataQuery/index.jsp" target="_blank">Central Meteorological Bureau's website</a>Central
+
                     system, which were the recorded diseases happening data that were collected from the government agency and the weather data from the <span><a href="http://e-service.cwb.gov.tw/HistoryDataQuery/index.jsp" target="_blank">Central Meteorological Bureau's website</a></span> that corresponded to the data forward. Then, we combined these two and then deleted the recorded data that had no weather data to match with.
                    Meteorological Bureau's website that corresponded to the data forward. Then, we combined these two, and then deleted the recorded data that had no weather data to match with.
+
                </p>
 +
                <p>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Finally, our data contain weather feature and the Label. The weather feature is a two-dimensional array which has 14 days x 11 features. The 11 features include relative humidity, rainfall and the maximum, minimum, average of the temperature and air pressure.
 +
The Label contains two classes, negative(diseases that would not happen) and positive(diseases that would happen).
 
                 </p>
 
                 </p>
  
 
                 <img src="https://static.igem.org/mediawiki/2017/e/e3/Dsocm_flow2.png" width="60%" style="display:block; margin:auto;">
 
                 <img src="https://static.igem.org/mediawiki/2017/e/e3/Dsocm_flow2.png" width="60%" style="display:block; margin:auto;">
 +
                <h4>Figure 1: Data process flow chart</h4>
  
 
                 <p>
 
                 <p>
Line 253: Line 390:
  
 
                 <h1>Convolutional Neural Network(CNN)</h1>
 
                 <h1>Convolutional Neural Network(CNN)</h1>
 
+
                <img src="https://static.igem.org/mediawiki/2017/7/71/CNN_flowchart_2.png" width="60%" style="display:block; margin:auto;">
 +
                <h4>Figure 2: This is an overview of our model,it contains the convolution layer, max pooling layer and multi full connection layer.</h4>
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Convolutional Neural Network(CNN) is a powerful network that can automatically recognize the patterns of the features in a period of time, just like the favorite weather change to the fungal diseases in the life cycle. In
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Convolutional Neural Network(CNN) is a powerful network that can automatically recognize the patterns of the features in a period of time, just like the favorite weather change to the fungal diseases in the life cycle. In other words, it means it can recognize the pathogen life cycle in a specific time and mark them out because life cycle’s patterns can happen everywhere over a period of time. For example, our model used weather data for the past two weeks as the model input, and the weather patterns can be found in this 14-day data.
                    other words, it means it can recognize the pathogen life cycle in a specific time and mark them out because life cycle’s patterns can happen everywhere over a period of time. For example, our model used weather data for tje past two
+
                    weeks as the model input, and the weather patterns can be found in this 14-day data.
+
 
                 </p>
 
                 </p>
  
Line 265: Line 401:
 
                 </p>
 
                 </p>
  
                 <p>一個公式</p>
+
                 <div class="latex">$$ f(x_i) = max(x_i) $$</div>
  
 
                 <p>
 
                 <p>
Line 271: Line 407:
 
                     =5 and so on.
 
                     =5 and so on.
 
                 </p>
 
                 </p>
 
+
                <p>
                 <p>一張圖</p>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Because the max pooling output is the two-dimensional tensor, we will flatten the max pooling output to one dimensional tensor for further full connection layer.
 +
                 </p>
 +
<p>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;After flattening the max pooling layer, we will use the full connection layer to classified the max pooling result, full connection layer is a basic Neural Network layer that that can switch the max pooling layer output into the high dimensional space and then classified them into two classes, negative(diseases that would not happen) and positive(diseases that would happen).
 +
                </p>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Finally, after using the full connection layer to convert the max pooling output into high dim space, the model then classified them into two classes, negative(disease will not happen) and positive(disease will happen).
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Finally, after using the full connection layer to convert the max pooling output into high dim space, the model then classified them into two classes, negative(diseases that didn’t happen) and positive(diseases that haven’t
 +
                    happened).
 
                 </p>
 
                 </p>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;However, the network output was (a number that) hard to be realized by humans , so we used the softmax function to transform it into the diseases happening probability. Here is the formula of the softmax function:
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;However, the network output was (a number that) hard to be realized by humans , so we used the <b>softmax function</b> to transform it into the diseases happening probability. Here is the formula of the
 +
                    softmax function:
 
                 </p>
 
                 </p>
  
                 <p>一個公式</p>
+
                 <div class="latex">$$ {\displaystyle \sigma (\mathbf {z} )_{j}={\frac {e^{z_{j}}}{\sum _{k=1}^{K}e^{z_{k}}}}} $$</div>
  
                 <p>
+
                 <p style="margin-bottom: 50px;">
 
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;So on, the model output can be easily realized by machines and humans, and next let’s talk about how we defined our model cost function and the optimizer.
 
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;So on, the model output can be easily realized by machines and humans, and next let’s talk about how we defined our model cost function and the optimizer.
 
                 </p>
 
                 </p>
  
 
                 <img src="https://static.igem.org/mediawiki/2017/0/05/Dsocm_flow3.png" width="60%" style="display:block; margin:auto;">
 
                 <img src="https://static.igem.org/mediawiki/2017/0/05/Dsocm_flow3.png" width="60%" style="display:block; margin:auto;">
 +
                <h4>Figure 3: Accuracy flow chart</h4>
  
 
                 <p>
 
                 <p>
Line 296: Line 439:
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Suppose His a cross-entropy function, y'i is the real label, and the yiis the network prediction output, the cost is
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Suppose H is a cross-entropy function, y'<sub>i</sub> is the real label, and the y<sub>i</sub> the network prediction output, the cost is
 
                 </p>
 
                 </p>
  
                 <p>一個公式</p>
+
                 <div class="latex">$$ H_{y'}(y) = - \sum_{i} y_{i}^{'}log(y_{i}) $$</div>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Parameters of the Neural Network were being optimized by Adam optimizer, which is a moat commonly used ways to optimize the network.
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Parameters of the Neural Network were being optimized by Adam optimizer, which is a most commonly used ways to optimize the network.
 
                 </p>
 
                 </p>
  
                 <p>
+
                 <p style="margin-bottom: 20px;">
 
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;After training, the model will be tested by independence test data, and the result was shown below.
 
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;After training, the model will be tested by independence test data, and the result was shown below.
 
                 </p>
 
                 </p>
  
                 <p>一張圖片</p>
+
                 <img src="https://static.igem.org/mediawiki/2017/c/ce/Dsocm_photo17.PNG" width="80%" style="display:block; margin:auto;">
 +
                <h4>Figure 4: Accuracy score = 82.5%</h4>
  
 
                 <p>
 
                 <p>
Line 320: Line 464:
  
 
         <!----------------------------------------------------------------------------->
 
         <!----------------------------------------------------------------------------->
         <img id="spgnm_pic" class="subtitle" src="https://static.igem.org/mediawiki/2017/0/0e/HpAg_study.png">
+
         <div class="subtitle">
 +
            <h6>Spore Germination Modeling</h6></div>
  
 
         <div class="dsocm_spgnm">
 
         <div class="dsocm_spgnm">
 
             <div class="show2">
 
             <div class="show2">
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;To predict the occurrence of diseases more effectively and timely, we made the spore germination rate. According to the researches, we first fitted out a linear equation for the spore germination rate base on the humidity
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;To predict the occurrence of diseases more effectively and timely, we made the spore germination rate. According to the researchers, we first fitted out a linear equation for the spore germination rate base on the humidity and a cubic equation for the spore germination rate base on the temperature. By multiplying them, we get the general model for the spore germination. At last, the model was combined to IoTtalk to apply peptides into farmlands.
                    and a cubic equation for the spore germination rate base on the temperature. By multiply them, we get the general model for the spore germination. At last, the model was combined to IoTtalk to apply peptides into farmlands.
+
  
 
                 </p>
 
                 </p>
Line 341: Line 485:
  
 
             <div class="hide2">
 
             <div class="hide2">
 +
 +
                <h1>Spore Germination Rate</h1>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Because we found out that the humidity affected the spore germination the most and the temperature secondly. Therefore, we wanted to found the relationships between temperature with the spore germination rate and the relationships
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Because we found out that the humidity and temperature affected the spore germination the most. Therefore, we wanted to build a general model for the spore germination rate based on temperature or humidity with different fungal species that can fit in every fungal species.
                    between humidity with the spore germination rate in different fungal species and tried to make a general model that can fit in every fungal species.
+
 
                 </p>
 
                 </p>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;At first, we used the spore germination data on the research papers to fit out the functions. It found out that the spore germination rate based on temperature can be fitted in a cubic equation.
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;At first, we used the spore germination data on the research papers to fit out the functions.
 
                 </p>
 
                 </p>
  
                 <p>一個方程式</p>
+
                 <div>
 +
                    <h2>Spore germination based on temperature <sup>[1]</sup>:</h2>
  
                <p>
+
                    <p>
                    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;And the spore germination rate based on the humidity can be fitted in a linear equation.
+
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(1) Myce-liophthora thermophila's spore germination rate based on temperature
                </p>
+
                    </p>
  
                <p>一個方程式</p>
+
                    <div class="latex">$$ y = 0.0004x^3 - 0.0132x^2 + 4.0447x -24.746 $$</div>
  
                <p>
+
                    <img src="https://static.igem.org/mediawiki/2017/e/e0/Dsocm_photo1.png" width="60%" style="display: block; margin: auto;">
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Later, we considered the two equations as independent events and multiplied them to form our general fungal spore germination model. Then we conducted experiments to determine the coefficients and verification.
+
                     <h4>Figure 5: Spore germination based on temperature.</h4>
                </p>
+
  
                <p>一個方程式</p>
+
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(2) Aspergillus niger’s spore germination rate based on temperature:
 +
                    </p>
  
                <h1>Method:</h1>
+
                    <div class="latex">$$ y = 0.0326x^3 - 4.0593x^2 + 160.2x - 1943.7 $$</div>
  
                 <h3>1. Fix humidity, change temperature:</h3>
+
                    <img src="https://static.igem.org/mediawiki/2017/2/21/Dsocm_photo2.png" width="60%" style="display: block; margin: auto;">
 +
                    <h4>Figure 6: Aspergillus niger’s spore germination rate based on temperature.</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(3) P. oryzae’s spore germination rate based on temperature<sup>[2]</sup>:
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ y = 0.06x^3 - 4.389x^2 + 106.86x - 774.66 $$</div>
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/f/f5/Dsocm_photo3.png" width="60%" style="display: block; margin: auto;">
 +
                    <h4>Figure 7: P. oryzae’s spore germination rate based on temperature.</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(4) Diplodia corticola’s spore germination rate based on temperature<sup>[3]</sup>:
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ y = -0.0041x^3 + 0.0169x^2 + 7.1531x - 32.071 $$</div>
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/8/8c/Dsocm_photo4.png" width="60%" style="display: block; margin: auto;">
 +
                    <h4>Figure 8: Diplodia corticola’s spore germination rate based on temperature.</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;As results, we found out that the spore germination rate based on temperature can be fitted in a cubic equation
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ f_1(x) = ax^3 + bx^2 + cx + d $$</div>
 +
                 </div>
 +
 
 +
                <div style="margin-top: 50px;">
 +
                    <h2>Spore germination based on relative humidity:</h2>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(1) Aspergillus niger’s spore germination rate based on relative humidity<sup>[4]</sup>:
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ y = 352.38x - 254.14 $$</div>
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/7/7d/Dsocm_photo5.png" width="60%" style="display: block; margin: auto;">
 +
                    <h4>Figure 9: Aspergillus niger’s spore germination rate based on relative humidity</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(2) Pseudocercospora’s spore germination rate based on relative humidity<sup>[5]</sup>:
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ y = 0.1071x - 10.165 $$</div>
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/9/9e/Dsocm_photo6.png" width="60%" style="display: block; margin: auto;">
 +
                    <h4>Figure 10: Pseudocercospora’s spore germination rate based on relative humidity.</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;As results, we found out the spore germination rate based on the relative humidity in a linear equation.
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ f_2(x) = ax+b $$</div>
 +
 
 +
                    <h2 style="margin-top: 50px;">The General model of spore germination rate:</h2>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Later, we considered the two equations as independent events and multiplied them to form our general fungal spore germination model. Then we conducted experiments to determine the coefficients and verification.
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ f_1 \times f_2 $$</div>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In other words, Just need the temperature and humidity conditions, we can calculate the spore germination in that environment.
 +
                    </p>
 +
                </div>
 +
 
 +
                <div style="margin-top: 50px;">
 +
                    <h2>Experiment Method:</h2>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The experiment part was divided into two part:
 +
                    </p>
 +
                    <h3>Table 1: Overall experimental disign</h3>
 +
                    <img src="https://static.igem.org/mediawiki/2017/d/d9/Spgm_photo.PNG" width="60%" style="display: block; margin: auto;">
 +
 
 +
 
 +
                    <div class="sub_text">
 +
                        <p style="margin-top: 20px;">
 +
                            (1) Fix humidity, change temperature:
 +
                        </p>
 +
 
 +
 
 +
                        <p>
 +
                            &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Remove the spores from fungal plate by Distillation-Distillation H<sub>2</sub>O to made the spore suspension solution (2x10<sup>5</sup>particles/mL) and 2% glucose solution mixing with equal volume,
 +
                            then put the solution in concave glass slides. Later, put it into temperature and humidity control box. The humidity was set at 100%, the temperature was set from 10 to 30 Celsius degrees and been tested once every 5 degrees.
 +
                        </p>
 +
 
 +
                        <p>
 +
                            (2) Fix temperature, change humidity:
 +
                        </p>
 +
 
 +
                        <p>
 +
                            &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Remove the spores from fungal plate by Distillation-Distillation H2O to made the spore suspension solution(2x10<sup>5</sup>particles/mL)) and 2% glucose solution mixing with equal volume, then
 +
                            put the solution in concave glass slides. Later, put it into temperature and humidity control box. The temperature was fixed at 25 Celsius degrees, the humidity was set from 80% to 100% and been tested once every 5 percents.
 +
 
 +
                        </p>
 +
 
 +
                        <p>
 +
                            (3)Independent event validation:
 +
                        </p>
 +
 
 +
                        <p>
 +
                            &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;We compared the result of our formula with the value of the spore germination in the reality to verify temperature and relative humidity are independent event to spore germination. We randomly choose the conditions of temperature
 +
                            23 and 13 Celsius, relative humidity 97% and 80%.
 +
                        </p>
 +
 
 +
 
 +
                    </div>
 +
                </div>
 +
 
 +
                <h1>Experiment Result:</h1>
 +
                <p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The standard of spore germination</p>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Washed the spores from fungal plate by water(DDW), then put the spore (suspension) water in concave glass slides. Later, put it into temperature and humidity control box. The humidity was set at 100%, the temperature was
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;following pictures are the Illustration of whether the spore germinate or not.
                    set from 10 to 30 Celsius degrees and been tested once every 5 degrees.
+
 
                 </p>
 
                 </p>
  
                 <h3>2. Fix temperature, change humidity</h3>
+
                 <div style="position: relative; left: 10%; margin-top: 50px;">
 +
                    <img src="https://static.igem.org/mediawiki/2017/f/f2/Dsocm_photo7.jpeg" width="40%">
 +
                    <img src="https://static.igem.org/mediawiki/2017/7/75/Dsocm_photo8.png" width="40%">
 +
                </div>
 +
                <h4 style="margin-bottom: 0 !important;">Figure 11: (Left)10 Celsius degrees and 100% relative humidity (9 hours): Spore had not germinated</h4>
 +
                <h4 style=" margin-top: 0;">Figure 12: (Right)25 Celsius degrees and 100% relative humidity (9 hours): Spore had germinated</h4>
 +
 
 +
 
 +
 
 +
                <div class="sub_text">
 +
                    <p>
 +
                        (1) Fixed relative humidity in 100%:
 +
                    </p>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The effective range of temperature is 10 to 30 Celsius degrees.
 +
                    </p>
 +
 
 +
                    <h3>Table 2: The spore germination rate of fixing the relative humidity and changing the temperature in different hours.</h3>
 +
                    <img src="https://static.igem.org/mediawiki/2017/b/bb/Dsocm_photo11.PNG" width="80%" style="display: block; margin: auto;">
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/2/2d/Dsocm_photo12.png" width="60%" style="display: block; margin: 20px auto 0 auto;">
 +
                    <h4>Figure 13: <i>Botrytis cinerea</i>’s spore germination rate based on temperature.(9 hours)</h4>
 +
 
 +
                    <p style="margin-top: 50px;">
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<i>Botrytis cinerea</i>’s spore germination rate based on temperature:
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ f_1 = ( -0.0625x_1^3 + 2.9974x_1^2 - 37.865x_1 + 141.68) $$</div>
 +
 
 +
                    <p>
 +
                        (2) Fixed temperature in 20 Celsius
 +
                    </p>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The effective range of relative humidity is 70%~100%, we use linear translation to let the spore germination in 100% relative humidity is also 100%.
 +
                    </p>
 +
 
 +
                    <h3>Table 3: </h3>
 +
                    <img src="https://static.igem.org/mediawiki/2017/5/59/Dsocm_photo9.png" width="80%" style="display: block; margin: auto;">
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/9/97/Dsocm_photo10.png" width="60%" style="display: block; margin: 20px auto 0 auto;">
 +
                    <h4>Figure 14: <i>Botrytis cinerea</i>’s spore germination rate based on relative humidity(9 hours)</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<i>Botrytis cinerea</i>’s spore germination rate based on relative humidity:
 +
                    </p>
 +
 
 +
                    <div class="latex">$$ f_2 = 316.88x - 216.88 $$</div>
 +
 
 +
                    <p>(3) Independent event validation</p>
 +
 
 +
                    <h3>Table 4: The results of our independent event validation.</h3>
 +
                    <img src="https://static.igem.org/mediawiki/2017/a/a4/Dsocm_photo14.png" width="60%" style="display: block; margin: auto;">
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/4/4e/Dsocm_photo15.png" width="60%" style="display: block; margin: 20px auto 0 auto;">
 +
                    <h4>Figure 15: The spore germination under the condition of 23 Celsius degrees and 97% relative humidity (9 hours) </h4>
 +
 
 +
                    <p style="margin-bottom: 20px;">
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In the conditions of 23 Celsius degrees and 97% relative humidity (9 hours) the spore germination rate by experiment is 92.45%. And according to our formula that the spore germination is 86.84%.
 +
                    </p>
 +
 
 +
                    <img src="https://static.igem.org/mediawiki/2017/9/9b/Dsocm_photo16.png" width="60%" style="display: block; margin: auto;">
 +
                    <h4>Figure 16: The spore germination under the condition of 13 Celsius degrees and 80% relative humidity (9 hours).</h4>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In the condition of 13 Celsius degrees and 80% relative humidity (9 hours) the spore germination rate by experiment is 5.41%. And according to our formula that the spore germination is 6.84%.
 +
                    </p>
 +
 
 +
                    <p>
 +
                        &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The result can proof that our formula is quite precise.
 +
                    </p>
 +
 
 +
                    <h3>Table 5: The results of comparing the spore germination rate of our experiment and formula.</h3>
 +
                    <img src="https://static.igem.org/mediawiki/2017/8/83/Compares.PNG" width="60%" style="display: block; margin: auto;">
 +
 
 +
 
 +
                </div>
 +
 
  
 +
                <h2 style="margin-top: 50px;">Conclusion: </h2>
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Washed the spores from fungal plate by water(DDW), then put the spore (suspension) water in concave glass slides. Later, put it into temperature and humidity control box. Temperature was fixed at 25 Celsius degrees, the humidity
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Our final formula is:
                    was set from 80% to 100% and been tested once every 5 percents.
+
 
                 </p>
 
                 </p>
 
+
                 <div class="latex2">$$ f_1 \times f_2 = [(-0.0625x_1^3 + 2.9974x_1^2 - 37.865x_1 + 141.8)\div 100]\times [(316.88x_2 - 216.88)\div 100] $$</div>
                 <p>一個圖表</p>
+
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;By modeling combined with the IoT system, two aspects of disease prediction are provided. Respectively:
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;where x<sub>1</sub> is temperature and x<sub>2</sub> is relative humidity.
 
                 </p>
 
                 </p>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;IoT immediately detect temperature and humidity (per hour detection), and we put the temperature and humidity condition that IoT detected into the formula we developed to calculate spore germination. If the spore germination
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;According our experiment result, we prove that temperature and relative humidity is an independent event and our formula is quite precise.
                    rate exceeds certain value, we will inform the user that the disease may occur.(How to develop see below)
+
 
 
                 </p>
 
                 </p>
 +
 +
 +
                <h1>Combined with IoTtalk</h1>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;If both of the predicted results are above the threshold, the user will be advised to spray the pesticide. The data send back to IoT system simultaneously and the user can decide whether to automatically spray.
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;IoT immediately detects temperature and humidity (per hour detection), and we put the temperature and humidity condition that IoT detected into the formula we developed to calculate spore germination. If the spore germination rate exceeds a certain value, we will inform the user that the disease may occur.
 
                 </p>
 
                 </p>
  
 
                 <p>
 
                 <p>
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In other words, Just need the temperature and humidity conditions, we can calculate the spore germination in that environment.
+
                     &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;If both of the predicted results are above the threshold, the user will be advised to spray the pesticide. The data sent back to IoT system simultaneously and the user can decide whether to automatically spray.
 
                 </p>
 
                 </p>
  
Line 404: Line 743:
 
                 </p>
 
                 </p>
  
                 <h1>Result</h1>
+
                 <h1>Reference</h1>
 
+
                <h3>Formula fit result:</h3>
+
 
+
                <h3>Data found from paper:</h3>
+
 
+
                <p>一張圖片</p>
+
 
+
 
+
  
                 <img src="#" width="60%" style="display:block; margin:auto;">
+
                 <p>[1]<small>H. HASSOUNI1, et al. COMPARATIVE SPORE GERMINATION OF FILAMENTOUS FUNGI ON SOLID STATE FERMENTATION UNDER DIFFERENT CULTURE CONDITIONS. 2006</small></p>
 +
                <p>[2]<small>Reference:Shivani R, Top 6 Factors Influencing Disease Cycle of Rice | Plant Disease</small></p>
 +
                <p>[3]<small>Reference:José R. Úrbez-Torres, Effect of Temperature on Conidial Germination of Botryosphaeriaceae Species Infecting Grapevines, 2010</small></p>
 +
                <p>[4]<small>Awad M. Abdel-Rahim, Hassan A. Arbab, Factors affecting spore germination in Aspergillus niger, 1985</small></p>
 +
                <p>[5]<small>L. H. Jacome, W. Schuh, and R. E. Stevenson Effect of Temperature and Relative Humidity on Germination and Germ Tube Development of Mycosphaerella fijiensis var. difformis, 1991</small></p>
  
 
                 <div><img class="hide2_pic" src="https://static.igem.org/mediawiki/2017/c/cb/Ptp_show.png" style="display:block; margin:auto;"></div>
 
                 <div><img class="hide2_pic" src="https://static.igem.org/mediawiki/2017/c/cb/Ptp_show.png" style="display:block; margin:auto;"></div>

Latest revision as of 02:48, 2 November 2017

navigation

NCTU_Formosa: Disease Occurrence Model
Prediction Modeling

     In order to cure the fungal diseases in reality, we used the method called Convolutional Neural Network(CNN) to catch the weather patterns that was hard to be recognized by humans. With this disease occurrence model, we can get the daily possibility of disease occurrences. Also, we set up a warning and auto-spraying system by combining the model with the IoTtalk to apply peptides into farmlands.

  1. Datasets--the past fungal diseases data and the weather data
  2. Convolutional Neural Network
  3. Softmax function
  4. How to define our model cost function and the optimizer
  5. Connection to the IoTtalk system

     For the disease occurrence prediction modeling, NCTU_Formosa built a model predicting the diseases relating to everyday weather information that was based on the neural network. There were two kinds of data in our predict system, which were the recorded diseases happening data that were collected from the government agency and the weather data from the Central Meteorological Bureau's website that corresponded to the data forward. Then, we combined these two and then deleted the recorded data that had no weather data to match with.

     Finally, our data contain weather feature and the Label. The weather feature is a two-dimensional array which has 14 days x 11 features. The 11 features include relative humidity, rainfall and the maximum, minimum, average of the temperature and air pressure. The Label contains two classes, negative(diseases that would not happen) and positive(diseases that would happen).

Figure 1: Data process flow chart

     In reality, the spore germination needs some specific weather conditions. Also, the plants may be influenced by the change of weather and become weak. However, this relationship is hard to be recognized by humans. Therefore, we used Convolutional Neural Network(CNN) as our method to catch the weather patterns.

Convolutional Neural Network(CNN)

Figure 2: This is an overview of our model,it contains the convolution layer, max pooling layer and multi full connection layer.

     Convolutional Neural Network(CNN) is a powerful network that can automatically recognize the patterns of the features in a period of time, just like the favorite weather change to the fungal diseases in the life cycle. In other words, it means it can recognize the pathogen life cycle in a specific time and mark them out because life cycle’s patterns can happen everywhere over a period of time. For example, our model used weather data for the past two weeks as the model input, and the weather patterns can be found in this 14-day data.

     After the CNN layer, the weather features then converted to weather change features, then a max pooling layer was added to filter some noises after the CNN layer. Weather patterns that caused diseases wouldn’t change in a short time, so the function of max pooling was to only return the maximum values in the filter.

$$ f(x_i) = max(x_i) $$

     For example: if the input array is [2,5,1,7,0,4] and the max pooling filter size is 2. When the filter step(distance the filter moves) is 1, the first max pooling output will be max(2,5)=5, and the second output is max(5,1) =5 and so on.

     Because the max pooling output is the two-dimensional tensor, we will flatten the max pooling output to one dimensional tensor for further full connection layer.

     After flattening the max pooling layer, we will use the full connection layer to classified the max pooling result, full connection layer is a basic Neural Network layer that that can switch the max pooling layer output into the high dimensional space and then classified them into two classes, negative(diseases that would not happen) and positive(diseases that would happen).

     Finally, after using the full connection layer to convert the max pooling output into high dim space, the model then classified them into two classes, negative(diseases that didn’t happen) and positive(diseases that haven’t happened).

     However, the network output was (a number that) hard to be realized by humans , so we used the softmax function to transform it into the diseases happening probability. Here is the formula of the softmax function:

$$ {\displaystyle \sigma (\mathbf {z} )_{j}={\frac {e^{z_{j}}}{\sum _{k=1}^{K}e^{z_{k}}}}} $$

     So on, the model output can be easily realized by machines and humans, and next let’s talk about how we defined our model cost function and the optimizer.

Figure 3: Accuracy flow chart

     We chose the cross-entropy as the network cost function because it performed well that the exclusion classification mission. Here this the formula of cross-entropy.

     Suppose H is a cross-entropy function, y'i is the real label, and the yi the network prediction output, the cost is

$$ H_{y'}(y) = - \sum_{i} y_{i}^{'}log(y_{i}) $$

     Parameters of the Neural Network were being optimized by Adam optimizer, which is a most commonly used ways to optimize the network.

     After training, the model will be tested by independence test data, and the result was shown below.

Figure 4: Accuracy score = 82.5%

     For the use in reality, the network would connect to our IoTtalk system, giving the user daily warning from the diseases happening.

Spore Germination Modeling

     To predict the occurrence of diseases more effectively and timely, we made the spore germination rate. According to the researchers, we first fitted out a linear equation for the spore germination rate base on the humidity and a cubic equation for the spore germination rate base on the temperature. By multiplying them, we get the general model for the spore germination. At last, the model was combined to IoTtalk to apply peptides into farmlands.

  1. Fitted out the equation with researches and get the general spore germination rate model.
  2. Conducted experiments as verification and determine the coefficients.
  3. Result

Spore Germination Rate

     Because we found out that the humidity and temperature affected the spore germination the most. Therefore, we wanted to build a general model for the spore germination rate based on temperature or humidity with different fungal species that can fit in every fungal species.

     At first, we used the spore germination data on the research papers to fit out the functions.

Spore germination based on temperature [1]:

     (1) Myce-liophthora thermophila's spore germination rate based on temperature

$$ y = 0.0004x^3 - 0.0132x^2 + 4.0447x -24.746 $$

Figure 5: Spore germination based on temperature.

     (2) Aspergillus niger’s spore germination rate based on temperature:

$$ y = 0.0326x^3 - 4.0593x^2 + 160.2x - 1943.7 $$

Figure 6: Aspergillus niger’s spore germination rate based on temperature.

     (3) P. oryzae’s spore germination rate based on temperature[2]:

$$ y = 0.06x^3 - 4.389x^2 + 106.86x - 774.66 $$

Figure 7: P. oryzae’s spore germination rate based on temperature.

     (4) Diplodia corticola’s spore germination rate based on temperature[3]:

$$ y = -0.0041x^3 + 0.0169x^2 + 7.1531x - 32.071 $$

Figure 8: Diplodia corticola’s spore germination rate based on temperature.

     As results, we found out that the spore germination rate based on temperature can be fitted in a cubic equation

$$ f_1(x) = ax^3 + bx^2 + cx + d $$

Spore germination based on relative humidity:

     (1) Aspergillus niger’s spore germination rate based on relative humidity[4]:

$$ y = 352.38x - 254.14 $$

Figure 9: Aspergillus niger’s spore germination rate based on relative humidity

     (2) Pseudocercospora’s spore germination rate based on relative humidity[5]:

$$ y = 0.1071x - 10.165 $$

Figure 10: Pseudocercospora’s spore germination rate based on relative humidity.

     As results, we found out the spore germination rate based on the relative humidity in a linear equation.

$$ f_2(x) = ax+b $$

The General model of spore germination rate:

     Later, we considered the two equations as independent events and multiplied them to form our general fungal spore germination model. Then we conducted experiments to determine the coefficients and verification.

$$ f_1 \times f_2 $$

     In other words, Just need the temperature and humidity conditions, we can calculate the spore germination in that environment.

Experiment Method:

     The experiment part was divided into two part:

Table 1: Overall experimental disign

(1) Fix humidity, change temperature:

     Remove the spores from fungal plate by Distillation-Distillation H2O to made the spore suspension solution (2x105particles/mL) and 2% glucose solution mixing with equal volume, then put the solution in concave glass slides. Later, put it into temperature and humidity control box. The humidity was set at 100%, the temperature was set from 10 to 30 Celsius degrees and been tested once every 5 degrees.

(2) Fix temperature, change humidity:

     Remove the spores from fungal plate by Distillation-Distillation H2O to made the spore suspension solution(2x105particles/mL)) and 2% glucose solution mixing with equal volume, then put the solution in concave glass slides. Later, put it into temperature and humidity control box. The temperature was fixed at 25 Celsius degrees, the humidity was set from 80% to 100% and been tested once every 5 percents.

(3)Independent event validation:

     We compared the result of our formula with the value of the spore germination in the reality to verify temperature and relative humidity are independent event to spore germination. We randomly choose the conditions of temperature 23 and 13 Celsius, relative humidity 97% and 80%.

Experiment Result:

     The standard of spore germination

     following pictures are the Illustration of whether the spore germinate or not.

Figure 11: (Left)10 Celsius degrees and 100% relative humidity (9 hours): Spore had not germinated

Figure 12: (Right)25 Celsius degrees and 100% relative humidity (9 hours): Spore had germinated

(1) Fixed relative humidity in 100%:

     The effective range of temperature is 10 to 30 Celsius degrees.

Table 2: The spore germination rate of fixing the relative humidity and changing the temperature in different hours.

Figure 13: Botrytis cinerea’s spore germination rate based on temperature.(9 hours)

     Botrytis cinerea’s spore germination rate based on temperature:

$$ f_1 = ( -0.0625x_1^3 + 2.9974x_1^2 - 37.865x_1 + 141.68) $$

(2) Fixed temperature in 20 Celsius

     The effective range of relative humidity is 70%~100%, we use linear translation to let the spore germination in 100% relative humidity is also 100%.

Table 3:

Figure 14: Botrytis cinerea’s spore germination rate based on relative humidity(9 hours)

     Botrytis cinerea’s spore germination rate based on relative humidity:

$$ f_2 = 316.88x - 216.88 $$

(3) Independent event validation

Table 4: The results of our independent event validation.

Figure 15: The spore germination under the condition of 23 Celsius degrees and 97% relative humidity (9 hours)

     In the conditions of 23 Celsius degrees and 97% relative humidity (9 hours) the spore germination rate by experiment is 92.45%. And according to our formula that the spore germination is 86.84%.

Figure 16: The spore germination under the condition of 13 Celsius degrees and 80% relative humidity (9 hours).

     In the condition of 13 Celsius degrees and 80% relative humidity (9 hours) the spore germination rate by experiment is 5.41%. And according to our formula that the spore germination is 6.84%.

     The result can proof that our formula is quite precise.

Table 5: The results of comparing the spore germination rate of our experiment and formula.

Conclusion:

     Our final formula is:

$$ f_1 \times f_2 = [(-0.0625x_1^3 + 2.9974x_1^2 - 37.865x_1 + 141.8)\div 100]\times [(316.88x_2 - 216.88)\div 100] $$

     where x1 is temperature and x2 is relative humidity.

     According our experiment result, we prove that temperature and relative humidity is an independent event and our formula is quite precise.

Combined with IoTtalk

     IoT immediately detects temperature and humidity (per hour detection), and we put the temperature and humidity condition that IoT detected into the formula we developed to calculate spore germination. If the spore germination rate exceeds a certain value, we will inform the user that the disease may occur.

     If both of the predicted results are above the threshold, the user will be advised to spray the pesticide. The data sent back to IoT system simultaneously and the user can decide whether to automatically spray.

     When the germination rate exceeds the threshold, we will advise the user to spray the corresponding peptide.

Reference

[1]H. HASSOUNI1, et al. COMPARATIVE SPORE GERMINATION OF FILAMENTOUS FUNGI ON SOLID STATE FERMENTATION UNDER DIFFERENT CULTURE CONDITIONS. 2006

[2]Reference:Shivani R, Top 6 Factors Influencing Disease Cycle of Rice | Plant Disease

[3]Reference:José R. Úrbez-Torres, Effect of Temperature on Conidial Germination of Botryosphaeriaceae Species Infecting Grapevines, 2010

[4]Awad M. Abdel-Rahim, Hassan A. Arbab, Factors affecting spore germination in Aspergillus niger, 1985

[5]L. H. Jacome, W. Schuh, and R. E. Stevenson Effect of Temperature and Relative Humidity on Germination and Germ Tube Development of Mycosphaerella fijiensis var. difformis, 1991

Untitled Document