243 |
|
|
244 |
/* Tests for F distribution */ |
/* Tests for F distribution */ |
245 |
|
|
246 |
|
/* p(x, nu1, nu2) := betaI(1 / (1 + nu2 / nu1 * x), nu1 / 2, nu2 / 2) */ |
247 |
|
|
248 |
TEST(gsl_cdf_Fdist_P, (0.0, 1.2, 1.3), 0.0, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (0.0, 1.2, 1.3), 0.0, TEST_TOL6); |
249 |
TEST(gsl_cdf_Fdist_P, (1e-100, 1.2, 1.3), 6.98194275525039002e-61, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (1e-100, 1.2, 1.3), 6.98194275525039002e-61, TEST_TOL6); |
250 |
TEST(gsl_cdf_Fdist_P, (0.001, 1.2, 1.3), 1.10608485860238564e-2, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (0.001, 1.2, 1.3), 1.10608485860238564e-2, TEST_TOL6); |
291 |
TEST(gsl_cdf_Fdist_P, (1000.0, 1.2, 500), 1.0, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (1000.0, 1.2, 500), 1.0, TEST_TOL6); |
292 |
TEST(gsl_cdf_Fdist_P, (10000.0, 1.2, 500), 1.0, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (10000.0, 1.2, 500), 1.0, TEST_TOL6); |
293 |
|
|
|
/* [(p(x, nu1, nu2) := betaI(1 / (1 + nu2 / nu1 * x), nu1 / 2, nu2 / 2))] */ |
|
294 |
TEST(gsl_cdf_Fdist_P, (0.0, 200, 500), 0.0, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (0.0, 200, 500), 0.0, TEST_TOL6); |
295 |
TEST(gsl_cdf_Fdist_P, (1e-100, 200, 500), 0.0, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (1e-100, 200, 500), 0.0, TEST_TOL6); |
296 |
TEST(gsl_cdf_Fdist_P, (0.001, 200, 500), 4.09325080403669893e-251, TEST_TOL6); |
TEST(gsl_cdf_Fdist_P, (0.001, 200, 500), 4.09325080403669893e-251, TEST_TOL6); |
308 |
|
|
309 |
/* Tests for gamma distribution */ |
/* Tests for gamma distribution */ |
310 |
|
|
311 |
|
/* p(x, a, b) := gammaP(b, x / a) */ |
312 |
|
|
313 |
TEST(gsl_cdf_gamma_P, (0.0, 1, 1), 0.0, TEST_TOL6); |
TEST(gsl_cdf_gamma_P, (0.0, 1, 1), 0.0, TEST_TOL6); |
314 |
TEST(gsl_cdf_gamma_P, (1e-100, 1, 1), 1e-100, TEST_TOL6); |
TEST(gsl_cdf_gamma_P, (1e-100, 1, 1), 1e-100, TEST_TOL6); |
315 |
TEST(gsl_cdf_gamma_P, (0.001, 1, 1), 9.99500166625008332e-4 , TEST_TOL6); |
TEST(gsl_cdf_gamma_P, (0.001, 1, 1), 9.99500166625008332e-4 , TEST_TOL6); |
325 |
TEST(gsl_cdf_gamma_P, (1000.0, 1, 1), 1e0, TEST_TOL6); |
TEST(gsl_cdf_gamma_P, (1000.0, 1, 1), 1e0, TEST_TOL6); |
326 |
TEST(gsl_cdf_gamma_P, (10000.0, 1, 1), 1e0, TEST_TOL6); |
TEST(gsl_cdf_gamma_P, (10000.0, 1, 1), 1e0, TEST_TOL6); |
327 |
|
|
328 |
|
TEST(gsl_cdf_gamma_P, (0.0, 10, 1), 0.0, TEST_TOL6); |
329 |
|
TEST(gsl_cdf_gamma_P, (1e-100, 10, 1), 1e-101, TEST_TOL6); |
330 |
|
TEST(gsl_cdf_gamma_P, (0.001, 10, 1), 9.99950001666625001e-5 , TEST_TOL6); |
331 |
|
TEST(gsl_cdf_gamma_P, (0.01, 10, 1), 9.99500166625008332e-4, TEST_TOL6); |
332 |
|
TEST(gsl_cdf_gamma_P, (0.1, 10, 1), 9.95016625083194643e-3, TEST_TOL6); |
333 |
|
TEST(gsl_cdf_gamma_P, (0.325, 10, 1), 3.19775501686939529e-2, TEST_TOL6); |
334 |
|
TEST(gsl_cdf_gamma_P, (1.0, 10, 1), 9.51625819640404268e-2, TEST_TOL6); |
335 |
|
TEST(gsl_cdf_gamma_P, (1.5, 10, 1), 1.39292023574942193e-1, TEST_TOL6); |
336 |
|
TEST(gsl_cdf_gamma_P, (2.0, 10, 1), 1.81269246922018141e-1, TEST_TOL6); |
337 |
|
TEST(gsl_cdf_gamma_P, (10.0, 10, 1), 6.32120558828557678e-1, TEST_TOL6); |
338 |
|
TEST(gsl_cdf_gamma_P, (20.0, 10, 1), 8.64664716763387308e-1, TEST_TOL6); |
339 |
|
TEST(gsl_cdf_gamma_P, (100.0, 10, 1), 9.99954600070237515e-1, TEST_TOL6); |
340 |
|
TEST(gsl_cdf_gamma_P, (1000.0, 10, 1), 1e0, TEST_TOL6); |
341 |
|
TEST(gsl_cdf_gamma_P, (10000.0, 10, 1), 1e0, TEST_TOL6); |
342 |
|
|
343 |
|
TEST(gsl_cdf_gamma_P, (0.0, 10, 17), 0e0, TEST_TOL6); |
344 |
|
TEST(gsl_cdf_gamma_P, (1e-100, 10, 17), 0e0, TEST_TOL6); |
345 |
|
TEST(gsl_cdf_gamma_P, (0.001, 10, 17), 2.81119174040422844e-83, TEST_TOL6); |
346 |
|
TEST(gsl_cdf_gamma_P, (0.01, 10, 17), 2.80880324651985887e-66, TEST_TOL6); |
347 |
|
TEST(gsl_cdf_gamma_P, (0.1, 10, 17), 2.78502998087492130e-49, TEST_TOL6); |
348 |
|
TEST(gsl_cdf_gamma_P, (0.325, 10, 17), 1.37283653245125844e-40, TEST_TOL6); |
349 |
|
TEST(gsl_cdf_gamma_P, (1.0, 10, 17), 2.55811932292544243e-32, TEST_TOL6); |
350 |
|
TEST(gsl_cdf_gamma_P, (1.5, 10, 17), 2.40420441175422372e-29, TEST_TOL6); |
351 |
|
TEST(gsl_cdf_gamma_P, (2.0, 10, 17), 3.05092926217898577e-27, TEST_TOL6); |
352 |
|
TEST(gsl_cdf_gamma_P, (10.0, 10, 17), 1.094920130378183e-15, TEST_TOL6); |
353 |
|
TEST(gsl_cdf_gamma_P, (20.0, 10, 17), 5.60605096173161688e-11, TEST_TOL6); |
354 |
|
TEST(gsl_cdf_gamma_P, (100.0, 10, 17), 2.70416097848011280e-2, TEST_TOL6); |
355 |
|
TEST(gsl_cdf_gamma_P, (1000.0, 10, 17), 1.000000000000000000e0, TEST_TOL6); |
356 |
|
TEST(gsl_cdf_gamma_P, (10000.0, 10, 17), 1.000000000000000000e0, TEST_TOL6); |
357 |
|
|
358 |
exit (gsl_test_summary()); |
exit (gsl_test_summary()); |
359 |
} |
} |