statis.m
上传用户:sdcdgc2008
上传日期:2014-04-13
资源大小:365k
文件大小:2k
开发平台:

Matlab

  1. % statit the feature
  2. clc
  3. clear all
  4. base_path = 'E:SOMBPdata_sourcefeature02';                      % minmax([I1';I2';I3';I4';I5';I6';I7';I8'])
  5. I1 = load([base_path,'apen.txt']);           % 1 >> 2        1.2563    1.8150        
  6. I2 = load([base_path,'kc.txt']);             % 0 >> 1        0.6299    0.9611
  7. I3 = load([base_path,'mir.txt']);            % 0 >> 5        0.4802    4.0544
  8. I4 = load([base_path,'asm.txt']);            % 0 >> 0.01    0.0004    0.0037
  9. I5 = load([base_path,'idm.txt']);            % 0 >> 0.5       0.0663    0.2906
  10. I6 = load([base_path,'cont.txt']);           % 0 >> 500     12.2600  337.5440
  11. I7 = load([base_path,'ent.txt']);            % 5 >> 9        5.8675    7.9858
  12. I8 = load([base_path,'gln.txt']);            % 10 >> 60      13.4204   53.2048
  13. sample_num = size(I1,1);
  14. negative_num = 75;
  15. positive_num = 55;
  16. feature_num = 8;
  17. figure;
  18. 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];
  19. T = ones(1,sample_num);
  20. for(i=1:8)
  21.     T = T.*P(i,:);
  22. end
  23. plot(T(1:positive_num),'r');
  24. hold on
  25. plot(T(1+positive_num:sample_num),'b');
  26. T = ([I1';I2';I3';I4';I5';I6';I7';I8']);
  27. FeatureName = ['ape','kc1','mir','asm','idm','cont','ent','gln'];
  28. figure;
  29. for(i=1:feature_num)
  30.     subplot(2,4,i)
  31.     plot(T(i,1:positive_num),'r');
  32.     hold on;
  33.     plot(T(i,1+positive_num:sample_num),'b');
  34.     title(FeatureName((i-1)*3+1:i*3));
  35.     [h(i),significance(i),ci] = ttest2(T(i,1:positive_num),T(i,1+positive_num:sample_num),0.025)
  36. end