commonmcode/metaPlot_2grp.m at master · blab-lab/commonmcode · 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
function metaPlot_2grp(meas, fld1, fld2, measName, y_label, colors, varargin)
%% Input options and configuratoins
fontSize = 12;
if ~isempty(fsic(varargin, '--font-size'))
    fontSize = varargin{fsic(varargin, '--font-size') + 1};
end

if ~isempty(fsic(varargin, 'fontSize'))
    fontSize = varargin{fsic(varargin, 'fontSize') + 1};
end

bPaired = ~isempty(fsic(varargin, '--paired'));

if ~isempty(fsic(varargin, '--fld-labels'))
    lbl_fld1 = varargin{fsic(varargin, '--fld-labels') + 1};
    lbl_fld2 = varargin{fsic(varargin, '--fld-labels') + 2};
else
    lbl_fld1 = fld1;
    lbl_fld2 = fld2;
end

%%
wndW = 500;
wndH = 360;
if ~isempty(fsic(varargin, 'wndW'))
    wndW = varargin{fsic(varargin, 'wndW') + 1};
end
if ~isempty(fsic(varargin, 'wndH'))
    wndH = varargin{fsic(varargin, 'wndH') + 1};
end

figure('name', measName, 'Position', [200, 200, wndW, wndH]);
set(gca, 'FontSize', fontSize)
if ~isempty(fsic(varargin, 'zeroLine'))
    plot([0, 3], [0, 0], '-', 'Color', [0.5, 0.5, 0.5]);
    hold on;
end

boxplot([meas.(fld1), meas.(fld2)], ...
        [1 * ones(size(meas.(fld1))), 2 * ones(size(meas.(fld2)))]);
hold on;
% boxplot([meas.(fld1), NaN NaN], [1 * ones(size(meas.(fld1))), 2, 2], 'Color', colors.(fld1));
% hold on;
% boxplot([NaN, NaN, meas.(fld1)], [1, 1, 2 * ones(size(meas.(fld1)))], 'Color', colors.(fld2));

if isempty(fsic(varargin, 'noIndPts'))
    plot(repmat(0.8, 1, numel(meas.(fld1))), meas.(fld1), 'o', 'Color', colors.(fld1));
    plot(repmat(1.8, 1, numel(meas.(fld2))), meas.(fld2), 'o', 'Color', colors.(fld2));
end

XLim = [0, 3];
if ~isempty(fsic(varargin, 'wndH'))
    XLim = varargin{fsic(varargin, 'XLim') + 1};
end

set(gca, 'XLim', XLim, 'XTick', [1, 2], 'XTickLabel', {lbl_fld1, lbl_fld2});

plot(1.2, mean(meas.(fld1)), 'o', 'Color', colors.(fld1));
plot([1.2, 1.2], mean(meas.(fld1)) + [-1, 1] * ste(meas.(fld1)), '-', 'Color', colors.(fld1));
plot(2.2, mean(meas.(fld2)), 'o', 'Color', colors.(fld2));
plot([2.2, 2.2], mean(meas.(fld2)) + [-1, 1] * ste(meas.(fld2)), '-', 'Color', colors.(fld2));

cohend = abs(cohen_d(meas.(fld2), meas.(fld1)));
xs = get(gca, 'XLim'); ys = get(gca, 'YLim');


if bPaired
    text(xs(1) + 0.5 * range(xs), ys(2) - 0.05 * range(ys), ...
         sprintf('N=%d', numel(meas.(fld1))), 'FontSize', fontSize);
    [h_t, p_t] = ttest(meas.(fld2), meas.(fld1));
    text(xs(1) + 0.5 * range(xs), ys(2) - 0.10 * range(ys), sprintf('Paired t-test: p=%.3e', p_t), 'FontSize', fontSize);
else
    text(xs(1) + 0.5 * range(xs), ys(2) - 0.05 * range(ys), ...
         sprintf('%s: N=%d; %s: N=%d', fld1, numel(meas.(fld1)), fld2, numel(meas.(fld2))), 'FontSize', fontSize);
    [h_t, p_t] = ttest2(meas.(fld2), meas.(fld1));
    text(xs(1) + 0.5 * range(xs), ys(2) - 0.10 * range(ys), sprintf('t-test: p=%.5e', p_t), 'FontSize', fontSize);
end

text(xs(1) + 0.5 * range(xs), ys(2) - 0.15 * range(ys), sprintf('Cohen''s d=%.5f', cohend), 'FontSize', fontSize);
[p_rs] = ranksum(meas.(fld2), meas.(fld1));
% text(xs(1) + 0.035 * range(xs), ys(2) - 0.15 * range(ys), sprintf('Rank-sum: p = %.5f', p_rs));

ylabel(y_label);

for i0 = 1 : 2
    if i0 == 1
        fld = fld1;
    else
        fld = fld2;
    end
    fprintf('%s: N = %d, mean = %f, SD = %f, SE = %f\n', fld, numel(meas.(fld)), mean(meas.(fld)), std(meas.(fld)), ste(meas.(fld)));
end
return