Difference between revisions of "Team:UFSCar-Brasil/part1.html"

Line 34: Line 34:
 
<p>This curve obtained r² = 0.95457, what means that more than 95% of data are modeled by this function. Using the <b>MatLab</b> software to function fitting. We also find the value of equation constants (Appendix Table 2).</p>
 
<p>This curve obtained r² = 0.95457, what means that more than 95% of data are modeled by this function. Using the <b>MatLab</b> software to function fitting. We also find the value of equation constants (Appendix Table 2).</p>
 
The surface equation is:
 
The surface equation is:
<p><b>z=1.07+0.06x-0.044y-0.04x²+0.0009y²+0.003xy</b>  (Eq. 1)</p>
+
 
 +
<h4 class="ui center aligned header"><b><b>z=1.07+0.06x-0.044y-0.04x²+0.0009y²+0.003xy</b>  (Eq. 1)</h4>
 +
 
 
<p>From known curve was possible to simulate a surface that models most likely behavior of the three variables in question.</p>
 
<p>From known curve was possible to simulate a surface that models most likely behavior of the three variables in question.</p>
 
<img class="ui centered  400x400 image" src="https://static.igem.org/mediawiki/2015/1/16/UFSCar-Brasil_Fitted_surface.png">
 
<img class="ui centered  400x400 image" src="https://static.igem.org/mediawiki/2015/1/16/UFSCar-Brasil_Fitted_surface.png">

Revision as of 15:42, 16 September 2015

Plasmolysis

How long Bug Shoo could be stored?

Fitting Data

Starting from the experimental data presented in Table 1 in the appendix, a 3D graph (Figure 1) was constructed. It was used to better system understanding, which is the relationship between absorbance, time and PEG 6000 concentration.

Figure 1: Experimental plots of absorbance versus time versus PEG concentration.

This graph was used to find best-fitted surface to experimental points. To this, parameter r² was used, because it shows the percentage of points that fit in a certain curve. Thus, we observed the following equation to explain this behaviour:

z= \(z_0\)+ax+by+cx²+dy²+fxy.

This curve obtained r² = 0.95457, what means that more than 95% of data are modeled by this function. Using the MatLab software to function fitting. We also find the value of equation constants (Appendix Table 2).

The surface equation is:

z=1.07+0.06x-0.044y-0.04x²+0.0009y²+0.003xy (Eq. 1)

From known curve was possible to simulate a surface that models most likely behavior of the three variables in question.

Figure 2: Fitted surface obtained using the expression (1).

Results and Discussion

With the ultimate goal to find the main points of the system, meaning those points where his plane tangent in the surface has gradient value null. But first, we analyzed (using the Hessian matrix) on these important points to make possible to know whether these are the maximum points, minimum or saddle point of the function.

The calculation of the Hessian matrix is given by:

$$ H[z(x,y)] = \begin{bmatrix} \frac{\partial^2 z}{\partial x^2} & \frac{\partial^2 z}{\partial z,\partial y} \\ \frac{\partial^2 z}{\partial z,\partial y} & \frac{\partial^2 z}{\partial y^2} \end{bmatrix} $$

thus to find the value of the Hessian of z(x,y) we obtained the partial derivatives of the second order, whose values are below:

(3) $$ \frac{\partial^2 z}{\partial x^2} = 2c = - 0,08774 $$ (4) $$ \frac{\partial^2 z}{\partial y^2} = 2d = - 0,00178 $$ (5) $$ \frac{\partial^2 z}{\partial x,\partial y} = \frac{\partial^2 z}{\partial y,\partial x} = f = 0,00319 $$

with the known partial derivatives, we have the Hessian:

(6) $$ H[z(x,y)]=(-0,08774)(0,00178)-(0,00319)^2=-0,0001663533 $$

As H[f(x,y)]<0, we know that in this function there is a saddle point and there are no minimum or maximum values (Guidorizzi, 1997). Saddle point is understood as the point where the slope is zero, but it is not a maximum or minimum value. In it there is the highest elevation in one direction and lowest in the perpendicular direction.

Now, known the point we want to find, we take the partial derivatives of the first order function and made them equal to zero. To know:

(7) $$ \frac{dz}{dx}=a+2cx+fy=0 $$ (8) $$ \frac{dz}{dy}=b+2dy+fx=0 $$

Therefore:

(9) $$ -\frac{a+2cx}{f}=y=\frac{-b+fx}{2d} $$

replacing the values that are found in Table 2 of the appendix, we have that x=(2da-bf)/(f²-4cd)=1,5 and thus y=22.2 and z=0.63 . Then it is possible to concluded that for optimal absorbance point (z = 0.63), we know that it is necessary 22.2% PEG where it gets 1.5 weeks of plasmolysis.

Graph 3: Extrapolated fitted surface for extreme values of PEG concentrations and time..
Is possible to observe that for very long time values the absorbance decrease for zero very fast and the time do not pass longer than 8 weeks. The PEG concentrations also increase the absorbance values for both extreme concentrations, when approximate of zero and pass of 30% of PEG the absorbance value explodes.

Conclusions

Through these analyzes we could first find one function that suited in r² = 95%, which is very satisfactory, and validates our whole discussion about the obtained surface.

Then, through the Hessian we found a saddle point, which is also particularly relevant biologically. Because the absorbance measures the metabolic activity of bacteria, and how we want that they are plasmolysed, this value should be low, but not zero, because if so, they would be dead and would not serve for our propose. As it is a saddle point, we know that it is not null. Finally, the point was found through gradient of the function and we concluded that for the minimum value of absorbance the best PEG concentration is 22,2%.

Appendix


List of tables


Table 1: Experimental dates obtained in microbiology lab

Our amazing sponsors!

Our amazing sponsors!