printSampleStats

PURPOSE ^

printSampleStats Print out sample statistics for multiple samples

SYNOPSIS ^

function printSampleStats(samples,commonModel,sampleNames,fileName)

DESCRIPTION ^

printSampleStats Print out sample statistics for multiple samples

 printSampleStats(samples,commonModel,sampleNames,fileName)

INPUTS
 samples       Samples to plot
 commonModel   COBRA model structure
 sampleNames   Names of samples

OPTIONAL INPUT
 fileName      Name of file to generate (Default = print to command
               window)

 Markus Herrgard

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function printSampleStats(samples,commonModel,sampleNames,fileName)
0002 %printSampleStats Print out sample statistics for multiple samples
0003 %
0004 % printSampleStats(samples,commonModel,sampleNames,fileName)
0005 %
0006 %INPUTS
0007 % samples       Samples to plot
0008 % commonModel   COBRA model structure
0009 % sampleNames   Names of samples
0010 %
0011 %OPTIONAL INPUT
0012 % fileName      Name of file to generate (Default = print to command
0013 %               window)
0014 %
0015 % Markus Herrgard
0016 
0017 if (nargin > 3)
0018   fid = fopen(fileName,'w');
0019 else
0020   fid = 1;
0021 end
0022 
0023 [sampleModes,sampleMeans,sampleStds,sampleMedians] = calcSampleStats(samples);
0024 
0025 fprintf(fid,'Rxn\t');
0026 if (isfield(commonModel,'subSystems'))
0027   fprintf(fid,'Subsystem\t');
0028 end
0029 
0030 for i = 1:length(sampleNames)
0031     fprintf(fid,'%s-mode\t',sampleNames{i});
0032 end
0033 for i = 1:length(sampleNames)
0034     fprintf(fid,'%s-mean\t',sampleNames{i});
0035 end
0036 for i = 1:length(sampleNames)
0037     fprintf(fid,'%s-median\t',sampleNames{i});
0038 end
0039 for i = 1:length(sampleNames)
0040     fprintf(fid,'%s-std\t',sampleNames{i});
0041 end
0042 fprintf(fid,'\n');
0043 
0044 for i = 1:length(commonModel.rxns)
0045     fprintf(fid,'%s\t',commonModel.rxns{i});
0046     if (isfield(commonModel,'subSystems'))
0047       fprintf(fid,'%s\t',commonModel.subSystems{i});
0048     end
0049     %for j = 1:length(samples)
0050         fprintf(fid,'%8.6f\t',sampleModes(i,:));
0051         %end
0052    %for j = 1:length(samples)
0053         fprintf(fid,'%8.6f\t',sampleMeans(i,:));
0054         %end
0055     %for j = 1:length(samples)
0056         fprintf(fid,'%8.6f\t',sampleMedians(i,:));
0057         %end
0058     %for j = 1:length(samples)
0059         fprintf(fid,'%8.6f\t',sampleStds(i,:));
0060         %end
0061     fprintf(fid,'\n');
0062 end
0063 if (fid > 1)
0064   fclose(fid);
0065 end

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