0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015 close all;
0016 clear all;
0017 goIncludeToolbox;
0018
0019
0020 parameter.MAX_LEVEL = 8;
0021 parameter.REPEAT_TIMES = 100;
0022 parameter.STRIDE_SEED_NUMBER = 10000;
0023 parameter.REWIRE_PROB_RING = 0.3;
0024 parameter.ANALYSIS_RECURRENTION = false;
0025
0026
0027 REWIRE = [0.01 0.03 0.05 0.1 0.15 0.20 0.4 0.5];
0028
0029
0030 loadStruct(parameter, 'base');
0031
0032
0033 ceNetwork = bioNetwork('./CE_ExpData/C_elegans_Matrix.mat', 'mat');
0034 ceNetwork = readLabel(ceNetwork, './CE_ExpData/C_elegans_Label.txt');
0035 ceNetwork = setIO(ceNetwork, ceNetwork.inputList, ceNetwork.outputList);
0036 rand('seed', 0);
0037 ceNetwork = levelAnalysis(ceNetwork, MAX_LEVEL);
0038
0039
0040
0041 [meanNetworks, allNetworks, stdNetworks] = ...
0042 makeNetAndAnalysis(ceNetwork, parameter);
0043
0044 loadStruct(meanNetworks, 'base');
0045 loadStruct(allNetworks, 'base');
0046 loadStruct(stdNetworks, 'base');
0047 clear meanNetworks allNetworks stdNetworks
0048
0049
0050 cxNetwork = bioNetwork('./PcbExpData/CXGraph_AllRevised.mat', 'mat');
0051 cxNetwork = readLabel(cxNetwork, './PcbExpData/CXLabel_AllRevised.txt');
0052 cxNetwork = setIO(cxNetwork, cxNetwork.inputList, cxNetwork.outputList);
0053 rand('seed', 0);
0054 cxNetwork = levelAnalysis(cxNetwork, MAX_LEVEL);
0055
0056
0057 tmpMeanNumEdge = full(mean(sum(ceNetwork.matrix')));
0058
0059 for iRewire = 1:numel(REWIRE)
0060 smallWireNet{iRewire} = cell(1, REPEAT_TIMES);
0061 for jNet = 1:REPEAT_TIMES
0062 seedNumber = jNet * STRIDE_SEED_NUMBER;
0063 smallWireNet{iRewire}{jNet} = smallWorldNetwork(ceNetwork.size, ...
0064 REWIRE(iRewire), seedNumber, tmpMeanNumEdge);
0065 smallWireNet{iRewire}{jNet} = setIO(smallWireNet{iRewire}{jNet}, ...
0066 ceNetwork.inputNumber, ceNetwork.outputNumber);
0067 smallWireNet{iRewire}{jNet} = ...
0068 levelAnalysis(smallWireNet{iRewire}{jNet}, MAX_LEVEL);
0069 end
0070
0071 meanSmallNet{iRewire} = meanFieldInCells(smallWireNet{iRewire});
0072 stdSmallNet{iRewire} = stdFieldInCells(smallWireNet{iRewire});
0073
0074 smallWireNet{iRewire} = [];
0075 end
0076 clear smallWireNet
0077
0078
0079 for iNet = 1:REPEAT_TIMES
0080 seedNumber = iNet * STRIDE_SEED_NUMBER;
0081 rand('seed', seedNumber);
0082 randIoCeNet{iNet} = ...
0083 bioNetwork('./CE_ExpData/C_elegans_Matrix.mat', 'mat');
0084 randIoCeNet{iNet} = setIO(randIoCeNet{iNet}, ...
0085 ceNetwork.inputNumber, ceNetwork.outputNumber);
0086 randIoCeNet{iNet} = levelAnalysis(randIoCeNet{iNet}, MAX_LEVEL);
0087 end
0088 meanRandIoCeNet = meanFieldInCells(randIoCeNet);
0089 stdRandIoCeNet = stdFieldInCells(randIoCeNet);
0090
0091
0092 for iNet = 1:REPEAT_TIMES
0093 seedNumber = iNet * STRIDE_SEED_NUMBER;
0094 rand('seed', seedNumber);
0095 sepIoCeNet{iNet} = ...
0096 bioNetwork('./CE_ExpData/C_elegans_Matrix.mat', 'mat');
0097 sepIoCeNet{iNet} = ...
0098 readLabel(sepIoCeNet{iNet}, './CE_ExpData/C_elegans_Label.txt');
0099
0100 sepIoCeNet{iNet} = setIoByModule(sepIoCeNet{iNet}, ...
0101 numel(sepIoCeNet{iNet}.inputList), numel(sepIoCeNet{iNet}.outputList));
0102 sepIoCeNet{iNet} = levelAnalysis(sepIoCeNet{iNet}, MAX_LEVEL);
0103 end
0104 meanSepIoCeNet = meanFieldInCells(sepIoCeNet);
0105 stdSepIoCeNet = stdFieldInCells(sepIoCeNet);
0106
0107
0108 revIoCeNet = bioNetwork('./CE_ExpData/C_elegans_Matrix.mat', 'mat');
0109 revIoCeNet = readLabel(revIoCeNet, './CE_ExpData/C_elegans_Label.txt');
0110 revIoCeNet = ...
0111 setIO(revIoCeNet, revIoCeNet.outputList, revIoCeNet.inputList);
0112 rand('seed', 0);
0113 revIoCeNet = levelAnalysis(revIoCeNet, MAX_LEVEL);