Team:Paris/Modeling/Implementation

From 2008.igem.org

(Difference between revisions)
(Parameters Finder Programs)
(Parameters Finder for our Example)
Line 18: Line 18:
We just show hereby the annoted program ''' find_FP.m ''' that is used to estimate, for instance, the parameters in :
We just show hereby the annoted program ''' find_FP.m ''' that is used to estimate, for instance, the parameters in :
-
* <span style="color:#0000FF;">&#131;5( [''FlhDC''], 0 ) = ''&beta;<sub>24</sub> * &#131;<sub>hill</sub>''( [''FlhDC''], ''K<sub>2</sub>'', ''n<sub>2</sub>'' )</span> and  
+
* '''&#131;5( [''FlhDC''], 0 ) = ''&beta;<sub>24</sub> * &#131;<sub>hill</sub>''( [''FlhDC''], ''K<sub>2</sub>'', ''n<sub>2</sub>'' )''' and  
-
* <span style="color:#0000FF;">&#131;5( 0, [''FliA''] ) = ''&beta;<sub>25</sub> * &#131;<sub>hill</sub>''( [''FliA''], ''K<sub>8</sub>'', ''n<sub>8</sub>'' )</span>
+
* '''&#131;5( 0, [''FliA''] ) = ''&beta;<sub>25</sub> * &#131;<sub>hill</sub>''( [''FliA''], ''K<sub>8</sub>'', ''n<sub>8</sub>'' )'''
== All Algorithms ==
== All Algorithms ==

Revision as of 19:47, 29 October 2008

Implementation


We use Matlab for all implementations.

Contents

Parameters Finder Programs

the data

The experimental data consist typically in two tables, X_data (various concentrations of the transcription factor) and Y_data (corresponding output values).

  • controlling X_data : thanks to the prior characterization of the inductible promoters that control the transcription factor concentrations, we can deduce from the Inv_f1.m and Inv_f2.m functions the necessary concentrations of aTc and arabinose to introduce in the medium to get the wanted concentrations of transcription factor.
  • getting Y_data : the linear conversion between the fluorescence of GFP at maturation and its concentration gives us directly the expected data.

Parameters Finder for our Example

We just show hereby the annoted program find_FP.m that is used to estimate, for instance, the parameters in :

  • ƒ5( [FlhDC], 0 ) = β24 * ƒhill( [FlhDC], K2, n2 ) and
  • ƒ5( 0, [FliA] ) = β25 * ƒhill( [FliA], K8, n8 )

All Algorithms

↓ Prior for Characterization ↑
↓ Parameters Finders ↑
↓ The Global Model ↑


<Back - to "Workflow on an Example"|