goPlotHorz

PURPOSE ^

This script produces the figure 5b. It plot the horizontal propagation versus

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 This script produces the figure 5b. It plot the horizontal propagation versus
 propagation levels in C. Elegans neural network, small-world network, random
 network, regular network and Drosophila's CX network. The curves are based on
 fitting Gompertz function. The errorbar is standard deviation.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % This script produces the figure 5b. It plot the horizontal propagation versus
0002 % propagation levels in C. Elegans neural network, small-world network, random
0003 % network, regular network and Drosophila's CX network. The curves are based on
0004 % fitting Gompertz function. 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 FIG_SIZE = [0 0 1200 600];
0013 lineWidth = 3;
0014 markerSize = 10;
0015 markerScale = 1.5;
0016 
0017 % fitting with Gompertz function
0018 cfit.CE = fitGompertz(0:(MAX_LEVEL - 1), ceNetwork.HorzProp);
0019 cfit.Rand = fitGompertz(0:(MAX_LEVEL - 1), meanRandNet.HorzProp);
0020 cfit.Regular = fitGompertz(0:(MAX_LEVEL - 1), meanRegularNet.HorzProp);
0021 cfit.Ring = fitGompertz(0:(MAX_LEVEL - 1), meanRingNet.HorzProp);
0022 cfit.CX = fitGompertz(0:(MAX_LEVEL - 1), cxNetwork.HorzProp);
0023 
0024 % save cfit object to each network
0025 ceNetwork.cfitHorzProp = cfit.CE;
0026 meanRandNet.cfitHorzProp = cfit.Rand;
0027 meanRegularNet.cfitHorzProp = cfit.Regular;
0028 meanRingNet.cfitHorzProp = cfit.Ring;
0029 cxNetwork.cfitHorzProp = cfit.CX;
0030 
0031 % plot marker and errorbar of each networks
0032 figure('name', 'Horizontal Propagation', 'outerPosition', FIG_SIZE, 'color', 'w');
0033 hold on;
0034 tmpP = linspace(0, MAX_LEVEL - 1, 100);
0035 plot(0:(MAX_LEVEL - 1), ceNetwork.HorzProp, 'Marker', netMarker('CE'), ...
0036     'Color', netColor('CE'), 'MarkerSize', markerSize, ...
0037     'MarkerFaceColor', netMarkerFace('CE'), ...
0038     'LineWidth', lineWidth, 'LineStyle', 'none');
0039 errorbar(0:(MAX_LEVEL - 1), meanRandNet.HorzProp, ...
0040     stdRandNet.HorzProp, 'Marker', netMarker('ER'), ...
0041     'Color', netColor('ER'), 'MarkerSize', markerSize * markerScale, ...
0042     'MarkerFaceColor', netMarkerFace('ER'), ...
0043     'LineWidth', lineWidth, 'LineStyle', 'none');
0044 errorbar(0:(MAX_LEVEL - 1), meanRegularNet.HorzProp, ...
0045     stdRegularNet.HorzProp, 'Marker', netMarker('RL'), ...
0046     'Color', netColor('RL'), 'MarkerSize', markerSize, ...
0047     'MarkerFaceColor', netMarkerFace('RL'), ...
0048     'LineWidth', lineWidth, 'LineStyle', 'none');
0049 errorbar(0:(MAX_LEVEL - 1), meanRingNet.HorzProp, ...
0050     stdRingNet.HorzProp, 'Marker', netMarker('SW'), ...
0051     'Color', netColor('SW'), 'MarkerSize', markerSize, ...
0052     'MarkerFaceColor', netMarkerFace('SW'), ...
0053     'LineWidth', lineWidth, 'LineStyle', 'none');
0054 plot(0:(MAX_LEVEL - 1), cxNetwork.HorzProp, 'Marker', netMarker('CX'), ...
0055     'Color', netColor('CX'), 'MarkerSize', markerSize * markerScale, ...
0056     'MarkerFaceColor', netMarkerFace('CX'), ...
0057     'LineWidth', lineWidth, 'LineStyle', 'none');
0058 
0059 % plot fitting curve of each network
0060 plot(tmpP, cfit.CE(tmpP), 'Color', netColor('CE'), 'LineWidth', 4);
0061 plot(tmpP, cfit.Rand(tmpP), 'Color', netColor('ER'), 'LineWidth', 4);
0062 plot(tmpP, cfit.Regular(tmpP), 'Color', netColor('RL'), 'LineWidth', 4);
0063 plot(tmpP, cfit.Ring(tmpP), 'Color', netColor('SW'), 'LineWidth', 4);
0064 plot(tmpP, cfit.CX(tmpP), 'Color', netColor('CX'), 'LineWidth', 4);
0065 plot([2 2], [0 cxNetwork.cfitHorzProp(2)], '--k', 'lineWidth', 3);
0066 set(gca, 'FontSize', 18, 'lineWidth', 4, 'FontWeight', 'bold');
0067 
0068 xlim([0 4]);
0069 ylim([0 1.1]);
0070 xlabel('Propagation level');
0071 ylabel('Horizontal propagation');
0072 set(gca, 'xTick', [0 1 2 3 4], 'yTick', [0 0.4 0.8 1])
0073 box off;

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