findBlockedReaction determines those reactions which cannot carry any flux in the given simulation conditions. BlockedReaction = findBlockedReaction(model) INPUT model COBRA model structure OUTPUT blockedReactions List of blocked reactions Ines Thiele 02/09
0001 function blockedReactions = findBlockedReaction(model) 0002 %findBlockedReaction determines those reactions which cannot carry any 0003 %flux in the given simulation conditions. 0004 % 0005 % BlockedReaction = findBlockedReaction(model) 0006 % 0007 %INPUT 0008 % model COBRA model structure 0009 % 0010 %OUTPUT 0011 % blockedReactions List of blocked reactions 0012 % 0013 % 0014 % Ines Thiele 02/09 0015 0016 0017 tol = 1e-10; 0018 blockedReactions =[]; 0019 [minMax(:,1),minMax(:,2)] = fluxVariability(model,0); 0020 cnt = 1; 0021 for i=1:length(minMax) 0022 if (minMax(i,2) < tol && minMax(i,2) > -tol && minMax(i,1) < tol && minMax(i,1) > -tol) 0023 blockedReactions(cnt) = model.rxns(i); 0024 cnt = cnt + 1; 0025 end 0026 end