"
Team:Alberta/Modeling/DualDeletionTest
From 2009.igem.org
(Difference between revisions)
Rpgguardian (Talk | contribs) |
Rpgguardian (Talk | contribs) |
||
| Line 1: | Line 1: | ||
% Remember to run start.m; | % Remember to run start.m; | ||
b=0.01; | 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'); | + | deltamodel = changeRxnBounds(model, model.rxns(a), 0, 'b'); |
| - | deltasolution = optimizeCbModel(deltamodel); | + | deltasolution = optimizeCbModel(deltamodel); |
if deltasolution.f<b; | if deltasolution.f<b; | ||
v(a)=model.grRules(a); | v(a)=model.grRules(a); | ||
end | end | ||
| - | end | + | end |
| - | v = v'; % Copy v and model.rxns into a txt file to compare resulting data | + | v = v'; % Copy v and model.rxns into a txt file to compare resulting data |
Latest 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
