www.pudn.com > SOM_Bp_HybridNetwork_matlab_emulator.rar > 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