0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017 close all;
0018 clear all;
0019 goIncludeToolbox;
0020
0021
0022 parameter.MAX_LEVEL = 8;
0023 parameter.REPEAT_TIMES = 100;
0024 parameter.STRIDE_SEED_NUMBER = 10000;
0025 parameter.REWIRE_PROB_RING = 0.3;
0026 parameter.ANALYZE_RECURRENT = true;
0027
0028
0029 loadStruct(parameter, 'base');
0030
0031
0032 ceNetwork = bioNetwork('./CE_ExpData/C_elegans_Matrix.mat', 'mat');
0033 ceNetwork = readLabel(ceNetwork, './CE_ExpData/C_elegans_Label.txt');
0034 ceNetwork = setIO(ceNetwork, ceNetwork.inputList, ceNetwork.outputList);
0035 rand('seed', 0);
0036 ceNetwork = levelAnalysis(ceNetwork, MAX_LEVEL, ANALYZE_RECURRENT);
0037
0038
0039
0040 ceNetwork.NumConnect = full(sum(sum(ceNetwork.matrix)));
0041 ceNetwork.meanNumEdge = full(mean(sum(ceNetwork.matrix')));
0042
0043
0044
0045
0046 randProb = ceNetwork.NumConnect / (ceNetwork.size^2 - ceNetwork.size);
0047 randNetwork = cell(1, REPEAT_TIMES);
0048 for iTime = 1:REPEAT_TIMES
0049 seedNumber = iTime * STRIDE_SEED_NUMBER;
0050 randNetwork{iTime} = randomNetwork(ceNetwork.size, randProb, seedNumber);
0051 randNetwork{iTime} = setIO(randNetwork{iTime}, ...
0052 ceNetwork.inputNumber, ceNetwork.outputNumber);
0053 randNetwork{iTime} = ...
0054 levelAnalysis(randNetwork{iTime}, MAX_LEVEL, ANALYZE_RECURRENT);
0055 end
0056
0057 meanRandNet = meanFieldInCells(randNetwork);
0058 stdRandNet = stdFieldInCells(randNetwork);
0059
0060
0061 smallRingNet = cell(1, REPEAT_TIMES);
0062 for iTime = 1:REPEAT_TIMES
0063 seedNumber = iTime * STRIDE_SEED_NUMBER;
0064 smallRingNet{iTime} = smallWorldNetwork(ceNetwork.size, ...
0065 REWIRE_PROB_RING, seedNumber, ceNetwork.meanNumEdge);
0066 smallRingNet{iTime} = setIO(smallRingNet{iTime}, ...
0067 ceNetwork.inputNumber, ceNetwork.outputNumber);
0068 smallRingNet{iTime} = ...
0069 levelAnalysis(smallRingNet{iTime}, MAX_LEVEL, ANALYZE_RECURRENT);
0070 end
0071
0072 meanRingNet = meanFieldInCells(smallRingNet);
0073 stdRingNet = stdFieldInCells(smallRingNet);
0074
0075
0076 cxNetwork = bioNetwork('./PcbExpData/CXGraph_AllRevised.mat', 'mat');
0077 cxNetwork = readLabel(cxNetwork, './PcbExpData/CXLabel_AllRevised.txt');
0078 cxNetwork = setIO(cxNetwork, cxNetwork.inputList, cxNetwork.outputList);
0079 rand('seed', 0);
0080 cxNetwork = levelAnalysis(cxNetwork, MAX_LEVEL, ANALYZE_RECURRENT);