| (12 intermediate revisions by 2 users not shown) | |||
| Line 131: | Line 131: | ||
<div class="col-xs-12 text-center menu-1 menu-wrap"> | <div class="col-xs-12 text-center menu-1 menu-wrap"> | ||
<ul> | <ul> | ||
| − | <li | + | <li><a href="index.html">Home</a></li> |
<li class="has-dropdown"> | <li class="has-dropdown"> | ||
<a href="Description">Project</a> | <a href="Description">Project</a> | ||
<ul class="dropdown" style="margin-top:-5px"> | <ul class="dropdown" style="margin-top:-5px"> | ||
<li ><a href="Description">Project</a></li> | <li ><a href="Description">Project</a></li> | ||
| − | <li><a href=" | + | <li><a href="Demostrate">Wet Lab</a></li> |
<li><a href="Model">Modeling</a></li> | <li><a href="Model">Modeling</a></li> | ||
<li><a href="Hardware">Device and Software</a></li> | <li><a href="Hardware">Device and Software</a></li> | ||
| Line 142: | Line 142: | ||
</li> | </li> | ||
<li><a href="Human_Practices">Human Practices</a></li> | <li><a href="Human_Practices">Human Practices</a></li> | ||
| − | + | <li><a href="Team">Our team</a></li> | |
| − | + | <li><a href="Collaborations">Collaborations</a></li> | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
</ul> | </ul> | ||
</div> | </div> | ||
| Line 157: | Line 151: | ||
<!-- </div> --> | <!-- </div> --> | ||
</nav> | </nav> | ||
| + | |||
| + | <header id="fh5co-header" class="fh5co-cover js-fullheight" role="banner" style="background-image: url(https://2019.igem.org/wiki/images/a/a4/T--SBS_NY--mainlight.pdf);" data-stellar-background-ratio="0.5" style="style="margin-top:-30px""> | ||
| + | <div class="overlay"></div> | ||
| + | <div class="container"> | ||
| + | <div class="row"> | ||
| + | <div class="col-md-12 text-center"> | ||
| + | <div class="display-t js-fullheight"> | ||
| + | <div class="display-tc js-fullheight animate-box" data-animate-effect="fadeIn"> | ||
| + | <h1> Modeling<em> our </em>Project</h1> | ||
| + | </div> | ||
| + | </div> | ||
| + | </div> | ||
| + | </div> | ||
| + | </div> | ||
| + | </header> | ||
| − | <div id="doc" class="markdown-body container-fluid comment-enabled" data-hard-breaks="true" style="position: relative; color: #f7f7f7; margin-top: 2.3cm;;"><h1 id="Modeling" style="color: #f7f7f7;"><a class="anchor hidden-xs" href="#Modeling" title="Modeling" style="color: #f7f7f7;"><span class="octicon octicon-link"></span></a>Modeling</h1> | + | <div id="doc" class="markdown-body container-fluid comment-enabled" data-hard-breaks="true" style="position: relative; color: #f7f7f7; margin-top: 2.3cm;;"><h1 id="Modeling" style="color: #f7f7f7;"><a class="anchor hidden-xs" href="#Modeling" title="Modeling" style="color: #f7f7f7;"><span class="octicon octicon-link"></span></a>Modeling our Project</h1> |
| − | <h2 id="Basic-Idea" style="color: #f7f7f7;"><a class="anchor hidden-xs" href="#Basic-Idea" title="Basic-Idea"><span class="octicon octicon-link"></span></a>Basic Idea</h2> | + | <h2 id="Basic-Idea" style="color: #f7f7f7; font-size: 140%;"><a class="anchor hidden-xs" href="#Basic-Idea" title="Basic-Idea"><span class="octicon octicon-link"></span></a>Basic Idea</h2> |
<p>Here we provide a simple but efficient mathematical model to interpret our data. In this model, we tried to exclude the influence of cross interactions and get high accurate result.</p> | <p>Here we provide a simple but efficient mathematical model to interpret our data. In this model, we tried to exclude the influence of cross interactions and get high accurate result.</p> | ||
| − | <h2 id="Assumption" style="color: #f7f7f7;"><a class="anchor hidden-xs" href="#Assumption" title="Assumption" style="color: #f7f7f7;"><span class="octicon octicon-link"></span></a>Assumption</h2> | + | <h2 id="Assumption" style="color: #f7f7f7; font-size: 140%;" ><a class="anchor hidden-xs" href="#Assumption" title="Assumption" style="color: #f7f7f7;"><span class="octicon octicon-link"></span></a>Assumption</h2> |
<p>Any abstract model requires proper assumptions to approximate real system. Here are our basic assumptions.</p><ol> | <p>Any abstract model requires proper assumptions to approximate real system. Here are our basic assumptions.</p><ol> | ||
| Line 200: | Line 209: | ||
I_i=\sum_j(I_{max})_{ij}H_{ij}(x_j) | I_i=\sum_j(I_{max})_{ij}H_{ij}(x_j) | ||
</script></span><br> | </script></span><br> | ||
| − | Here <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-19-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-334" style="width: 1.81em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.501em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.346em, 1001.501em, 2.585em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-335"><span class="msubsup" id="MathJax-Span-336"><span style="display: inline-block; position: relative; width: 1.449em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.932em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-337" style="font-family: STIXGeneral-Italic;">𝐻<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.881em;"><span class="texatom" id="MathJax-Span-338"><span class="mrow" id="MathJax-Span-339"><span class="mi" id="MathJax-Span-340" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span class="mi" id="MathJax-Span-341" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.372em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.253em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub></math></span></span><script type="math/tex" id="MathJax-Element-19">H_{ij}</script></span> is the responsing curve of <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-20-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>j</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-342" style="width: 1.346em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.087em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.191em, 1001.087em, 2.534em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-343"><span class="msubsup" id="MathJax-Span-344"><span style="display: inline-block; position: relative; width: 1.036em; height: 0px;"><span style="position: absolute; clip: rect(3.205em, 1000.364em, 4.341em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-345" style="font-family: STIXGeneral-Italic;">𝑗</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -4.336em; left: 0.416em;"><span class="texatom" id="MathJax-Span-346"><span class="mrow" id="MathJax-Span-347"><span class="mi" id="MathJax-Span-348" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑡</span><span class="mi" id="MathJax-Span-349" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">ℎ</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.309em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.378em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>j</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math></span></span><script type="math/tex" id="MathJax-Element-20">j^{th}</script></span> ion to <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-21-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-350" style="width: 1.294em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.036em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.191em, 1001.036em, 2.327em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-351"><span class="msubsup" id="MathJax-Span-352"><span style="display: inline-block; position: relative; width: 0.984em; height: 0px;"><span style="position: absolute; clip: rect(3.205em, 1000.261em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-353" style="font-family: STIXGeneral-Italic;">𝑖</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -4.336em; left: 0.313em;"><span class="texatom" id="MathJax-Span-354"><span class="mrow" id="MathJax-Span-355"><span class="mi" id="MathJax-Span-356" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑡</span><span class="mi" id="MathJax-Span-357" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">ℎ</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.059em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.128em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math></span></span><script type="math/tex" id="MathJax-Element-21">i^{th}</script></span> vector.</p><h2 id="Model"style="color: #f7f7f7;"><a class="anchor hidden-xs" href="#Model" title="Model"><span class="octicon octicon-link"></span></a>Model</h2><p>Based on the analysis in assumption, the model system is clear. For <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-22-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-358" style="width: 1.294em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.036em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.191em, 1001.036em, 2.327em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-359"><span class="msubsup" id="MathJax-Span-360"><span style="display: inline-block; position: relative; width: 0.984em; height: 0px;"><span style="position: absolute; clip: rect(3.205em, 1000.261em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-361" style="font-family: STIXGeneral-Italic;">𝑖</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -4.336em; left: 0.313em;"><span class="texatom" id="MathJax-Span-362"><span class="mrow" id="MathJax-Span-363"><span class="mi" id="MathJax-Span-364" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑡</span><span class="mi" id="MathJax-Span-365" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">ℎ</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.059em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.128em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math></span></span><script type="math/tex" id="MathJax-Element-22">i^{th}</script></span> the detector should be:<br> | + | Here <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-19-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-334" style="width: 1.81em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.501em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.346em, 1001.501em, 2.585em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-335"><span class="msubsup" id="MathJax-Span-336"><span style="display: inline-block; position: relative; width: 1.449em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.932em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-337" style="font-family: STIXGeneral-Italic;">𝐻<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.881em;"><span class="texatom" id="MathJax-Span-338"><span class="mrow" id="MathJax-Span-339"><span class="mi" id="MathJax-Span-340" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span class="mi" id="MathJax-Span-341" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.372em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.253em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub></math></span></span><script type="math/tex" id="MathJax-Element-19">H_{ij}</script></span> is the responsing curve of <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-20-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>j</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-342" style="width: 1.346em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.087em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.191em, 1001.087em, 2.534em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-343"><span class="msubsup" id="MathJax-Span-344"><span style="display: inline-block; position: relative; width: 1.036em; height: 0px;"><span style="position: absolute; clip: rect(3.205em, 1000.364em, 4.341em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-345" style="font-family: STIXGeneral-Italic;">𝑗</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -4.336em; left: 0.416em;"><span class="texatom" id="MathJax-Span-346"><span class="mrow" id="MathJax-Span-347"><span class="mi" id="MathJax-Span-348" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑡</span><span class="mi" id="MathJax-Span-349" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">ℎ</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.309em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.378em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>j</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math></span></span><script type="math/tex" id="MathJax-Element-20">j^{th}</script></span> ion to <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-21-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-350" style="width: 1.294em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.036em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.191em, 1001.036em, 2.327em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-351"><span class="msubsup" id="MathJax-Span-352"><span style="display: inline-block; position: relative; width: 0.984em; height: 0px;"><span style="position: absolute; clip: rect(3.205em, 1000.261em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-353" style="font-family: STIXGeneral-Italic;">𝑖</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -4.336em; left: 0.313em;"><span class="texatom" id="MathJax-Span-354"><span class="mrow" id="MathJax-Span-355"><span class="mi" id="MathJax-Span-356" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑡</span><span class="mi" id="MathJax-Span-357" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">ℎ</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.059em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.128em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math></span></span><script type="math/tex" id="MathJax-Element-21">i^{th}</script></span> vector.</p> |
| + | <h2 id="Model"style="color: #f7f7f7; font-size: 140%;"><a class="anchor hidden-xs" href="#Model" title="Model"><span class="octicon octicon-link"></span></a>Model</h2> | ||
| + | <p>Based on the analysis in assumption, the model system is clear. For <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><span class="MathJax" id="MathJax-Element-22-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math>" role="presentation" style="position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-358" style="width: 1.294em; display: inline-block;"><span style="display: inline-block; position: relative; width: 1.036em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.191em, 1001.036em, 2.327em, -999.997em); top: -2.167em; left: 0em;"><span class="mrow" id="MathJax-Span-359"><span class="msubsup" id="MathJax-Span-360"><span style="display: inline-block; position: relative; width: 0.984em; height: 0px;"><span style="position: absolute; clip: rect(3.205em, 1000.261em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-361" style="font-family: STIXGeneral-Italic;">𝑖</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -4.336em; left: 0.313em;"><span class="texatom" id="MathJax-Span-362"><span class="mrow" id="MathJax-Span-363"><span class="mi" id="MathJax-Span-364" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑡</span><span class="mi" id="MathJax-Span-365" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">ℎ</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span></span><span style="display: inline-block; width: 0px; height: 2.172em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -0.059em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 1.128em;"></span></span></nobr><span class="MJX_Assistive_MathML" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mi>i</mi><mrow class="MJX-TeXAtom-ORD"><mi>t</mi><mi>h</mi></mrow></msup></math></span></span><script type="math/tex" id="MathJax-Element-22">i^{th}</script></span> the detector should be:<br> | ||
<span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><div class="MathJax_Display" style="text-align: center;"><span class="MathJax" id="MathJax-Element-23-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"><msub><mi>I</mi><mi>i</mi></msub><mo>=</mo><munder><mo>&#x2211;</mo><mi>j</mi></munder><mo stretchy="false">(</mo><msub><mi>I</mi><mrow class="MJX-TeXAtom-ORD"><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><msub><mo stretchy="false">)</mo><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>x</mi><mi>j</mi></msub><mo stretchy="false">)</mo></math>" role="presentation" style="text-align: center; position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-366" style="width: 11.418em; display: inline-block;"><span style="display: inline-block; position: relative; width: 9.403em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.346em, 1009.352em, 4.032em, -999.997em); top: -2.477em; left: 0em;"><span class="mrow" id="MathJax-Span-367"><span class="msubsup" id="MathJax-Span-368"><span style="display: inline-block; position: relative; width: 0.777em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.519em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-369" style="font-family: STIXGeneral-Italic;">𝐼<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.467em;"><span class="mi" id="MathJax-Span-370" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-371" style="font-family: STIXGeneral-Regular; padding-left: 0.313em;">=</span><span class="munderover" id="MathJax-Span-372" style="padding-left: 0.313em;"><span style="display: inline-block; position: relative; width: 1.294em; height: 0px;"><span style="position: absolute; clip: rect(2.843em, 1001.191em, 4.6em, -999.997em); top: -3.975em; left: 0em;"><span class="mo" id="MathJax-Span-373" style="font-family: STIXSizeOneSym; vertical-align: -0.514em;">∑</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; clip: rect(3.412em, 1000.261em, 4.341em, -999.997em); top: -2.838em; left: 0.519em;"><span class="mi" id="MathJax-Span-374" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-375" style="font-family: STIXGeneral-Regular;">(</span><span class="msubsup" id="MathJax-Span-376"><span style="display: inline-block; position: relative; width: 1.81em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.519em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-377" style="font-family: STIXGeneral-Italic;">𝐼<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.467em;"><span class="texatom" id="MathJax-Span-378"><span class="mrow" id="MathJax-Span-379"><span class="mi" id="MathJax-Span-380" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑚</span><span class="mi" id="MathJax-Span-381" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑎</span><span class="mi" id="MathJax-Span-382" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑥</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="msubsup" id="MathJax-Span-383"><span style="display: inline-block; position: relative; width: 0.881em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.313em, 4.29em, -999.997em); top: -3.975em; left: 0em;"><span class="mo" id="MathJax-Span-384" style="font-family: STIXGeneral-Regular;">)</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.313em;"><span class="texatom" id="MathJax-Span-385"><span class="mrow" id="MathJax-Span-386"><span class="mi" id="MathJax-Span-387" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span class="mi" id="MathJax-Span-388" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="msubsup" id="MathJax-Span-389"><span style="display: inline-block; position: relative; width: 1.449em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.932em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-390" style="font-family: STIXGeneral-Italic;">𝐻<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.881em;"><span class="texatom" id="MathJax-Span-391"><span class="mrow" id="MathJax-Span-392"><span class="mi" id="MathJax-Span-393" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span class="mi" id="MathJax-Span-394" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-395" style="font-family: STIXGeneral-Regular;">(</span><span class="msubsup" id="MathJax-Span-396"><span style="display: inline-block; position: relative; width: 0.881em; height: 0px;"><span style="position: absolute; clip: rect(3.36em, 1000.519em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-397" style="font-family: STIXGeneral-Italic;">𝑥</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.571em;"><span class="mi" id="MathJax-Span-398" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-399" style="font-family: STIXGeneral-Regular;">)</span></span><span style="display: inline-block; width: 0px; height: 2.482em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -1.747em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 3.003em;"></span></span></nobr><span class="MJX_Assistive_MathML MJX_Assistive_MathML_Block" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML" display="block"><msub><mi>I</mi><mi>i</mi></msub><mo>=</mo><munder><mo>∑</mo><mi>j</mi></munder><mo stretchy="false">(</mo><msub><mi>I</mi><mrow class="MJX-TeXAtom-ORD"><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><msub><mo stretchy="false">)</mo><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>x</mi><mi>j</mi></msub><mo stretchy="false">)</mo></math></span></span></div><script type="math/tex; mode=display" id="MathJax-Element-23"> | <span class="mathjax"><span class="MathJax_Preview" style="color: inherit;"></span><div class="MathJax_Display" style="text-align: center;"><span class="MathJax" id="MathJax-Element-23-Frame" tabindex="0" data-mathml="<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"><msub><mi>I</mi><mi>i</mi></msub><mo>=</mo><munder><mo>&#x2211;</mo><mi>j</mi></munder><mo stretchy="false">(</mo><msub><mi>I</mi><mrow class="MJX-TeXAtom-ORD"><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><msub><mo stretchy="false">)</mo><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>x</mi><mi>j</mi></msub><mo stretchy="false">)</mo></math>" role="presentation" style="text-align: center; position: relative;"><nobr aria-hidden="true"><span class="math" id="MathJax-Span-366" style="width: 11.418em; display: inline-block;"><span style="display: inline-block; position: relative; width: 9.403em; height: 0px; font-size: 121%;"><span style="position: absolute; clip: rect(1.346em, 1009.352em, 4.032em, -999.997em); top: -2.477em; left: 0em;"><span class="mrow" id="MathJax-Span-367"><span class="msubsup" id="MathJax-Span-368"><span style="display: inline-block; position: relative; width: 0.777em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.519em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-369" style="font-family: STIXGeneral-Italic;">𝐼<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.467em;"><span class="mi" id="MathJax-Span-370" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-371" style="font-family: STIXGeneral-Regular; padding-left: 0.313em;">=</span><span class="munderover" id="MathJax-Span-372" style="padding-left: 0.313em;"><span style="display: inline-block; position: relative; width: 1.294em; height: 0px;"><span style="position: absolute; clip: rect(2.843em, 1001.191em, 4.6em, -999.997em); top: -3.975em; left: 0em;"><span class="mo" id="MathJax-Span-373" style="font-family: STIXSizeOneSym; vertical-align: -0.514em;">∑</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; clip: rect(3.412em, 1000.261em, 4.341em, -999.997em); top: -2.838em; left: 0.519em;"><span class="mi" id="MathJax-Span-374" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-375" style="font-family: STIXGeneral-Regular;">(</span><span class="msubsup" id="MathJax-Span-376"><span style="display: inline-block; position: relative; width: 1.81em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.519em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-377" style="font-family: STIXGeneral-Italic;">𝐼<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.467em;"><span class="texatom" id="MathJax-Span-378"><span class="mrow" id="MathJax-Span-379"><span class="mi" id="MathJax-Span-380" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑚</span><span class="mi" id="MathJax-Span-381" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑎</span><span class="mi" id="MathJax-Span-382" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑥</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="msubsup" id="MathJax-Span-383"><span style="display: inline-block; position: relative; width: 0.881em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.313em, 4.29em, -999.997em); top: -3.975em; left: 0em;"><span class="mo" id="MathJax-Span-384" style="font-family: STIXGeneral-Regular;">)</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.313em;"><span class="texatom" id="MathJax-Span-385"><span class="mrow" id="MathJax-Span-386"><span class="mi" id="MathJax-Span-387" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span class="mi" id="MathJax-Span-388" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="msubsup" id="MathJax-Span-389"><span style="display: inline-block; position: relative; width: 1.449em; height: 0px;"><span style="position: absolute; clip: rect(3.153em, 1000.932em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-390" style="font-family: STIXGeneral-Italic;">𝐻<span style="display: inline-block; overflow: hidden; height: 1px; width: 0.054em;"></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.881em;"><span class="texatom" id="MathJax-Span-391"><span class="mrow" id="MathJax-Span-392"><span class="mi" id="MathJax-Span-393" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑖</span><span class="mi" id="MathJax-Span-394" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span></span></span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-395" style="font-family: STIXGeneral-Regular;">(</span><span class="msubsup" id="MathJax-Span-396"><span style="display: inline-block; position: relative; width: 0.881em; height: 0px;"><span style="position: absolute; clip: rect(3.36em, 1000.519em, 4.135em, -999.997em); top: -3.975em; left: 0em;"><span class="mi" id="MathJax-Span-397" style="font-family: STIXGeneral-Italic;">𝑥</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span><span style="position: absolute; top: -3.82em; left: 0.571em;"><span class="mi" id="MathJax-Span-398" style="font-size: 70.7%; font-family: STIXGeneral-Italic;">𝑗</span><span style="display: inline-block; width: 0px; height: 3.98em;"></span></span></span></span><span class="mo" id="MathJax-Span-399" style="font-family: STIXGeneral-Regular;">)</span></span><span style="display: inline-block; width: 0px; height: 2.482em;"></span></span></span><span style="display: inline-block; overflow: hidden; vertical-align: -1.747em; border-left-width: 0px; border-left-style: solid; width: 0px; height: 3.003em;"></span></span></nobr><span class="MJX_Assistive_MathML MJX_Assistive_MathML_Block" role="presentation"><math xmlns="http://www.w3.org/1998/Math/MathML" display="block"><msub><mi>I</mi><mi>i</mi></msub><mo>=</mo><munder><mo>∑</mo><mi>j</mi></munder><mo stretchy="false">(</mo><msub><mi>I</mi><mrow class="MJX-TeXAtom-ORD"><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><msub><mo stretchy="false">)</mo><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>H</mi><mrow class="MJX-TeXAtom-ORD"><mi>i</mi><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>x</mi><mi>j</mi></msub><mo stretchy="false">)</mo></math></span></span></div><script type="math/tex; mode=display" id="MathJax-Element-23"> | ||
I_i=\sum_j(I_{max})_{ij}H_{ij}(x_j) | I_i=\sum_j(I_{max})_{ij}H_{ij}(x_j) | ||
| Line 450: | Line 461: | ||
</script> | </script> | ||
| + | |||
| + | <footer id="fh5co-footer" role="contentinfo" class="fh5co-section"> | ||
| + | <div class="container"> | ||
| + | <div class="row row-pb-md"> | ||
| + | <div class="col-md-4 fh5co-widget"> | ||
| + | <h4>Our Goals</h4> | ||
| + | <p>We aim to devise a way to use our MerR-like heavy metal biosensors to detect heavy metal, which will address the issue of soil/water pollution, especially under urgent conditions.</p> | ||
| + | </div> | ||
| + | <div class="col-md-2 col-md-push-1 fh5co-widget"> | ||
| + | <h4>Our Project</h4> | ||
| + | <ul class="fh5co-footer-links"> | ||
| + | <li><a href="index.html">Home</a></li> | ||
| + | <li><a href="Description">Project</a></li> | ||
| + | <li><a href="Model">Modeling</a></li> | ||
| + | <li><a href="Hardware">Device & Software</a></li> | ||
| + | </ul> | ||
| + | </div> | ||
| + | |||
| + | <div class="col-md-2 col-md-push-1 fh5co-widget"> | ||
| + | <h4>Human Practice</h4> | ||
| + | <ul class="fh5co-footer-links"> | ||
| + | <li><a href="Human_Practices">Human Practices</a></li> | ||
| + | <li><a href="Collaborations">Collaborations</a></li> | ||
| + | </ul> | ||
| + | </div> | ||
| + | |||
| + | <div class="col-md-4 col-md-push-1 fh5co-widget"> | ||
| + | <h4>Contact Information</h4> | ||
| + | <ul class="fh5co-footer-links"> | ||
| + | <li>1 Chapman Pkwy, <br> Stony Brook, New York 11790</li> | ||
| + | <li><a href="mailto:jiacheng.yin@sbs.org">jiacheng.yin@sbs.org</a></li> | ||
| + | </ul> | ||
| + | </div> | ||
| + | |||
| + | </div> | ||
| + | |||
| + | <div class="row copyright"> | ||
| + | <div class="col-md-12 text-center"> | ||
| + | <p> | ||
| + | <small class="block">2019 Stony Brook School iGEM team</small> | ||
| + | <small class="block">Template Designed by <a href="http://freehtml5.co/" target="_blank">FreeHTML5.co</a></small> | ||
| + | </p> | ||
| + | </div> | ||
| + | </div> | ||
| + | |||
| + | </div> | ||
| + | </footer> | ||
<div class="gototop js-top"> | <div class="gototop js-top"> | ||
Latest revision as of 01:45, 22 October 2019
<!DOCTYPE html>
Modeling our Project
Basic Idea
Here we provide a simple but efficient mathematical model to interpret our data. In this model, we tried to exclude the influence of cross interactions and get high accurate result.
Assumption
Any abstract model requires proper assumptions to approximate real system. Here are our basic assumptions.
The first assumption is natural. Consider the whole responsing process as a computation module which consists of input, processing layer, and output. Therefore, the input is the number of ions. Processing layer is to produce fluorence protein when the ions successfully bound with vector. Output is the number of fluorence protein. The only thing affect the output is the probability of ion binding, which is determined by input. In average, we have the following:
𝐼∝𝐻(𝑥) 𝐼 refers to intensity and 𝐻(𝑥) refers to the probability of binding. 𝑥 stands for the concentration of ion.
Here
The second assumption helps us to calculate the probability in assumption one, which is based on a simplified but well-established biochemical model. In short, Hill equation reflects the binding of ligands to macromolecules, as a function of the ligand concentration. It can be written as:
𝑃𝑎𝑐𝑡𝑖𝑣𝑒=𝐻𝑎(𝑥)=𝑥𝑛𝑘𝑛𝑎+𝑥𝑛𝑃𝑖𝑛ℎ𝑖𝑏𝑖𝑡𝑖𝑣𝑒=𝐻𝑖(𝑥)=𝑘𝑛𝑖𝑘𝑛𝑖+𝑥𝑛
To explain the physical meaning of those parameters, here we briefly introduce the model. Suppose an “on” state of processing layer requires𝑛 binding ions, and corresponding chemical reaction can be written as:
𝑃+𝑛𝐼⇌𝑛𝑃𝐼 𝐾𝑑=[𝑃][𝐼]𝑛[𝑛𝑃𝐼] [𝑛𝑃𝐼] over all states. To simplify, we only consider two states: all binding([𝑃] ) and no binding([𝑛𝑃𝐼] ).Therefore the probability is:
𝑃=[𝑛𝑃𝐼][𝑃]+[𝑛𝑃𝐼]=[𝑃][𝐼]𝑛𝐾𝑑[𝑃]+[𝑃][𝐼]𝑛𝐾𝑑=[𝐼]𝑛𝐾𝑑+[𝐼]𝑛=[𝐼]𝑛𝑘𝑛𝑎+[𝐼]𝑛 𝐾𝑑 is a positive constant, it can be rewritten in power form:𝐾𝑑=𝑘𝑛𝑎 . So the probability is the function of ion concentration. This is exactly the second assumption.
And the dissociation constant can be expressed as:
What we focus on is the probability of binding, which can be calculated by the proportion of
Since
By integrating assumption 1 and 2, it is easy to find out the function of intensity:
𝐼=𝐼𝑚𝑎𝑥[𝐼]𝑛𝑘𝑛𝑎+[𝐼]𝑛
𝐼𝑚𝑎𝑥 is the maximum intensity the system can get.
The third assumption ensures a linear model system. Since each kind of ion has its unique responsing curve to each detector, the final intensity will be the addition of intensity from all kinds of ion. This assumption ensure that the intensity for each ion will only depend on the concentration of their own and will not be influenced by other ions. Therefore, the intensity for𝑖𝑡ℎ detector should be:
𝐼𝑖=∑𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑥𝑗) 𝐻𝑖𝑗 is the responsing curve of 𝑗𝑡ℎ ion to 𝑖𝑡ℎ vector.
Here
Model
Based on the analysis in assumption, the model system is clear. For𝑖𝑡ℎ the detector should be:
𝐼𝑖=∑𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑥𝑗) 𝐻𝑖𝑗 is the responsing curve of 𝑗𝑡ℎ ion to 𝑖𝑡ℎ vector.
Here
To analyze our data, we should first determine the all coeffcients in the expression and then use those standard functions to calculate the actual concentrations in an unknown sample.
Determine Coefficients
During the experiment, we have tested how the detectors response to ions by a concentration gradient. To couple the experiment with model, it is necessary to do some transformation on fitting equation.
Generally,𝐼𝑗=(𝐼𝑚𝑎𝑥)𝐻(𝑥) can be rewritten as:
log(𝐼𝐼𝑚𝑎𝑥−𝐼)=𝑛log(𝑘𝑎)−𝑛log(𝐼) 𝐼𝑚𝑎𝑥 in this equation because this value determine the reprocessed data of output. Another question is, due to the large scale of our data, to ease the workload of proceesing such data. To meet the needs of these two question, define the ratio between output data and the maximum of all output data as the standard output. As following shows:
In this form, we can easily get a linear relation between our input concerntration and output. The question is how to find out
The elements in𝐼′𝑜𝑢𝑡𝑝𝑢𝑡 fit following equation:
log𝐼′𝑖max𝐼𝑜𝑢𝑡𝑝𝑢𝑡𝐼𝑚𝑎𝑥−𝐼′𝑖max𝐼𝑜𝑢𝑡𝑝𝑢𝑡=𝑛log𝑥𝑖−𝑛log𝑘 𝐼𝑚𝑎𝑥max𝐼𝑜𝑢𝑡𝑝𝑢𝑡 as a parameter 𝑃𝐼𝑚𝑎𝑥 . So the equation we actually simulate is following one:
log𝑦′𝑖𝑃𝐼𝑚𝑎𝑥−𝑦′𝑖=𝑛log𝑥𝑖−𝑛log𝑘
We define the value of
Use Mathematica, the following code is shown:
outputdata = {Output1, Output2, Output3, Output4, Output5, Output6, Output7}; Processeddata = outputdata/Max[outputdata] // N; data' = {{Log10[10^(-10)], Processeddata[[1]]}, {Log10[10^(-9)], Processeddata[[2]]}, {Log10[10^(-8)], Processeddata[[3]]}, {Log10[10^(-7)], Processeddata[[4]]}, {Log10[10^(-6)], Processeddata[[5]]}, {Log10[10^(-5)], Processeddata[[6]]}, {Log10[10^(-4)], Processeddata[[7]]}}; data = {{data'[[1, 1]], data'[[1, 2]]}, {data'[[2, 1]], data'[[2, 2]]}, {data'[[3, 1]], data'[[3, 2]]}, {data'[[4, 1]], data'[[4, 2]]}, {data'[[5, 1]], data'[[5, 2]]}, {data'[[6, 1]], data'[[6, 2]]}, {data'[[7, 1]], data'[[7, 2]]}}; solu = Flatten[ Solve[Log10[(y*PImax)/(1 - (y*PImax))] == n*x - n*logk, y]]; fitparameter = (FindFit[data, y /. solu, {PImax, logk, n}, x]) fit = y /. solu /. fitparameter; Show[ListPlot[data, PlotStyle -> Red], Plot[fit, {x, -11, 0}], PlotRange -> {0, 1}]Data Analysis
In last section, we successfully got statistics of each detector. Now they will be used to analyze an unknown sample.
Denote the concentration of each ion in the sample is𝑋𝑛 .
For each detector:
𝐼𝑖=∑𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑋𝑗)𝑓𝑜𝑟 𝑖=1,2,3,···,𝑛 {(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑥)} has been determined for all 𝑖,𝑗 . {𝐼𝑖} is the output of unknown sample in 𝑖𝑡ℎ detector.
Which
Now we have𝑛 equation for 𝑛 variables, it should determined the value of all variables. But unfortunately, this is not a linear system and more importantly, the technique we used to get linear form in last section cannot be transplanted here. A general way to solve such an nonlinear system is so-called “Netwon Iteration Method”.
New Form
First rewrite the model as:
∑𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑋𝑗)−𝐼𝑖=0
𝐹(𝑋)=⎛⎝⎜⎜⎜⎜𝐹1(𝑋)𝐹2(𝑋)···𝐹𝑛(𝑋)⎞⎠⎟⎟⎟⎟;𝑋=⎛⎝⎜⎜⎜⎜𝑥1𝑥2···𝑥𝑛⎞⎠⎟⎟⎟⎟𝐹𝑖(𝑋)=∑𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑥𝑗)−𝐼𝑖,𝑓𝑜𝑟 𝑖=1,2,3,···,𝑛∴𝐹(𝑋)=0
𝐽(𝑋)=⎛⎝⎜⎜⎜⎜⎜⎜∂𝐹1(𝑋)∂𝑥1∂𝐹2(𝑋)∂𝑥1···∂𝐹𝑛(𝑋)∂𝑥1∂𝐹1(𝑋)∂𝑥1∂𝐹2(𝑋)∂𝑥1···∂𝐹𝑛(𝑋)∂𝑥1············∂𝐹1(𝑋)∂𝑥𝑛∂𝐹2(𝑋)∂𝑥𝑛···∂𝐹𝑛(𝑋)∂𝑥𝑛⎞⎠⎟⎟⎟⎟⎟⎟∂𝐹𝑖(𝑋)∂𝑥𝑗=∂∂𝑥𝑗(∑𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗𝐻𝑖𝑗(𝑥𝑗)−𝐼𝑖)=(𝐼𝑚𝑎𝑥)𝑖𝑗∂∂𝑥𝑗𝐻𝑖𝑗(𝑥𝑗)
𝑥𝑛+1=𝑥𝑛−𝐽−1(𝑥𝑛)𝐹(𝑥𝑛)
lim𝑛→+∞𝑥𝑛=⎛⎝⎜⎜⎜⎜𝑋1𝑋2···𝑋𝑛⎞⎠⎟⎟⎟⎟
Define:
Now calculate the Jacobian matrix:
Once Jacobian matrix is determined, take iteration:
Theoretically we have
This is exactly the newton iteration method.
Notice
The iteration method does not always works well, which is depending on the property of iteration functions. For some nonlinear system, the approximation could be useful in a very ting neighborhood. Such neighborhood determines the flexibility of intial value choosing. If the initial value is choosen far from the solution and the solving system is bad, then the iteration will be deficient. Therefore, to choose proper position and proper inital value to start iteration algorithm is critial to final result.
How to choose initial value
To get a good initial value for iteration, a proper range of solution should be guessed. Based on the biological property of detectors, the cross talk between target ion and untarget ion should be low. In extreme, the response of a wonderful detector to untarget ion should be a constant. First get rid of all posible background, and then solve the equation only with target ion.
That means to solve𝑥0𝑖 in following:
(𝐼𝑚𝑎𝑥)𝑖𝑖𝐻𝑖𝑖(𝑥𝑖)=𝐼𝑖−∑𝑖≠𝑗(𝐼𝑚𝑎𝑥)𝑖𝑗 𝑖=1,2,···,𝑛
For all
Then we got initial value for iteration:
𝑥0=⎛⎝⎜⎜⎜⎜𝑥01𝑥02···𝑥0𝑛⎞⎠⎟⎟⎟⎟
How to choose data?
In our test, the sample can be diluted to 10X or even 100X according to its quality. Therefore, we can choose the data that benefits to the data processing. From analysis, we want to make a robust approximation, which means a small perturbation will not lead to a big change of result. So the region with large derivative should be avoided. One thing should be noticed that since the inputs in algorithm are the detector results and outputs are deduced from response curves, the derivative of the inverse functions of the curve, not response curve, should be considered. Therefore, we choose the detector results with biggest change as the “proper” data. To define the “Change”, here considers the difference between adjacent data points.
Example
Here we provide an example for 2 variables:
H[n_, k_, x_, V_] := V*x^n/(k^n + x^n)(*Hill Function*); DH[n_, k_, x_, V_] := ( V*k^n n x^(-1 + n))/(k^n + x^ n)^2(*Derivative of Hill Function*); M = {{0.003,0.0007}} (*Initial Value,Determine by "Solve[5*x^2/((10^-3)^2+x^2)\\[Equal]4.5,x];Solve[7*x^2/((10^-3.5)^2+x^2)\[Equal]6,x]"*); For[i = 1, i < 10, i++, x1 = M[[i, 1]]; x2 = M[[i, 2]]; g1 = H[2, 10^(-3), x1, 5] + H[2, 10^(-1), x2, 4] - 4.5; g2 = H[2, 10^(-2), x1, 3] + H[2, 10^(-3.5), x2, 7] - 6; {x1d, y1d} = {x1, x2} - Inverse[{{DH[2, 10^(-3), x1, 5], DH[2, 10^(-1), x2, 4]}, {DH[2, 10^(-2), x1, 3], DH[2, 10^(-3.5), x2, 7]}}].{g1, g2}; AppendTo[M, {x1d, y1d}]](*Newton Method*); Print[M[[10]]]Result:{0.00299939,0.000679022}
Compared to original:{0.003,0.0007}
Good approximation.