findBlockedReaction

PURPOSE ^

findBlockedReaction determines those reactions which cannot carry any

SYNOPSIS ^

function blockedReactions = findBlockedReaction(model)

DESCRIPTION ^

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

Generated on Thu 21-Jun-2012 15:39:23 by m2html © 2003