testing_statistics.m

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

%==========================================================================
%
% 函数名：testing_statistics
% 函数介绍：求t检验的检验统计量t。
% 输入参数：a是基因i的某种表达。可能是原始表达，可能是扰动后的表达。
%          r1是normal样本的个数。
%          r2是tumor样本的个数。
% 输出参数：t是t检验的检验统计量。
%==========================================================================
function [ t ] = testing_statistics ( a, r1, r2 )


%将两类样本分开。
x = a ( 1 : r1 );%normal
y = a (  r1 + 1 : r1 + r2 );%tumor

r = r1 + r2 - 2;

%求样本均值和方差。
    m1  = mean ( x );%样本1的均值
    m2  = mean ( y );%样本2的均值
    
    s1  = (var ( x ))*(r1 - 1) / r;%样本1的方差。var（）通过除以n-1达到标准化，n为样本个数，所以这里要乘以n-1。
    s2  = (var ( y ))*(r2 - 1) / r;%样本2的方差。

%求检验统计量

t = ( m2 - m1 ) / sqrt ( (( r1 * s1 ) + ( r2 * s2 )) * (1 / r1 + 1 / r2) / r );

%end of function