mnl_seegview/plotCortexWeights.m at main · MultimodalNeuroimagingLab/mnl_seegview · 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
%
% Modified by HH 2021
%
function plotCortexWeights(locs, wts, threshFrac, wm)

    if nargin < 3 || isempty(threshFrac), threshFrac = 0; end % abs(values) below threshFrac*wm deemed not significant

    wts(isnan(wts)) = 0;

    % this function plots colored dots in brain rendering
    if nargin < 4
        wm=max(abs(wts)); % scale to max if no max given
    end
    wts(wts > wm) = wm; wts(wts < -wm) = -wm; % clip weights
    pthresh = threshFrac*wm;

    hold on
    for k=1:size(locs,1)

        if any(isnan(locs(k, :))), continue; end

        if abs(wts(k)) < pthresh % insignificant
            plot3(locs(k,1),locs(k,2),locs(k,3),'.',...
            'MarkerSize',16,...
            'Color',.01*[1 1 1])

            plot3(locs(k,1),locs(k,2),locs(k,3),'.',...
            'MarkerSize',14,...
            'Color',.35*[1 1 1])

        elseif wts(k) >= pthresh % positively significant

            plot3(locs(k,1),locs(k,2),locs(k,3),'.',...
            'MarkerSize',30*abs(wts(k))/wm+23,...
            'Color',.01*[1 1 1])

            hold on, plot3(locs(k,1),locs(k,2),locs(k,3),'.',...
            'MarkerSize',30*abs(wts(k))/wm+20,...
            'Color',.99*[1 1-wts(k)/wm 1-wts(k)/wm])

        else % negatively significant
            plot3(locs(k,1),locs(k,2),locs(k,3),'.',...
            'MarkerSize',30*abs(wts(k))/wm+23,...
            'Color',.01*[1 1 1])

            plot3(locs(k,1),locs(k,2),locs(k,3),'.',...
            'MarkerSize',30*abs(wts(k))/wm+20,...
            'Color',.99*[1+wts(k)/wm 1+wts(k)/wm 1])
        end
    end
    hold off

end