cii.m

来自「分类信息指数是基于信噪比的计算变量所含分类信息的算法」· M 代码 · 共 27 行

%==========================================================================

% 函数名：cii
% 函数介绍： 求基因的分类信息指数。
% 输入参数： x是基因表达谱数据的标准化表达量。
%           r1是第一类样本的个数。
% 输出参数：information_index是基因的分类信息指数。
%          row是基因表达谱数据中基因的个数。

%==========================================================================

function [information_index, row] = cii ( x, r1 )
[row, column] = size ( x );    %row是基因个数；column是总样本个数
information_index = zeros ( 1, row );
for i = 1 : row  %gene循环
    a1 = x( i, 1 : r1);   %基因i，第一组样本的表达量。
    a2 = x( i, r1 + 1 : column ); %基因i，第二组样本的表达量。
    %样本均值
    u1 = mean ( a1 );
    u2 = mean ( a2 );
    %样本标准差
    s1 = std ( a1 );  %标准差（=方差开平方）
    s2 = std ( a2 );
    information_index ( i ) = abs( u1 - u2 ) / (2 * (s1 + s2)) + (1/2) * log2((s1^2 + s2^2) / (2 * s1 * s2));
end
% end