checkFormulaValidty

PURPOSE ^

Assesses whether metabolites in model are likely to be documented in

SYNOPSIS ^

function [dbBool, noDocMets]=checkFormulaValidty(model)

DESCRIPTION ^

 Assesses whether metabolites in model are likely to be documented in
 databases.
 
 INPUT
 Model structure array
 
 OUTPUTS
 dbBool    A boolian vector where the number of rows is equal to the number of
           metabolites in model. Contains a logical 1 in rows for metabolites that
           are likely to be documented in databases and a logical 0
           elsewhere.
 
 noDocMets A cell array containing a list of all metabolites in model that
           are not likely to be documented in any database.

 Ronan Fleming

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [dbBool, noDocMets]=checkFormulaValidty(model)
0002 % Assesses whether metabolites in model are likely to be documented in
0003 % databases.
0004 %
0005 % INPUT
0006 % Model structure array
0007 %
0008 % OUTPUTS
0009 % dbBool    A boolian vector where the number of rows is equal to the number of
0010 %           metabolites in model. Contains a logical 1 in rows for metabolites that
0011 %           are likely to be documented in databases and a logical 0
0012 %           elsewhere.
0013 %
0014 % noDocMets A cell array containing a list of all metabolites in model that
0015 %           are not likely to be documented in any database.
0016 %
0017 % Ronan Fleming
0018 %
0019 
0020 [nMets, nRxns] = size(model.S);
0021 
0022 dbBool = ones(nMets,1);
0023 
0024 for n = 1:nMets;
0025 
0026     if ~isempty(strfind(model.mets{n},'_hs'))
0027         dbBool(n)=0;
0028     end
0029     if ~isempty(strfind(model.mets{n},'hs_'))
0030         dbBool(n)=0;
0031     end
0032     if ~isempty(strfind(model.mets{n},'retn'))
0033         dbBool(n)=0;
0034     end
0035     if ~isempty(strfind(model.mets{n},'vitd'))
0036         dbBool(n)=0;
0037     end
0038     if ~isempty(strfind(model.mets{n},'rtotal'))
0039         dbBool(n)=0;
0040     end
0041     if ~isempty(strfind(model.mets{n},'deg'))
0042         dbBool(n)=0;
0043     end
0044     if ~isempty(strfind(model.metFormulas{n},'R'))
0045         dbBool(n)=0;
0046     end
0047     if ~isempty(strfind(model.metFormulas{n},'FULLR'))
0048         dbBool(n)=0;
0049     end
0050     if ~isempty(strfind(model.metFormulas{n},'X'))
0051         dbBool(n)=0;
0052     end
0053 
0054 end
0055 
0056 noDocMets = model.mets(~dbBool, 1);

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