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