findOrphanRxns find all orphan reactions in model (reactions with no associated genes), not including exchange rxns orphans = findOrphanRxns(model) INPUT model a COBRA model with GPRs OUTPUT orphans all orphan reactions in the model Jeff Orth 4/15/09
0001 function orphans = findOrphanRxns(model) 0002 %findOrphanRxns find all orphan reactions in model (reactions with no 0003 %associated genes), not including exchange rxns 0004 % 0005 % orphans = findOrphanRxns(model) 0006 % 0007 %INPUT 0008 % model a COBRA model with GPRs 0009 % 0010 %OUTPUT 0011 % orphans all orphan reactions in the model 0012 % 0013 % Jeff Orth 4/15/09 0014 0015 rxns = find(strcmp('',model.grRules)); 0016 [selExc,selUpt] = findExcRxns(model,true,false); 0017 orphans = model.rxns(setdiff(rxns,find(selExc))); 0018 0019 0020