ttest.cpp

此程序为VC++常用数值算法这本书附赠的光盘中包含了本书中全部的源代码

void ttest(double data1[],int n1,double data2[],int n2,double& t,double& prob)
{
	double ave1,ave2,var1,var2;
    avevar(data1, n1, ave1, var1);
    avevar(data2, n2, ave2, var2);
    int df = n1 + n2 - 2;
    double var = ((n1 - 1) * var1 + (n2 - 1) * var2) / df;
    t = (ave1 - ave2) / sqrt(var * (1.0 / n1 + 1.0 / n2));
    prob = betai(0.5 * df, 0.5, df / (df + t * t));
}