calcSampleDifference

PURPOSE ^

calcSampleDifference Calculate a flux difference sample between two conditions

SYNOPSIS ^

function [sampleDiff,sampleRatio] = calcSampleDifference(sample1,sample2,nPts)

DESCRIPTION ^

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

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