Upload New File

Data Analysis/analyzeInPhaseMat.m

+function [dataFitsInd, dataFitsPooled, dataFused] = analyzeInPhaseMat(pathIn)
+
+if ~exist('pathIn', 'var')
+    % user interaction to get folder that should be worked on
+    pathIn = uigetdir;
+end
+filei_recursive = dir([pathIn '/**/ResultsNew.mat']);
+nFiles = length(filei_recursive);
+
+% list of control variables
+listControlVar = {'labelingFraction', 'saltConcentration', 'pH', 'drugConcentration'};
+
+% reset counters
+counterFused = 0;
+counterFitsInd = 1;
+for k = 1 : nFiles
+    tic
+    
+    % path to the current results file
+    currFilePath = fullfile(filei_recursive(k).folder, filei_recursive(k).name);
+    
+    % diplay current folder and filename
+    disp(['Folder ', num2str(k), ' of ', num2str(nFiles)])
+    disp(['Current Path: ', currFilePath])
+    
+    % load the new results file
+    clear('Results')
+    load(currFilePath, 'Results');
+    
+    % check for data with no detected emulsion droplets and remove
+    fieldNames1 = fieldnames(Results);
+    for kk = 1 : numel(fieldNames1)
+        fieldNames2 = fieldnames(Results.(fieldNames1{kk}));
+        for kkk = 1 : length(fieldNames2)
+            indx(kk,kkk) = isempty(Results.(fieldNames1{kk}).(fieldNames2{kkk}));
+        end
+    end
+    
+    indx = sum(indx') > 0;
+    Results = rmfield(Results, fieldNames1(indx));
+    
+    
+    % prepare output of all data files fused into one struct
+    fieldi = fieldnames(Results);
+    for kk = 1 : length(fieldi)
+        currName = ['Experiment_',  num2str(kk + counterFused, '%04d')];
+        dataFused.(currName) = Results.(fieldi{kk});
+    end
+    counterFused = counterFused + length(fieldi);
+    
+    % scan through list of control variables
+    for kk = 1 : length(listControlVar)
+        controlVars{kk} = unique(cell2mat(structCrawler(Results, listControlVar{kk})'));
+        nAssays(kk) = length(controlVars{kk});
+    end
+    indx = nAssays > 1;
+    
+    % individual temperatures are fitted
+    % if any other control variable was changed in one experiment it will
+    % be also separated and fitted individually for each temperature
+    if  sum(nAssays) == length(listControlVar)
+        temperatures = unique(cell2mat(structCrawler(Results, 'temperature')'));
+        for kk = 1 : length(temperatures)
+            currResults = sortPDExperiments(Results, 'temperature', temperatures(kk));
+            dataFitsInd(counterFitsInd) = fitLinearConcSeries_V03_wBW(currResults);
+            counterFitsInd = counterFitsInd + 1;
+        end
+    else
+        for kk = 1 : nAssays(indx)
+            currVar = controlVars{indx};
+            currResults = sortPDExperiments(Results, listControlVar{indx}, currVar(kk));
+            temperatures = unique(cell2mat(structCrawler(currResults, 'temperature')'));
+            for kkk = 1 : length(temperatures)
+                currResultsT = sortPDExperiments(currResults, 'temperature', temperatures(kkk));
+                dataFitsInd(counterFitsInd) = fitLinearConcSeries_V03_wBW(currResultsT);
+                counterFitsInd = counterFitsInd + 1;
+            end
+        end
+        
+    end
+    toc
+    
+end
+
+% aggregate all datasets of same control variables and fit them pooled
+for k = 1 : length(listControlVar)
+    controlVars{k} = unique(cell2mat(structCrawler(dataFused, listControlVar{k})'));
+    nAssays(k) = length(controlVars{k});
+end
+
+counterFitsFused = 1;
+currVar1 = unique(cell2mat(structCrawler(dataFused, listControlVar{1})'));
+for k = 1 : length(currVar1)
+    currResults1 = sortPDExperiments(dataFused, listControlVar{1}, currVar1(k));
+    currVar2 = unique(cell2mat(structCrawler(currResults1, listControlVar{2})'));
+    for kk = 1 : length(currVar2)
+        currResults2 = sortPDExperiments(currResults1, listControlVar{2}, currVar2(kk));
+        currVar3 = unique(cell2mat(structCrawler(currResults2, listControlVar{3})'));
+        for kkk = 1 : length(currVar3)
+            currResults3 = sortPDExperiments(currResults2, listControlVar{3}, currVar3(kkk));
+            temperatures = unique(cell2mat(structCrawler(currResults3, 'temperature')'));
+            for kkkk = 1 : length(temperatures)
+                currResults4 = sortPDExperiments(currResults3, 'temperature', temperatures(kkkk));
+                dataFitsPooled(counterFitsFused) = fitLinearConcSeries_V03_wBW(currResults4);
+                counterFitsFused = counterFitsFused + 1;
+            end
+        end
+    end
+end
+
+end
\ No newline at end of file