Skip to content
Snippets Groups Projects
Commit fcbe0ff4 authored by jiglesia's avatar jiglesia
Browse files

Update fitLinearConcSeries_V03_wBW.m

parent 5ad69c8f
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Data Analysis/fitLinearConcSeries_V03_wBW.m
+ 7
− 7
View file @ fcbe0ff4
...@@ -144,12 +144,12 @@ end ...@@ -144,12 +144,12 @@ end
indxSort = proteinConcentrations >= min(unqProteinC(~dataExclude)); indxSort = proteinConcentrations >= min(unqProteinC(~dataExclude));
% determine cout via linear fit to volume fraction via fluorescence % determine c_out via linear fit to volume fraction via fluorescence
try try
if moreThanOne if moreThanOne
[FitVinVtot, GoF] = fit(proteinConcentrations(indxSort), volumeFraction(indxSort), linFit_Eq ,... [FitVinVtot, GoF] = fit(proteinConcentrations(indxSort), volumeFraction(indxSort), linFit_Eq ,...
'StartPoint', [0.01, -0.01], 'Robust', 'Bisquare'); % 'Lower', [0, -1], 'Upper', [1, 0] 'StartPoint', [0.01, -0.01], 'Robust', 'Bisquare'); % 'Lower', [0, -1], 'Upper', [1, 0]
% cout via formula form mass conservation % c_out via formula form mass conservation
currA = FitVinVtot.a; currA = FitVinVtot.a;
currB = FitVinVtot.b; currB = FitVinVtot.b;
currConf = confint(FitVinVtot, 0.95); currConf = confint(FitVinVtot, 0.95);
...@@ -158,11 +158,11 @@ try ...@@ -158,11 +158,11 @@ try
currCout = - currB / currA; currCout = - currB / currA;
currCout_err = sqrt( (-1/currA)^2 * currB_err^2 + ... currCout_err = sqrt( (-1/currA)^2 * currB_err^2 + ...
(currB/(currA^2))^2 * currA_err^2); (currB/(currA^2))^2 * currA_err^2);
% cin via formula form mass conservation % c_in via formula form mass conservation
currCin = (1 - currB ) / currA; currCin = (1 - currB ) / currA;
currCin_err = sqrt( (-1/currA)^2 * currB_err^2 + ... currCin_err = sqrt( (-1/currA)^2 * currB_err^2 + ...
((currB-1)/(currA^2))^2 * currA_err^2); ((currB-1)/(currA^2))^2 * currA_err^2);
% this would mean a negative cout % this would mean a negative c_out
if currB > 0 if currB > 0
currA = NaN; currA = NaN;
currA_err = NaN; currA_err = NaN;
...@@ -209,7 +209,7 @@ try ...@@ -209,7 +209,7 @@ try
'Lower', [0, -1], 'Upper', [1, 0], 'StartPoint', [0.01, -0.01],... 'Lower', [0, -1], 'Upper', [1, 0], 'StartPoint', [0.01, -0.01],...
'Robust', 'Bisquare', 'Weights', weightsBF); 'Robust', 'Bisquare', 'Weights', weightsBF);
% cout via formula form mass conservation % c_out via formula form mass conservation
currABF = FitVinVtot_BF.a; currABF = FitVinVtot_BF.a;
currBBF = FitVinVtot_BF.b; currBBF = FitVinVtot_BF.b;
currConf = confint(FitVinVtot_BF, 0.66); currConf = confint(FitVinVtot_BF, 0.66);
...@@ -218,7 +218,7 @@ try ...@@ -218,7 +218,7 @@ try
currCoutBF = - currBBF / currABF; currCoutBF = - currBBF / currABF;
currCoutBF_err = sqrt( (-1/currABF)^2 * currBBF_err^2 + ... currCoutBF_err = sqrt( (-1/currABF)^2 * currBBF_err^2 + ...
(currBBF/(currABF^2))^2 * currABF_err^2); (currBBF/(currABF^2))^2 * currABF_err^2);
% cin via formula form mass conservation % c_in via formula form mass conservation
currCinBF = (1 - currBBF ) / currABF; currCinBF = (1 - currBBF ) / currABF;
currCinBF_err = sqrt( (-1/currABF)^2 * currBBF_err^2 + ... currCinBF_err = sqrt( (-1/currABF)^2 * currBBF_err^2 + ...
((currBBF-1)/(currABF^2))^2 * currABF_err^2); ((currBBF-1)/(currABF^2))^2 * currABF_err^2);
...@@ -332,4 +332,4 @@ catch ...@@ -332,4 +332,4 @@ catch
dataOut.experimentDate = nan; dataOut.experimentDate = nan;
end end
end end
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment