0001 function [dbBool, noDocMets]=checkFormulaValidty(model)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
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);