sigImagesc

UNTITLED4 Summary of this function goes here
   shows data2 - data1 values
   only shows values that are significantly different (p < 0.05)
   data1 = 3D mat
   data2 = 3D mat. First 2 dim needs to be same size as those of data1
   labels = cell array of labels
   cLim = color bar limit

0001 function ax = sigImagesc(ax,data1,data2,labels,cLim)
0002 %UNTITLED4 Summary of this function goes here
0003 %   shows data2 - data1 values
0004 %   only shows values that are significantly different (p < 0.05)
0005 %   data1 = 3D mat
0006 %   data2 = 3D mat. First 2 dim needs to be same size as those of data1
0007 %   labels = cell array of labels
0008 %   cLim = color bar limit
0009 
0010 
0011 
0012 % find significant differences in connectivity matrix
0013 sigCon = false(size(data1,1),size(data1,2));
0014 for y = 1:size(data1,1)
0015     for x = 1:size(data1,2)
0016         group1 = data1(y,x,:);
0017         group1(isinf(group1)) = [];
0018         group1(isnan(group1)) = [];
0019         group2 = data2(y,x,:);
0020         group2(isinf(group2)) = [];
0021         group2(isnan(group2)) = [];
0022         if ~isempty(group1) || ~isempty(group2)
0023             h = ttest2(group1,group2);
0024             if ~isnan(h)
0025                 sigCon(y,x) = h;
0026             end
0027         end
0028     end
0029 end
0030 
0031 
0032 customCMap = blueWhiteRed(100);
0033 
0034 connectivityDiff = nanmean(data2,3) - nanmean(data1,3);
0035 connectivityDiff(~sigCon) = nan;
0036 imAlpha=ones(size(connectivityDiff));
0037 imAlpha(isnan(connectivityDiff))=0;
0038 imagesc(connectivityDiff,'AlphaData',imAlpha,cLim);
0039 set(gca,'color',0.5*[1 1 1]);
0040 xtickangle(60);
0041 colormap(customCMap);
0042 s4Pos = get(ax,'position');
0043 colorbar;
0044 set(ax,'Position',s4Pos);
0045 title('Connectivity diff (pre-post)');
0046 xtickLoc = 0.5:(size(data1,1)-0.5);
0047 set(gca,'XTick',xtickLoc);
0048 set(gca,'YTick',(1:size(data1,1)));
0049 set(gca,'XTickLabel',labels)
0050 set(gca,'YTickLabel',labels)
0051 set(gca,'TickLength',[0 0])
0052 set(gca, 'FontSize', 7);
0053 xlim([0.5 (size(data1,1)+0.5)]);
0054 ylim([0.5 (size(data1,1)+0.5)]);
0055 end
0056