choose.m

分类信息指数是基于信噪比的计算变量所含分类信息的算法

%==========================================================================
%
% 函数名：choose
% 函数介绍：选出分类分类信息指数大于阈值的基因（备选特征基因），输出他们的归一化表达量和基因的顺序号。
%          运行该函数时需要设定一个合适的阈值。阈值可以是柱图的突跳点，也可以根据经验自行设定或其它方法求得。
% 输入参数：c是所有基因的分类信息指数向量。
%          m是原始基因表达谱数据中总的基因个数。
%          data_stand是原始数据的归一化表达量。
% 输出参数：data_important是备选特征基因的原始表达量。
%          order是备选特征基因的顺序号。
%          c_important是备选特征基因的分类信息指数。
%          std_dataimportant是备选特征基因的归一化表达量。

%===================================================
function [ data_important, order, c_important, std_dataimportant ] = choose ( c, m, data_stand,data )

gene_data = data ( 2 : m + 1 , : ); %基因的原始表达量。
cii_threshold = input ('请输入基因分类信息指数阈值：');

p = 1;
for i = 1 : m
    if c ( i ) > cii_threshold            
        std_dataimportant ( p, : ) = data_stand ( i, : );
        data_important ( p , : ) = gene_data ( i , : );
        order ( p ) = i;
        c_important ( p ) = c ( i );
        p = p + 1;
    end
end