Team:Alberta/Modeling/DualDeletionTest

From 2009.igem.org

(Difference between revisions)
(New page: % Remember to run start.m; b=0.01; deltamodel = model; solution1 = optimizeCbModel(deltamodel, 'max', 'false', 'false'); for a=1:length(model.rxns); deltamodel = model; deltamodel = change...)
Line 1: Line 1:
-
% Remember to run start.m;
+
% Remember to run start.m;
-
b=0.01;
+
b=0.01;
-
deltamodel = model;
+
    deltamodel = model;
-
solution1 = optimizeCbModel(deltamodel, 'max', 'false', 'false');
+
    solution1 = optimizeCbModel(deltamodel, 'max', 'false', 'false');
-
for a=1:length(model.rxns);
+
    for a=1:length(model.rxns);
-
deltamodel = model;
+
    deltamodel = model;
deltamodel = changeRxnBounds(model, model.rxns(a), 0, 'b');
deltamodel = changeRxnBounds(model, model.rxns(a), 0, 'b');
deltasolution = optimizeCbModel(deltamodel);
deltasolution = optimizeCbModel(deltamodel);

Revision as of 20:38, 17 October 2009

% Remember to run start.m;
b=0.01;
    deltamodel = model;
    solution1 = optimizeCbModel(deltamodel, 'max', 'false', 'false');
    for a=1:length(model.rxns);
    deltamodel = model;

deltamodel = changeRxnBounds(model, model.rxns(a), 0, 'b'); deltasolution = optimizeCbModel(deltamodel);

   if deltasolution.f<b;
       v(a)=model.grRules(a);
   end

end v = v'; % Copy v and model.rxns into a txt file to compare resulting data