Welch's t-test: Difference between revisions
Content added Content deleted
No edit summary |
No edit summary |
||
Line 59: | Line 59: | ||
for (size_t x = 0; x < array2_size; x++) { |
for (size_t x = 0; x < array2_size; x++) { |
||
mean2 += array2[x]; |
mean2 += array2[x]; |
||
⚫ | |||
⚫ | |||
⚫ | |||
} |
} |
||
mean1 /= array1_size; |
mean1 /= array1_size; |
||
Line 68: | Line 71: | ||
for (size_t x = 0; x < array2_size; x++) { |
for (size_t x = 0; x < array2_size; x++) { |
||
variance2 += (array2[x]-mean2)*(array2[x]-mean2); |
variance2 += (array2[x]-mean2)*(array2[x]-mean2); |
||
⚫ | |||
⚫ | |||
⚫ | |||
} |
} |
||
if ((variance1 == 0.0) && (variance2 == 0.0)) { |
if ((variance1 == 0.0) && (variance2 == 0.0)) { |
||
Line 93: | Line 93: | ||
sum2 += (pow(h * i,a-1))/(sqrt(1-h * i)); |
sum2 += (pow(h * i,a-1))/(sqrt(1-h * i)); |
||
} |
} |
||
⚫ | |||
double denominator; |
|||
if ((a + 0.5) > 1755) { |
|||
denominator = 1.0; |
|||
for (unsigned int n = 1; n < 300000; n++) { |
|||
denominator /= (1+ a*0.5/(n*(a+0.5+n))); |
|||
} |
|||
denominator = (a+0.5)*denominator/(a*0.5); |
|||
} else { |
|||
denominator = tgammal(a)*1.77245385090551610/tgammal(a+0.5); |
|||
} |
|||
⚫ | |||
} |
} |
||
//------------------- |
//------------------- |