printFluxVector Print a flux vector with reaction labels printFluxVector(model,fluxData,nonZeroFlag,excFlag,sortCol,fileName,headerRow,formulaFlag) INPUTS model COBRA model structure fluxData Data matrix/vector (for example, solution.x) OPTIONAL INPUTS nonZeroFlag Only print nonzero rows (Default = false) excFlag Only print exchange fluxes (Default = false) sortCol Column used for sorting (-1, none; 0, labels; >0, data columns;) (Default = -1) fileName Name of output file (Default = []) headerRow Header (Default = []) formulaFlag Print reaction formulas (Default = false) Markus Herrgard 6/9/06
0001 function printFluxVector(model,fluxData,nonZeroFlag,excFlag,sortCol,fileName,headerRow,formulaFlag) 0002 %printFluxVector Print a flux vector with reaction labels 0003 % 0004 % printFluxVector(model,fluxData,nonZeroFlag,excFlag,sortCol,fileName,headerRow,formulaFlag) 0005 % 0006 %INPUTS 0007 % model COBRA model structure 0008 % fluxData Data matrix/vector (for example, solution.x) 0009 % 0010 %OPTIONAL INPUTS 0011 % nonZeroFlag Only print nonzero rows (Default = false) 0012 % excFlag Only print exchange fluxes (Default = false) 0013 % sortCol Column used for sorting (-1, none; 0, labels; >0, data 0014 % columns;) (Default = -1) 0015 % fileName Name of output file (Default = []) 0016 % headerRow Header (Default = []) 0017 % formulaFlag Print reaction formulas (Default = false) 0018 % 0019 % Markus Herrgard 6/9/06 0020 0021 if (nargin < 3) 0022 nonZeroFlag = false; 0023 end 0024 if (nargin < 4) 0025 excFlag = false; 0026 end 0027 if (nargin < 5) 0028 sortCol = -1; 0029 end 0030 if (nargin < 6) 0031 fileName = []; 0032 end 0033 if (nargin < 7) 0034 headerRow = []; 0035 end 0036 if (nargin < 8) 0037 formulaFlag = false; 0038 end 0039 0040 if (excFlag) 0041 selExchange = findExcRxns(model,true,false); 0042 labels = model.rxns(selExchange); 0043 fluxData = fluxData(selExchange,:); 0044 else 0045 labels = model.rxns; 0046 end 0047 0048 % Add reaction formulas 0049 if (formulaFlag) 0050 rxnNames = labels; 0051 formulas = printRxnFormula(model,labels,false,false); 0052 for i = 1:length(rxnNames) 0053 labels{i,2} = formulas{i}; 0054 end 0055 end 0056 0057 printLabeledData(labels,fluxData,nonZeroFlag,sortCol,fileName,headerRow)