goPlotChConnectCE

PURPOSE ^

This script produces the figure 11. It plot channel connectivity matrix of

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 This script produces the figure 11. It plot channel connectivity matrix of
 C. Elegans network with highlighted modules in propagation levels 0-3.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % This script produces the figure 11. It plot channel connectivity matrix of
0002 % C. Elegans network with highlighted modules in propagation levels 0-3.
0003 
0004 %   ---------
0005 %   Yen-Nan Lin, NTHU, 2010-2014, Matlab 2012a
0006 
0007 close all;
0008 
0009 NUM_FIG_COL = 4;
0010 FIG_SIZE = [0, 0, 1500, 300];
0011 MAX_PLOT_LEVEL = 4;
0012 REDUCE_COLORBAR_WIDTH = 0.01;
0013 
0014 % setting of figure
0015 figure('outerPosition', FIG_SIZE, 'color', 'w');
0016 panelObj = panel();
0017 panelObj.pack(1, NUM_FIG_COL);
0018 panelObj.fontsize = 15;
0019 panelObj.fontname = 'Arial';
0020 panelObj.fontweight = 'bold';
0021 panelObj.margin = [10 15 3 8];
0022 panelObj.de.margin = 3;
0023 panelObj.title('CE Network');
0024 panelObj.ylabel('Input nodes');
0025 panelObj.xlabel('Output nodes');
0026 
0027 % plot channel connectivity matrix of C. Elegans
0028 nInput = ceNetwork.inputNumber;
0029 nOutput = ceNetwork.outputNumber;
0030 for iLevel = 1:MAX_PLOT_LEVEL
0031     panelObj(1, iLevel).select();
0032     plotOpt.showXTick = true;
0033     if iLevel == 1
0034         plotOpt.showYTick = true;
0035     else
0036         plotOpt.showYTick = false;
0037     end
0038     myImagesc(ceNetwork.ChConnectMat{iLevel}, plotOpt);
0039     panelObj(1, iLevel).title(['Level ', int2str(iLevel - 1)]);
0040 
0041     % highlight modules
0042     for jModule = 1:ceNetwork.NumModule
0043         tmpInputIxList = find(ceNetwork.moduleId(1:nInput) == jModule);
0044         tmpOutputIxList = find(...
0045             ceNetwork.moduleId((nInput + 1):(nInput + nOutput)) == jModule);
0046         if isempty(tmpInputIxList) || isempty(tmpOutputIxList)
0047             continue;
0048         end
0049         xPosition = min(tmpOutputIxList) - 0.5;
0050         yPosition = min(tmpInputIxList) - 0.5;
0051         width = numel(tmpOutputIxList);
0052         height = numel(tmpInputIxList);
0053         rectangle('position', [xPosition, yPosition, width, height], ...
0054             'edgeColor', 'w', 'lineWidth', 4);
0055     end
0056 end

Generated on Thu 30-Jan-2014 00:00:07 by m2html © 2005