ReciPSIICOS/plot_simulation_histogram.m at master · bioelectric-interfaces/ReciPSIICOS · GitHub

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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
% Input:
% r -- array with all biases of size (number_of_methods x 2 x
%   number_of_snrs x 1 x number_of_simulations), where:
%   number_of_methods = 4 (ReciPSIICOS, wReciPSIICOS, LCMV, MNE),
%   2 corresponds to synchronous/asynchronous
%   among all snrs used you need to choose the desired index as noise_i
% var -- array with all spreading values, the same atructure as r

% BIG FIGURE
c = lines(8);
noise_i = 4;
mc = size(r, 5);

% SYNCHRONOUS SOURCES
figure
subplot(4,1,1)
histogram(r(1,1,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(1,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(1,1,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(1,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)

subplot(4,1,2)
histogram(r(2,1,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(6,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(2,1,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(6,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)

subplot(4,1,3)
histogram(r(3,1,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(2,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(3,1,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(2,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)

subplot(4,1,4)
histogram(r(4,1,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(4,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(4,1,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(4,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)


% ASYNCHRONOUS SOURCES

figure
subplot(4,1,1)
histogram(r(1,2,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(1,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(1,2,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(1,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)

subplot(4,1,2)
histogram(r(2,2,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(6,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(2,2,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(6,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)

subplot(4,1,3)
histogram(r(3,2,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(2,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(3,2,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(2,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)

subplot(4,1,4)
histogram(r(4,1,noise_i,1,:), 'BinWidth', 0.002, 'FaceColor', c(4,:), ...
    'EdgeColor', 'k', 'FaceAlpha', 0.3, 'EdgeAlpha', 0.3, 'LineWidth', 2)
hold on
histogram(var(4,1,noise_i,1,:), 'BinWidth', 0.001, 'FaceColor', c(4,:), ...
    'FaceAlpha', 0.8, 'LineWidth', 2)
xlim([0, 0.08])
ylim([0 mc/4])
xlabel('Meters, m')
ylabel('Number of simulations')
set(gca,'FontSize', 18)
legend({'Estimation bias', 'Point spreading'}, 'FontSize', 24)