* Bug #11891 fixed - Inaccurate Fisher ratio for one-way ANOVA
[scilab.git] / scilab / modules / statistics / tests / unit_tests / ftest.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
4 // Copyright (C) ????-2008 - INRIA
5 //
6 //  This file is distributed under the same license as the Scilab package.
7 // =============================================================================
8
9 // =============================================================================
10 // Tests for ftest() function
11 // =============================================================================
12
13 // Run with test_run('statistics', 'ftest', ['no_check_error_output']);
14
15 prec = 1.e-5;
16
17 samples = [46 55 54;
18          53 54 50;
19          49 58 51;
20          50 61 51;
21          46 52 49];
22
23 [f, p] = ftest(samples);
24
25 expected_p = 0.0055319;
26
27 assert_checktrue(norm(p-expected_p) < prec);
28
29 prec = 1.e-32;
30
31 T = [
32 107.8681568
33 107.8681465
34 107.8681572
35 107.8681785
36 107.8681446
37 107.8681903
38 107.8681526
39 107.8681494
40 107.8681616
41 107.8681587
42 107.8681519
43 107.8681486
44 107.8681419
45 107.8681569
46 107.8681508
47 107.8681672
48 107.8681385
49 107.8681518
50 107.8681662
51 107.8681424
52 107.8681360
53 107.8681333
54 107.8681610
55 107.8681477
56 ];
57
58 [f,p]=ftest([T,T*(1+1.e-6)]);
59 expected_p = 6.033347845799172588D-31;
60
61 assert_checktrue(norm(p-expected_p) < prec);