goPlotSummaryIO

PURPOSE ^

This script produces the figure 7c. It plot the vertical propagation versus

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 This script produces the figure 7c. It plot the vertical propagation versus
 horizontal propagation of C. Elegans neural network, small-world network,
 random IO CE network, separated IO CE network, reversed IO CE network and
 distant IO SW network in level 2. The errorbar is standard deviation.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % This script produces the figure 7c. It plot the vertical propagation versus
0002 % horizontal propagation of C. Elegans neural network, small-world network,
0003 % random IO CE network, separated IO CE network, reversed IO CE network and
0004 % distant IO SW network in level 2. The errorbar is standard deviation.
0005 
0006 %   ---------
0007 %   Yen-Nan Lin, NTHU, 2010-2014, Matlab 2012a
0008 
0009 close all;
0010 
0011 % plot setting
0012 xErrBarWidth = 0.015;
0013 yErrBarWidth = 0.015;
0014 textSize = 16;
0015 lineWidth = 3;
0016 markerSize = 10;
0017 markerScaleL = 1.5;
0018 markerScaleS = 1.2;
0019 
0020 FIG_SIZE = [0 0 600 650];
0021 figure('name', 'summary', 'outerPosition', FIG_SIZE, 'color', 'w');
0022 
0023 % plot marker of CE, rand IO CE, separated IO CE, reversed IO CE, distant IO SW
0024 % and SW networks
0025 hold on;
0026 plot(ceNetwork.HorzProp(3), ceNetwork.VertProp(3), ...
0027     'Marker', netMarker('CE'), 'Color', netColor('CE'), ...
0028     'MarkerFaceColor', netMarkerFace('CE'), ...
0029     'MarkerSize', markerSize, 'LineWidth', lineWidth);
0030 plot(meanRandIoCeNet.HorzProp(3), meanRandIoCeNet.VertProp(3), ...
0031     'Marker', netMarker('RandIoCE'), 'Color', netColor('RandIoCE'), ...
0032     'MarkerFaceColor', netMarkerFace('RandIoCE'), ...
0033     'MarkerSize', markerSize * markerScaleS, 'LineWidth', lineWidth);
0034 plot(meanSepIoCeNet.HorzProp(3), meanSepIoCeNet.VertProp(3), ...
0035     'Marker', netMarker('SepIoCE'), 'Color', netColor('SepIoCE'), ...
0036     'MarkerFaceColor', netMarkerFace('SepIoCE'), ...
0037     'MarkerSize', markerSize * markerScaleL, 'LineWidth', lineWidth);
0038 plot(revIoCeNet.HorzProp(3), revIoCeNet.VertProp(3), ...
0039     'Marker', netMarker('RevIoCE'), 'Color', netColor('RevIoCE'), ...
0040     'MarkerFaceColor', netMarkerFace('RevIoCE'), ...
0041     'MarkerSize', markerSize * markerScaleS, 'LineWidth', lineWidth);
0042 plot(meanRingNet.HorzProp(3), meanRingNet.VertProp(3), ...
0043     'Marker', netMarker('SW'), 'Color', netColor('SW'), ...
0044     'MarkerFaceColor', netMarkerFace('SW'), ...
0045     'MarkerSize', markerSize, 'LineWidth', lineWidth);
0046 plot(meanRingDistNet.HorzProp(3), meanRingDistNet.VertProp(3), ...
0047     'Marker', netMarker('DistSW'), 'Color', netColor('DistSW'), ...
0048     'MarkerFaceColor', netMarkerFace('DistSW'), ...
0049     'MarkerSize', markerSize, 'LineWidth', lineWidth);
0050 set(gca, 'FontSize', 18, 'lineWidth', 4, 'FontWeight', 'bold');
0051 xlabel('Horizontal propagation');
0052 ylabel('Vertical propagation');
0053 xlim([0, 1]);
0054 ylim([0.4, 1]);
0055 box off;
0056 
0057 % plot errorbar of rand IO CE, separated IO CE, SW, distant IO SW networks
0058 hold on;
0059 tmpHorz = cellfun(@(x) x.HorzProp(3), randIoCeNet);
0060 tmpVert = cellfun(@(x) x.VertProp(3), randIoCeNet);
0061 hErr = ploterr(mean(tmpHorz), mean(tmpVert), std(tmpHorz), std(tmpVert), ...
0062     'abshhx', xErrBarWidth, 'abshhy', yErrBarWidth);
0063 set(hErr, 'lineWidth', 2, 'Color', netColor('RandIoCE'));
0064 
0065 tmpHorz = cellfun(@(x) x.HorzProp(3), sepIoCeNet);
0066 tmpVert = cellfun(@(x) x.VertProp(3), sepIoCeNet);
0067 hErr = ploterr(mean(tmpHorz), mean(tmpVert), std(tmpHorz), std(tmpVert), ...
0068     'abshhx', xErrBarWidth, 'abshhy', yErrBarWidth);
0069 set(hErr, 'lineWidth', 2, 'Color', netColor('SepIoCE'));
0070 
0071 tmpHorz = cellfun(@(x) x.HorzProp(3), smallRingNet);
0072 tmpVert = cellfun(@(x) x.VertProp(3), smallRingNet);
0073 hErr = ploterr(mean(tmpHorz), mean(tmpVert), std(tmpHorz), std(tmpVert), ...
0074     'abshhx', xErrBarWidth, 'abshhy', yErrBarWidth);
0075 set(hErr, 'lineWidth', 2, 'Color', netColor('SW'));
0076 
0077 tmpHorz = cellfun(@(x) x.HorzProp(3), smallRingDistNet);
0078 tmpVert = cellfun(@(x) x.VertProp(3), smallRingDistNet);
0079 hErr = ploterr(mean(tmpHorz), mean(tmpVert), std(tmpHorz), std(tmpVert), ...
0080     'abshhx', xErrBarWidth, 'abshhy', yErrBarWidth);
0081 set(hErr, 'lineWidth', 2, 'Color', netColor('DistSW'));
0082 hold off;
0083 
0084 % plot the dash line of small-world networks with rewiring probability from
0085 % 0.01 to 0.5 (defined in goAnalysis.m 'REWIRE')
0086 hold on;
0087 tmpVert = [meanRegularNet.VertProp(3), ...
0088     cellfun(@(x) x.VertProp(3), meanSmallNet(1:6)), ...
0089     meanRingNet.VertProp(3), ...
0090     cellfun(@(x) x.VertProp(3), meanSmallNet(7:end)), ...
0091     meanRandNet.VertProp(3)];
0092 tmpHorz = [meanRegularNet.HorzProp(3), ...
0093     cellfun(@(x) x.HorzProp(3), meanSmallNet(1:6)), ...
0094     meanRingNet.HorzProp(3), ...
0095     cellfun(@(x) x.HorzProp(3), meanSmallNet(7:end)), ...
0096     meanRandNet.HorzProp(3)];
0097 plot(tmpHorz, tmpVert, ':', 'lineWidth', 3, 'Color', netColor('SW'));
0098 hold off;
0099 
0100 % annotation
0101 text(double(ceNetwork.HorzProp(3)), double(ceNetwork.VertProp(3)), ...
0102     '   CE', 'FontSize', textSize, 'FontWeight', 'bold');
0103 text(double(meanRingNet.HorzProp(3)), double(meanRingNet.VertProp(3)), ...
0104     '   SW', 'FontSize', textSize, 'FontWeight', 'bold');
0105 text(double(meanRandIoCeNet.HorzProp(3)), double(meanRandIoCeNet.VertProp(3)), ...
0106     {'   Random', '   I/O CE'}, 'FontSize', textSize, 'FontWeight', 'bold', ...
0107     'VerticalAlignment', 'bottom');
0108 text(double(meanSepIoCeNet.HorzProp(3)), double(meanSepIoCeNet.VertProp(3)), ...
0109     {'Separated  ', 'I/O CE  '}, 'FontSize', textSize, 'FontWeight', 'bold', ...
0110     'VerticalAlignment', 'bottom', 'HorizontalAlignment', 'right');
0111 text(double(revIoCeNet.HorzProp(3)), double(revIoCeNet.VertProp(3)), ...
0112     {'  Reversed', '  I/O CE'}, 'FontSize', textSize, 'FontWeight', 'bold', ...
0113     'VerticalAlignment', 'top');
0114 text(double(meanRingDistNet.HorzProp(3)), double(meanRingDistNet.VertProp(3) - 0.02), ...
0115     {'Distant', 'I/O SW'}, 'FontSize', textSize, 'FontWeight', 'bold', ...
0116     'VerticalAlignment', 'top', 'HorizontalAlignment', 'center');
0117 
0118 breakYaxis(2);
0119 set(gca, 'xtick', 0:0.2:1);

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