standardization.m

计算扰动后的T检验Matlab程序

%==========================================================================
%
% 函数名：standardization
% 函数介绍：对原始数据进行标准化。方法：（x－基因最小值）/（基因最大值－基因最小值）
%          计算每类样本的个数。
% 输入参数：x是原始基因表达谱数据。行是基因，列是样本。x的第一行是样本标记。0表示normal；1表示tumor。
% 输出参数：y是原始数据的标准化。
%          style_1是normal样本的个数。
%          style_2是tumor样本的个数。
%          gene_number是基因表达谱数据中的基因个数。
%==========================================================================

function [ y, style_1, style_2, gene_number  ] = standardization ( x )

style = x ( 1 , : );%找到样本标记行。
[ m, n ] = size ( x );
gene_number = m-1;%基因表达谱数据中共有基因m-1个。

%对基因表达谱数据标准化。
x = x ( 2 : m , : );%将样本标记行去掉。

for i = 1 : m - 1
    B = max (x (i, :));
    S = min (x (i, :));
    for j = 1 : n
        y(i, j) = (x (i, j) - S) / (B - S);   
    end
end 
%标准化结束

%计算每类样本的个数。
style_1 = 0;%normal的个数。
style_2 = 0;%tumor的个数。

for i = 1 : n
    if style (i) == 0
        style_1 = style_1 + 1;
    else
        style_2 = style_2 + 1;
    end
end
%计算样本个数结束。

%end of function