disturb_t.m

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

%==========================================================================
%
% 函数名：disturb_t
% 函数介绍：求扰动下的检验统计量，记做小写的t。
% 输入参数：data是基因i的原始表达。
%          r1是正常样本的个数。
%          r2是肿瘤样本的个数。
%          
% 输出参数：disturb_t是扰动实验中得到的所有的t。t是一个向量，向量大小等于扰动次数。
%          
% 说明：该函数设定扰动次数（实验次数）disturb_time = (r1+r2)!/(r1!*r2!)。该数值可根据需要自行设定，但是不能太小。
%==========================================================================
function [ disturb_t ] = disturb_t ( data, r1, r2 )

disturb_time = ( factorial(r1 + r2) ) / (  factorial(r1) * factorial(r2)  ) - 1;%扰动次数，即实验次数。减1是因为原始组合的t已经求过了。

disturb_combination ( 1: r1  ) = 1 : r1 ;%第一次扰动的组合情况。

for i = 1 : disturb_time 
    [ disturb_combination ] = combination ( disturb_combination, r1 + r2, r1 );    %生成{1,2,3...,n}的r—组合
   
    %扰动之后的基因i的表达
    disturb_data  = data ( disturb_combination  ) ;
    k = r1;%disturb_data的指针。
    nature = 1 : r1 + r2;%自然组合。
    
    for n = 1 : r1 + r2
        F = 0;
        for m = 1 : r1
            if  nature ( n ) == disturb_combination ( m ) 
                F = 1;
                break;
            end
        end
        
        if F == 0
            k = k + 1;
            disturb_data ( k ) = data ( nature ( n ) ) ;
        end
    end
    
    %求检验统计量t。
    [ t ] = testing_statistics ( disturb_data, r1, r2);
    disturb_t ( i ) = t;
    
end
%end