ttest.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 11 行

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));
}