Welch's t-test: Difference between revisions

Content added Content deleted
Line 76: Line 76:
);
);
const double a = DEGREES_OF_FREEDOM/2, x = DEGREES_OF_FREEDOM/(WELCH_T_STATISTIC*WELCH_T_STATISTIC+DEGREES_OF_FREEDOM);
const double a = DEGREES_OF_FREEDOM/2, x = DEGREES_OF_FREEDOM/(WELCH_T_STATISTIC*WELCH_T_STATISTIC+DEGREES_OF_FREEDOM);
const double N = 3000, h = x/N;//only about 5000 is needed for 15 digits precision
const double N = 3000, h = x/N;
double sum1 = 0.0, sum2 = 0.0;
double sum1 = 0.0, sum2 = 0.0;
for(unsigned int i = 0;i < N; i++) {
for(unsigned int i = 0;i < N; i++) {