leaddbs/ea_corrplot_gen.m at master · ebakstein/leaddbs · 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
function [R_upd,p_upd,R,p,f]=ea_corrplot(varargin)
% this simple function is a small wrapper for a corrplot figure.


X=varargin{1};
description=varargin{2};
labels=varargin{3};
if nargin>3
    handles=varargin{4};
end
if nargin>4
    if ~isempty(varargin{5})
        if ischar(varargin{5}) || isequal([1 3],size(varargin{5}))
            color=varargin{5};
        else
            groups=varargin{5};
        end
    else
        groups=ones(length(X),1);
    end
else
    groups=ones(length(X),1);
end

corrtype='Pearson';
try corrtype=varargin{6}; end

dim=size(X,2);

disp(description);


nnans=isnan(sum(X,2));
X(nnans,:)=[];
try groups(nnans)=[]; end

switch corrtype
    case {'permutation_spearman','permutation'}
        [R,p]=ea_permcorr(X(:,1),X(:,2),'spearman');
        R_upd=R; p_upd=p;
    case 'permutation_pearson'
        [R,p]=ea_permcorr(X(:,1),X(:,2),'pearson');
        R_upd=R; p_upd=p;
    otherwise
[R,p]=corr(X,'rows','pairwise','type',corrtype);
R_upd=R(2:end,1);
p_upd=p(2:end,1);
end


% support for nans


%labels=M.stats(1).ea_stats.atlases.names(lidx);


jetlist=lines;
%jetlist(groups,:);
for area=1:length(R_upd)
    %% plot areas:
    f=figure('color','w','name',description);
    if exist('color','var')
        scatter(X(:,1),X(:,area+1),[],'o','MarkerEdgeColor','w','MarkerFaceColor',color);
    else
        try
        scatter(X(:,1),X(:,area+1),[],'o','MarkerEdgeColor','w','MarkerFaceColor',jetlist(groups,:));
        catch
        scatter(X(:,1),X(:,area+1),[],jetlist(groups,:),'filled');

        end

    end

    h=lsline;
    set(h,'color','k');
    %axis square

    ax=gca;
% ax.YTick=[-20:20:90];
% ax.YLim=[-15,90];
% ax.XTick=[-20:10:90];
% ax.XLim=[25,65];
[~,fn]=fileparts(labels{area+1});
if length(fn)>4
    if strcmp(fn(end-3:end),'.nii')
        [~,fn]=fileparts(fn);
    end
end
    title([description,' (R=',sprintf('%.3f',R_upd(area)),', p=',sprintf('%.3f',p_upd(area)),').'],'FontSize',16,'FontName','Helvetica');
    xlabel(sub2space(labels{1}),'FontSize',16,'FontName','Helvetica');
    ylabel(labels{2},'FontSize',16,'FontName','Helvetica');
 %   spacing=mean([nanvar(X(:,1)),nanvar(X(:,area+1))]);
 %   xlim([nanmin(X(:,1))-spacing,nanmax(X(:,2))+spacing]);
 %   ylim([nanmin(X(:,area+1))-spacing,nanmax(X(:,area+1))+spacing]);
    if nargin>3
        if ~isempty(varargin{4})
        odir=get(handles.groupdir_choosebox,'String');
        ofname=[odir,description,'_',fn,'_',labels{1},'.png'];
        ea_screenshot(ofname);
        end
    end
end


function str=sub2space(str) % replaces subscores with spaces
str(str=='_')=' ';