-
Notifications
You must be signed in to change notification settings - Fork 0
/
automatedEntropyScript.m
52 lines (40 loc) · 1.24 KB
/
automatedEntropyScript.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
%% automatic entropy from color and FSC/SSC data
[fileNames,path] = uigetfile('E:\zon_lab\FACS\*.fcs','Multiselect','On');
if iscell(fileNames) == 1
numFiles = numel(fileNames);
else
numFiles = 1;
fileNames = {fileNames};
end
entropy = cell(numFiles,1);
score = cell(numFiles,1);
stdEntropy = cell(numFiles,1);
giniCoeff = cell(numFiles,1);
stdGiniCoeff = cell(numFiles,1);
rho = cell(numFiles,1);
sampleName = cell(numFiles,1);
multiWaitbar('analyzing files...',0);
for kk = 1:numFiles
[entropy{kk}, stdEntropy{kk}, giniCoeff{kk}, stdGiniCoeff{kk}, rho{kk}, sampleName{kk}, score{kk}] = autoEntropy(fullfile(path,fileNames{kk}),3);
multiWaitbar('analyzing files...', kk./numFiles);
end
multiWaitbar('CloseAll');
gini = [giniCoeff{:}]';
recombination = normalize_var(recomb,0,1);
modGini = gini.*recombination;
figure, hGiniRecomb = scatter(recomb,gini,50,'filled');
title('Recombation vs. Gini Coeff');
figure, hGini = stem(1:numFiles,gini);
title('Gini Coefficient');
xticks(1:numFiles);
xticklabels(sampleName);
xtickangle(45);
xlim([0, (numFiles + 1)]);
ylim([0 1]);
figure, hModGini = stem(1:numFiles, modGini);
title('Modified Gini');
xticks(1:numFiles);
xticklabels(sampleName);
xtickangle(45);
xlim([0, (numFiles + 1)]);
ylim([0,1]);