original_ranksum.m

来自「秩和检验的Matlab源码程序」· M 代码 · 共 22 行

%==========================================================================
%
% 函数名：original_ranksum
% 函数介绍：求初始排列的秩和。=检验统计量W。
% 输入参数：a是一个基因的归一化表达量。
%          m是正常样本的个数。
%          column是总的样本个数。
% 输出参数：r是初始秩和，就是算法中的检验统计量。
%
%==========================================================================
function [ r ] = original_ranksum ( a, column, m )
B = sort ( a );
r = 0;
for i = 1 : column
    for j = 1 : m
        if B ( i ) == a ( j )
            r = r + i;
        end
    end
end
%end