calcSampleDifference Calculate a flux difference sample between two conditions sampleDiff = calcSampleDifference(sample1,sample2,nPts) Selects randomly nPts flux vectors from sample1 and sample2 and calcutes the difference between the flux vectors INPUTS sample1 First flux sample sample2 Second flux sample nPts Number of flux difference profiles desired OUTPUTS sampleDiff sampleRatio Markus Herrgard 11/13/06
0001 function [sampleDiff,sampleRatio] = calcSampleDifference(sample1,sample2,nPts) 0002 %calcSampleDifference Calculate a flux difference sample between two conditions 0003 % 0004 % sampleDiff = calcSampleDifference(sample1,sample2,nPts) 0005 % 0006 % Selects randomly nPts flux vectors from sample1 and sample2 and calcutes 0007 % the difference between the flux vectors 0008 % 0009 %INPUTS 0010 % sample1 First flux sample 0011 % sample2 Second flux sample 0012 % nPts Number of flux difference profiles desired 0013 % 0014 %OUTPUTS 0015 % sampleDiff 0016 % sampleRatio 0017 % 0018 % Markus Herrgard 11/13/06 0019 0020 [nFlux1,nSample1] = size(sample1); 0021 [nFlux2,nSample2] = size(sample2); 0022 0023 if (nFlux1 ~= nFlux2) 0024 error('Samples have different numbers of rxns'); 0025 else 0026 select1 = randint(nPts,1,nSample1)+1; 0027 select2 = randint(nPts,1,nSample2)+1; 0028 sampleDiff = sample2(:,select2) - sample1(:,select1); 0029 sampleRatio = sample2(:,select2)./sample1(:,select1); 0030 end 0031