findSubSysGen lists the subsystem that a reaction occurs in encoded by a gene. Returns list of subsystems. If multiple reactions are associated with gene, subsystem of first occurance will be listed. [GenSubSystem] = findSubSysGen(model) INPUT model COBRA model structure OUTPUT GenSubSystem array listing genes and subsystmes Imes Thiele 10/09
0001 function [GenSubSystem] = findSubSysGen(model) 0002 %findSubSysGen lists the subsystem that a reaction occurs in encoded by a 0003 %gene. Returns list of subsystems. If multiple reactions are associated 0004 %with gene, subsystem of first occurance will be listed. 0005 % 0006 % [GenSubSystem] = findSubSysGen(model) 0007 % 0008 %INPUT 0009 % model COBRA model structure 0010 % 0011 %OUTPUT 0012 % GenSubSystem array listing genes and subsystmes 0013 % 0014 % Imes Thiele 10/09 0015 0016 GenSubSystem(:,1) = model.genes; 0017 for i = 1 : length(model.genes) 0018 tmp = find(model.rxnGeneMat(:,i)); 0019 GenSubSystem(i,2) = model.subSystems(tmp(1)); 0020 end