0001 function plotHistConv(model,samples,rxnNames,nSubSamples)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 nSkip = 10;
0015 nBin = 20;
0016
0017 [nRxns,nSamples] = size(samples);
0018
0019 [isInModel,rxnInd] = ismember(rxnNames,model.rxns);
0020 rxnInd = rxnInd(isInModel);
0021
0022 nPlotRxn = sum(isInModel);
0023 rxnNames = rxnNames(isInModel);
0024
0025 nCol = ceil(sqrt(nPlotRxn));
0026 nRow = ceil(nPlotRxn/nCol);
0027
0028 subSampleSize = floor(nSamples/nSubSamples);
0029
0030 clf
0031 for rxnID = 1:nPlotRxn
0032 subplot(nRow,nCol,rxnID);
0033 hold on
0034 maxx = -1e9;
0035 minx = 1e9;
0036 maxy = 0;
0037 for subID = 1:nSubSamples
0038 [n,x] = hist(samples(rxnInd(rxnID),1:nSkip:(subSampleSize*subID))',nBin);
0039 plot(x,n/sum(n));
0040 maxx = max([max(x) maxx]);
0041 minx = min([min(x) minx]);
0042 maxy = max([max(n/sum(n)) maxy]);
0043 end
0044 axis([minx maxx 0 maxy]);
0045 title(rxnNames{rxnID});
0046 end