Team:Groningen/Modelling/Characterization

From 2009.igem.org

(Difference between revisions)
(New page: {{Team:Groningen/Modelling/Header}} {{todo|TODO: Talk about the devices we have and in what way we want to characterize them.}} ==Uptake measurements== {|class="ourtable" style="float:ri...)
(Starting to write fitting function.)
Line 54: Line 54:
Below we list all results, which have been used for fitting all necessary parameters.
Below we list all results, which have been used for fitting all necessary parameters.
-
{{todo|TODO: List results.}}
+
{{todo|TODO: List results. Take conversion from nmol/mg and mg/ml to µM and Vc/Vs into account.}}
 +
 
 +
<html>
 +
<input type="button" value="Fit" onClick="alert(computeCost(arsenicModelConstants(),experiments));"/>
 +
<script type="text/javascript" src="/Team:Groningen/Modelling/Arsenic.js?action=raw"></script>
 +
<script type="text/javascript" src="/Team:Groningen/Modelling/Model.js?action=raw"></script>
 +
<script type="text/javascript">
 +
var experiments = [{constants:{Vc:0.0073e-3,Vs:(1.1-0.0073)*1e-3,beta4:0},AsT:10e-6*1e-3,
 +
                    AsinT:{1:101.917808219178,10:394.520547945205,20:723.287671232877,
 +
                        40:1111.23287671233,60:1229.58904109589}}];
 +
function computeCost(c,e) {
 +
  var cost = 0, n = 0, x0, xt, times;
 +
  for(var i in e) {
 +
    var nc = {};
 +
    for(var a in c) nc[a] = c[a];
 +
    for(var a in e[i].constants) nc[a] = e[i].constants[a];
 +
    times = [];
 +
    for(var t in e[i].AsinT) times.push(t);
 +
    x0 = arsenicModelInitialization(nc,e[i].AsT);
 +
    xt = simulate(x0,times,function(t,d){return arsenicModelGradient(nc,d);});
 +
    for(var j in xt.time) {
 +
      if (e[i].AsinT[xt.time[j]]) {
 +
        cost += Math.pow(e[i].AsinT[xt.time[j]]-xt.AsinT[j],2);
 +
        n++;
 +
      }
 +
    }
 +
  }
 +
  return Math.sqrt(cost/n);
 +
}
 +
</script>
 +
</html>

Revision as of 13:29, 28 September 2009

[http://2009.igem.org/Team:Groningen http://2009.igem.org/wiki/images/f/f1/Igemhomelogo.png]


TODO: Talk about the devices we have and in what way we want to characterize them.

Uptake measurements

Sampling scheme
Time (min)
0 10 20 40 60
As(III)exT(0)
(µM)
0 x
10 x x x x x
20 x
50 x
100 x

To efficiently look at both time and concentration dependent processes we took samples as in the table on the right. Below we list all results, which have been used for fitting all necessary parameters.

TODO: List results. Take conversion from nmol/mg and mg/ml to µM and Vc/Vs into account.