Add English imposed sentence in tests. test_run(statistics, [bug_7858, bug_13050])
[scilab.git] / scilab / modules / statistics / tests / nonreg_tests / bug_7858.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 //
8 // <-- CLI SHELL MODE -->
9 //
10 // <-- ENGLISH IMPOSED -->
11 //
12 // <-- Non-regression test for bug 7858 -->
13 //
14 // <-- Bugzilla URL -->
15 // http://bugzilla.scilab.org/show_bug.cgi?id=7858
16 //
17 // <-- Short Description -->
18 // variance and variancef can take now take the a priori mean as input, and it can be a scalar.
19 // =============================================================================
20 //--------------------------------------------------------------
21 // Variance
22 x = [
23 0.9    0.7
24 0.1    0.1
25 0.5    0.4 ];
26 orien = "r";
27 expectedM = [1 1];
28 expectedV = [1.07 1.26]/3;
29 // Voluntarily passing a bad a priori mean, to check that we get the same results with scalar
30 [Variance, Mean] = variance( x, orien, ones(mean(x,orien)) );
31 assert_checkalmostequal([Variance Mean], [expectedV expectedM]);
32 // With a scalar
33 [Variance, Mean] = variance( x, orien, 1 );
34 WARNING: variance: The significance of input argument #3 has been modified. Please refer to the variance help page.
35 assert_checkalmostequal([Variance Mean], [expectedV expectedM]);
36 orien = "c";
37 expectedM = [1; 1; 1];
38 expectedV = [0.05; 0.81; 0.305];
39 // Voluntarily passing a bad a priori mean, to check that we get the same results with scalar
40 [Variance, Mean] = variance( x, orien, ones(mean(x,orien)) );
41 assert_checkalmostequal([Variance Mean], [expectedV expectedM]);
42 // With a scalar
43 [Variance, Mean] = variance( x, orien, 1 );
44 WARNING: variance: The significance of input argument #3 has been modified. Please refer to the variance help page.
45 assert_checkalmostequal([Variance Mean], [expectedV expectedM]);
46 //--------------------------------------------------------------
47 // Variancef
48 x = [0.2113249 0.0002211 0.6653811; 0.7560439 0.9546254 0.6283918];
49 fre = [1 2 3; 3 4 3];
50 orien = "r";
51 refM = [1 1 1];
52 refV = [0.889522663062 0.593015108708 0.593015108708];
53 [v, m] = variancef( x, fre, orien, ones(meanf(x,fre,orien)) );
54 assert_checkalmostequal([v m], [refV refM]);
55 [v, m] = variancef( x, fre, orien, 1 );
56 assert_checkalmostequal([v m], [refV refM]);
57 orien = "c";
58 refM = [1; 1];
59 refV = [0.288922678414; 0.019966608736];
60 [v, m] = variancef( x, fre, orien, ones(meanf(x,fre,orien)) );
61 assert_checkalmostequal([v m], [refV refM]);
62 [v, m] = variancef( x, fre, orien, 1 );
63 assert_checkalmostequal([v m], [refV refM]);