-
Notifications
You must be signed in to change notification settings - Fork 92
/
drosDoBootstrap.m
65 lines (54 loc) · 2.44 KB
/
drosDoBootstrap.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
53
54
55
56
57
58
59
60
61
62
63
64
65
% DROSDOBOOTSTRAP Perform bootstrap sampling to assess significance of differences in rankings
% FORMAT
% DESC Perform bootstrap sampling to assess significance of differences in rankings
%
% COPYRIGHT : Antti Honkela, 2009
% SHEFFIELDML
demRunRankings
%drosLoadData
%load results/rankings
thresholds = [5000];
REPEATS = 100000;
tt = [20, 100, 250];
matrices = {};
data = {};
tfs = {'twi', 'mef2'};
for i=1:2,
tf = tfs{i};
rankings = {indrank.(tf), disimrank.(tf), oderank.(tf), corrrank.(tf), mutarank.(tf)};
[a, b] = drosBootstrapEvaluation(drosexp, drosinsitu, rankings, chip_validation.(tf), REPEATS, tt, 0);
a, mean(b)
matrices = [matrices, {a}];
data = [data, {b}];
[a, b] = drosBootstrapEvaluation(drosexp, drosinsitu, rankings, mutant_validation.(tf), REPEATS, tt, 0);
a, mean(b)
matrices = [matrices, {a}];
data = [data, {b}];
[a, b] = drosBootstrapEvaluation(drosexp, drosinsitu, rankings, chip_validation.(tf), REPEATS, tt, 1);
a, mean(b)
matrices = [matrices, {a}];
data = [data, {b}];
[a, b] = drosBootstrapEvaluation(drosexp, drosinsitu, rankings, mutant_validation.(tf), REPEATS, tt, 1);
a, mean(b)
matrices = [matrices, {a}];
data = [data, {b}];
end
%i = 2; [a, b] = drosBootstrapEvaluation(drosexp2, drosinsitu, {indids{i}, disimids{i}, corrids{i}, mutaids{i}}, validation.(tfnames{i}), 1000, [20 100 250], 1); a, mean(b)
%i = 2; [a, b] = drosBootstrapEvaluation(drosexp2, drosinsitu, {indids{i}, disimids{i}, corrids{i}, mutaids{i}}, validation.(tfnames{i}), 1000, [20 100 250], 0); a, mean(b)
%save /share/work/ahonkela/bootstrap_results3.mat matrices data
fprintf('\\multicolumn{18}{c}{Twist global ChIP-chip} \\\\\n')
drosPrintBootstrapMatrices(matrices{1}, 1:5)
fprintf('\\multicolumn{18}{c}{Twist global knock-outs} \\\\\n')
drosPrintBootstrapMatrices(matrices{2}, 1:5, 5)
fprintf('\\multicolumn{18}{c}{Twist focused ChIP-chip} \\\\\n')
drosPrintBootstrapMatrices(matrices{3}, 1:5)
fprintf('\\multicolumn{18}{c}{Twist focused knock-outs} \\\\\n')
drosPrintBootstrapMatrices(matrices{4}, 1:5, 5)
fprintf('\\multicolumn{18}{c}{Mef2 global ChIP-chip} \\\\\n')
drosPrintBootstrapMatrices(matrices{5}, 1:5)
fprintf('\\multicolumn{18}{c}{Mef2 global knock-outs} \\\\\n')
drosPrintBootstrapMatrices(matrices{6}, 1:5, 5)
fprintf('\\multicolumn{18}{c}{Mef2 focused ChIP-chip} \\\\\n')
drosPrintBootstrapMatrices(matrices{7}, 1:5)
fprintf('\\multicolumn{18}{c}{Mef2 focused knock-outs} \\\\\n')
drosPrintBootstrapMatrices(matrices{8}, 1:5, 5)