Welch's t-test: Difference between revisions

Line 42:
double calculate_Pvalue (const double *array1, const size_t array1_size, const double *array2, const size_t array2_size) {
if (array1_size <= 1) {
return 1.0;//not enough information to calculate p-value
}
if (array2_size <= 1) {
return 1.0;//not enough information to calculate p-value
}
double mean1 = 0.0, mean2 = 0.0;
Line 64:
}
if ((mean1 == mean2) && (variance1 == variance2)) {
return 1.0;//not worth calculating
}
if ((variance1 == 0.0) && (variance2 == 0.0)) {
return 1.0;
}
variance1 = variance1/(array1_size-1);
Line 87 ⟶ 90:
if ((a + 0.5) > 1755) {
denominator = 1.0;
for (unsigned int n = 1; n < 170100300000; n++) {
denominator *= (1 / (1+ a*0.5/(n*(a+0.5+n))));
}