Team:Alberta/Modeling/DualDeletionTest

From 2009.igem.org

Revision as of 20:38, 17 October 2009 by Rpgguardian (Talk | contribs)
% 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