Team:Alberta/Modeling/DualDeletionTest
From 2009.igem.org
(Difference between revisions)
Rpgguardian (Talk | contribs) (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...) |
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; | + | 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