* Bug #11891 fixed - Inaccurate Fisher ratio for one-way ANOVA
[scilab.git] / scilab / modules / statistics / tests / unit_tests / ftuneq.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 ftuneq() function
11 // =============================================================================
12
13 // Run with test_run('statistics', 'ftuneq', ['no_check_error_output']);
14
15 prec = 1.e-5;
16
17 sample1 = [46 55 54];
18 sample2 = [53 54];
19 sample3 = [50 49 58 51 50];
20 sample4 = [61 51 46 52];
21
22 [f, p] = ftuneq(sample1, sample2, sample3, sample4);
23
24 expected_p = 0.9602320;
25
26 assert_checktrue(norm(p-expected_p) < prec);
27
28 prec = 1.e-32;
29
30 T = [
31 107.8681568
32 107.8681465
33 107.8681572
34 107.8681785
35 107.8681446
36 107.8681903
37 107.8681526
38 107.8681494
39 107.8681616
40 107.8681587
41 107.8681519
42 107.8681486
43 107.8681419
44 107.8681569
45 107.8681508
46 107.8681672
47 107.8681385
48 107.8681518
49 107.8681662
50 107.8681424
51 107.8681360
52 107.8681333
53 107.8681610
54 107.8681477
55 ];
56
57 [f,p]=ftuneq(T,T*(1+1.e-6));
58 expected_p = 6.033347845799172588D-31;
59
60 assert_checktrue(norm(p-expected_p) < prec);