plotHistConv

PURPOSE ^

plotHistConv Plot convergence of sample histograms

SYNOPSIS ^

function plotHistConv(model,samples,rxnNames,nSubSamples)

DESCRIPTION ^

plotHistConv   Plot convergence of sample histograms

  plotHistConv(model,samples,rxnNames,nSubSamples)

INPUTS
 model         COBRA model structure
 samples       Sample fluxes
 rxnNames      List of reactions to plot
 nSubSamples   Number of sub samples

 Markus Herrgard 8/14/06

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function plotHistConv(model,samples,rxnNames,nSubSamples)
0002 %plotHistConv   Plot convergence of sample histograms
0003 %
0004 %  plotHistConv(model,samples,rxnNames,nSubSamples)
0005 %
0006 %INPUTS
0007 % model         COBRA model structure
0008 % samples       Sample fluxes
0009 % rxnNames      List of reactions to plot
0010 % nSubSamples   Number of sub samples
0011 %
0012 % Markus Herrgard 8/14/06
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

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