Team:Warwick/Modelling/Toolbox

From 2014.igem.org

(Difference between revisions)
 
(8 intermediate revisions not shown)
Line 36: Line 36:
         <!--CONTENT START-->
         <!--CONTENT START-->
 +
       
 +
<a href="https://2014.igem.org/Main_Page"> <img class = "headerImage" style = "width: 12%;"
 +
src="https://static.igem.org/mediawiki/2014/2/22/Logo2014v2.png"> </a>
          
          
         <a href="/Team:Warwick"> <img class = "headerImage" style = "width: 30%;" src="https://static.igem.org/mediawiki/2014/f/ff/RepliconLogoON.png"> </a>
         <a href="/Team:Warwick"> <img class = "headerImage" style = "width: 30%;" src="https://static.igem.org/mediawiki/2014/f/ff/RepliconLogoON.png"> </a>
-
         <!--<a href="/Team:Warwick/DMB">--><a href="/Team:Warwick"> <img class = "headerImage" style = "width: 12%;"  
+
         <a href="/Team:Warwick"> <img class = "headerImage" style = "width: 12%;"  
src="https://static.igem.org/mediawiki/2014/f/f6/Warwick_logo.png"> </a>
src="https://static.igem.org/mediawiki/2014/f/f6/Warwick_logo.png"> </a>
Line 52: Line 55:
             <li> <a href = "/Team:Warwick/Interlab"> INTERLAB </a> </li>
             <li> <a href = "/Team:Warwick/Interlab"> INTERLAB </a> </li>
             <li> <a href = "/Team:Warwick/Attributions"> ATTRIBUTIONS </a> </li>
             <li> <a href = "/Team:Warwick/Attributions"> ATTRIBUTIONS </a> </li>
-
             <li> <a href = "/Team:Warwick/Judging"> JUDGING </a> </li>
+
              
         </div>
         </div>
