www.gusucode.com > som-bp混合神经网络的matlab程序源码 > matlab_emulator/data_statis/statis.m
% statit the feature clc clear all base_path = 'E:\SOMBP\data_source\feature02\'; % minmax([I1';I2';I3';I4';I5';I6';I7';I8']) I1 = load([base_path,'apen.txt']); % 1 >> 2 1.2563 1.8150 I2 = load([base_path,'kc.txt']); % 0 >> 1 0.6299 0.9611 I3 = load([base_path,'mir.txt']); % 0 >> 5 0.4802 4.0544 I4 = load([base_path,'asm.txt']); % 0 >> 0.01 0.0004 0.0037 I5 = load([base_path,'idm.txt']); % 0 >> 0.5 0.0663 0.2906 I6 = load([base_path,'cont.txt']); % 0 >> 500 12.2600 337.5440 I7 = load([base_path,'ent.txt']); % 5 >> 9 5.8675 7.9858 I8 = load([base_path,'gln.txt']); % 10 >> 60 13.4204 53.2048 sample_num = size(I1,1); negative_num = 75; positive_num = 55; feature_num = 8; figure; P = [(I1'-ones(1,sample_num ));I2';(I3'./5);(I4'./0.01);I5'./0.5;(I6'./500);(I7'-ones(1,sample_num ).*5)./4;(I8'-ones(1,sample_num ).*10)./50]; T = ones(1,sample_num); for(i=1:8) T = T.*P(i,:); end plot(T(1:positive_num),'r'); hold on plot(T(1+positive_num:sample_num),'b'); T = ([I1';I2';I3';I4';I5';I6';I7';I8']); FeatureName = ['ape','kc1','mir','asm','idm','cont','ent','gln']; figure; for(i=1:feature_num) subplot(2,4,i) plot(T(i,1:positive_num),'r'); hold on; plot(T(i,1+positive_num:sample_num),'b'); title(FeatureName((i-1)*3+1:i*3)); [h(i),significance(i),ci] = ttest2(T(i,1:positive_num),T(i,1+positive_num:sample_num),0.025) end