plotmeangraph.m

function plotmeangraph(listname,bandidx,varargin)

loadpaths
param = finputcheck(varargin, {
    'group', 'string', [], 'crsdiag'; ...
    'groupnames', 'cell', {}, {'UWS','MCS-','MCS+','EMCS','LIS','CTRL'}; ...
    });

load(sprintf('%s/groupdata_%s.mat',filepath,listname));
load(sprintf('sortedlocs_%d.mat',size(allcoh,3)));

groupvar = subjlist.(param.group);
groups = unique(groupvar(~isnan(groupvar)));

bands = {
    'delta'
    'theta'
    'alpha'
    'beta'
    'gamma'
    };

plotqt = 0.7;

for g = 1:length(groups)
    fprintf('Averaging %d subjects in group %d.\n', sum(groupvar == groups(g)), groups(g));
    groupcoh(g,:,:) = squeeze(mean(allcoh(groupvar == groups(g),bandidx,:,:),1));
    threshcoh(g,:,:) = threshold_proportional(squeeze(groupcoh(g,:,:)),1-plotqt);
    for c = 1:size(threshcoh,2)
        groupdeg(g,c) = sum(threshcoh(g,c,:))/(size(threshcoh,2)-1);
    end
end

% erange = [min(nonzeros(threshcoh(:))) max(threshcoh(:))];
% vrange = [min(nonzeros(groupdeg(:))) max(groupdeg(:))];
erange = [0 0.4];
vrange = [0 0.4];

for g = size(groupcoh,1):-1:1
    minfo(g,:) = plotgraph3d(squeeze(groupcoh(g,:,:)),'plotqt',plotqt,'escale',erange,'vscale',vrange,'arcs','strength');
    
    camva(8);
    camtarget([-9.7975  -28.8277   41.8981]);
    campos([-1.7547    1.7161    1.4666]*1000);
    camzoom(1.25);
    fprintf('%s %s - number of modules: %d\n',param.groupnames{g},bands{bandidx},length(unique(minfo(g,:))));
    set(gcf,'Name',sprintf('%s %s',param.groupnames{g},bands{bandidx}));
    set(gcf,'InvertHardCopy','off');
    print(gcf,sprintf('%s/figures/meangraph_%s_%s.tif',filepath,param.groupnames{g},bands{bandidx}),'-dtiff','-r150');
    %     saveas(gcf,sprintf('figures/meangraph_%s_%s.fig',grouplist{g},bands{bandidx}));
    close(gcf);
end