<div id="secondaryMenu">
<div id="secondaryMenu">
<li> <a href = "/Team:Warwick/Modelling/Experiments"> EXPERIMENTS </a> </li>
<li> <a href = "/Team:Warwick/Modelling/Experiments"> EXPERIMENTS </a> </li>
-
<li> <a href = "/Team:Warwick/Modelling/System"> SYSTEM </a> </li>
+
<li> <a href = "/Team:Warwick/Modelling/System"> RESULTS </a> </li>
<li> <a href = "/Team:Warwick/Modelling/Toolbox"> <span>TOOLBOX</span> </a> </li>
<li> <a href = "/Team:Warwick/Modelling/Toolbox"> <span>TOOLBOX</span> </a> </li>
<li> <a href = "/Team:Warwick/Modelling/Ebola"> EBOLA </a> </li>
<li> <a href = "/Team:Warwick/Modelling/Ebola"> EBOLA </a> </li>
Line 73: Line 76:
<p>When executed on the raw tecan data in excel format it will create new sheets inside the file and ask for which plates were inside the technical replicate, leaving the space blank indicates no other technical replicates and moves on to the next biological replication of the experiment, leaving blank twice indicates end of plates and the program terminates. Each experiment is given a name within the program and put directly into the spreadsheet. This data is then in the perfect format for our data analyst who was able to import the spreadsheets directly and graph exactly what was needed.<br><br>
<p>When executed on the raw tecan data in excel format it will create new sheets inside the file and ask for which plates were inside the technical replicate, leaving the space blank indicates no other technical replicates and moves on to the next biological replication of the experiment, leaving blank twice indicates end of plates and the program terminates. Each experiment is given a name within the program and put directly into the spreadsheet. This data is then in the perfect format for our data analyst who was able to import the spreadsheets directly and graph exactly what was needed.<br><br>
-
Click <a href = "https://static.igem.org/mediawiki/2014/a/aa/Warwick_Protocol_booklet.pdf"> here </a> to download our TECAN Plate Reader analysis program.
+
Click <a href = "https://static.igem.org/mediawiki/2014/0/0d/Analyse.m"> here </a> to download our TECAN Plate Reader analysis program. You will also need to download <a href = "https://static.igem.org/mediawiki/2014/0/00/XlsColNum2Str.m"> this </a> function file (taken from <a href = "http://www.mathworks.co.uk/matlabcentral/fileexchange/15748-excel-column-number-to-column-name/content/xlsColNum2Str.m"> this page</a>) to be able to run it.
</p>
</p>
<p><b><ins><font size="4">DE Solvers</b></font></ins> <br>
<p><b><ins><font size="4">DE Solvers</b></font></ins> <br>
To help with the modelling we created several solvers that went above and beyond the MATLAB functions. We developed:
To help with the modelling we created several solvers that went above and beyond the MATLAB functions. We developed:
-
<br><br> - an ODE solver based on a fourth order Runge-Kutta method (RK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/a/aa/Warwick_Protocol_booklet.pdf"> here </a> to download our RK4 program.
+
<br><br> - an ODE solver based on a fourth order Runge-Kutta method (RK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/0/0d/RK4.m"> here </a> to download our RK4 program.
-
<br><br> - a Time Delay system solver using fourth order Runge-Kutta methods (TDRK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/a/aa/Warwick_Protocol_booklet.pdf"> here </a> to download our TDRK4 program.
+
<br><br> - a Time Delay system solver using fourth order Runge-Kutta methods (TDRK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/f/f2/TDRK4.m"> here </a> to download our TDRK4 program.
-
<br><br> - a Stochastic DE solver using fourth order Runge-Kutta methods (SRK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/a/aa/Warwick_Protocol_booklet.pdf"> here </a> to download our SRK4 program.
+
<br><br> - a Stochastic DE solver using fourth order Runge-Kutta methods (SRK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/3/34/SRK4.m"> here </a> to download our SRK4 program.
-
<br><br> - a Stochastic Time Delay solver using fourth order Runge-Kutta methods (STDRK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/a/aa/Warwick_Protocol_booklet.pdf"> here </a> to download our STDRK4 program.
+
<br><br> - a Stochastic Time Delay solver using fourth order Runge-Kutta methods (STDRK4).<br> Click <a href = "https://static.igem.org/mediawiki/2014/c/c4/STDRK4.m"> here </a> to download our STDRK4 program.
<br><br>
<br><br>
We hope that these solvers prove useful to future teams. It is our hope that by doing the difficult part of creating the solving mechanism for more complicated systems, future teams will be more willing to venture further afield and explore more mathematically diverse and realistic models.
We hope that these solvers prove useful to future teams. It is our hope that by doing the difficult part of creating the solving mechanism for more complicated systems, future teams will be more willing to venture further afield and explore more mathematically diverse and realistic models.

Latest revision as of 02:17, 18 October 2014

Programs



Feel free to download, edit and re-upload any of the MATLAB programs here to any website, or your own wiki, for your own purposes. All of the content here (as in all of IGEM) is open source.


TECAN Plate Reader
A serious issue that we faced whilst analysing our results was the speed and accuracy at which we could take the raw data from experiments and extract the information we needed. One of our most used pieces of equipment was the tecan plate reader. Each time a new reading came in, one of the team would laboriously sit for hours copying the raw text over to excel and then picking through sheets of numbers sorting out technical and biological replicates before graphing. Each time the same formulas would be rewritten and each time errors would be untraceable and set the whole process back to the beginning. This would have to be done with both Fluorescence and Optical Density making the task twice as long again. Seeing as one of the main constraints of this project is time we needed a solution to the problem that would improve efficiency and accuracy. Having been working on Matlab, the team soon got to grips with the built in excel read/write functions and so Analyse.m was born.

The idea for the program was to output an expected value and standard deviation over all technical replicates of an experiment, and then to take those figures and give an average over the bigger data set of all biological replicates, as well as the standard error for those. This would allow us to see how each individual batch reacted over time and to disregard outliers, but it would also give a uniformity in how the information was presented. Results would no longer be presented on the mood and whimsy of the analyser. Analyse.m does exactly that.

When executed on the raw tecan data in excel format it will create new sheets inside the file and ask for which plates were inside the technical replicate, leaving the space blank indicates no other technical replicates and moves on to the next biological replication of the experiment, leaving blank twice indicates end of plates and the program terminates. Each experiment is given a name within the program and put directly into the spreadsheet. This data is then in the perfect format for our data analyst who was able to import the spreadsheets directly and graph exactly what was needed.

Click here to download our TECAN Plate Reader analysis program. You will also need to download this function file (taken from this page) to be able to run it.

DE Solvers
To help with the modelling we created several solvers that went above and beyond the MATLAB functions. We developed:

- an ODE solver based on a fourth order Runge-Kutta method (RK4).
Click here to download our RK4 program.

- a Time Delay system solver using fourth order Runge-Kutta methods (TDRK4).
Click here to download our TDRK4 program.

- a Stochastic DE solver using fourth order Runge-Kutta methods (SRK4).
Click here to download our SRK4 program.

- a Stochastic Time Delay solver using fourth order Runge-Kutta methods (STDRK4).
Click here to download our STDRK4 program.

We hope that these solvers prove useful to future teams. It is our hope that by doing the difficult part of creating the solving mechanism for more complicated systems, future teams will be more willing to venture further afield and explore more mathematically diverse and realistic models.