Welch's t-test: Difference between revisions

Content added Content deleted

Inline

@@ Line 69: / Line 69: @@
 	variance2 = variance2/(array2_size-1);
 	const double WELCH_T_STATISTIC = (mean1-mean2)/sqrt(variance1/array1_size+variance2/array2_size);
-	const double DEGREES_OF_FREEDOM = pow((variance1/array1_size+variance2/array2_size),2.0)//numerator
+	const double DEGREES_OF_FREEDOM = pow((variance1/array1_size+variance2/array2_size),2.0)
 	 /
 	(
@@ Line 84: / Line 84: @@
            sum2 += (pow(h * i,a-1))/(sqrt(1-h * i));
 	}
-	return ((h / 6.0) * ((pow(x,a-1))/(sqrt(1-x)) + 4.0 * sum1 + 2.0 * sum2))/(tgammal(a)*1.77245385090551610/tgammal(a+0.5));//heavily simplified to reduce function calls and speed computation
+	return ((h / 6.0) * ((pow(x,a-1))/(sqrt(1-x)) + 4.0 * sum1 + 2.0 * sum2))/(tgammal(a)*1.77245385090551610/tgammal(a+0.5));
 }
 //-------------------