findGeneIDs Find gene numbers in a model geneID = findGeneIDs(model,geneList) INPUTS model COBRA model structure geneList List of genes OUTPUT geneID List of gene IDs corresponding to geneList Jeff Orth 7/16/09
0001 function geneID = findGeneIDs(model,geneList) 0002 %findGeneIDs Find gene numbers in a model 0003 % 0004 % geneID = findGeneIDs(model,geneList) 0005 % 0006 %INPUTS 0007 % model COBRA model structure 0008 % geneList List of genes 0009 % 0010 %OUTPUT 0011 % geneID List of gene IDs corresponding to geneList 0012 % 0013 % Jeff Orth 7/16/09 0014 0015 if (iscell(geneList)) 0016 [tmp,geneID] = ismember(geneList,model.genes); 0017 else 0018 geneID = find(strcmp(model.genes,geneList)); 0019 if (isempty(geneID)) 0020 geneID = 0; 0021 end 0022 if (length(geneID) > 1) 0023 geneID = geneID(1); 0024 end 0025 end