0001 function printSampleStats(samples,commonModel,sampleNames,fileName)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
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
0050 fprintf(fid,'%8.6f\t',sampleModes(i,:));
0051
0052
0053 fprintf(fid,'%8.6f\t',sampleMeans(i,:));
0054
0055
0056 fprintf(fid,'%8.6f\t',sampleMedians(i,:));
0057
0058
0059 fprintf(fid,'%8.6f\t',sampleStds(i,:));
0060
0061 fprintf(fid,'\n');
0062 end
0063 if (fid > 1)
0064 fclose(fid);
0065 end