* Bug #7858 fixed - Statistics: variance and variancef
[scilab.git] / scilab / modules / statistics / tests / unit_tests / variancef.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 x = [0.2113249 0.0002211 0.6653811; 0.7560439 0.9546254 0.6283918];
10 fre = [1 2 3; 3 4 3];
11 refM = 0.63623244375;
12 refV = 0.096057419504;
13 [v, m] = variancef(x, fre);
14 assert_checkalmostequal([v m], [refV refM]);
15 refM = [0.61986415 0.636490633333 0.64688645];
16 refV = [0.074179697240 0.242903351429 0.000410462494];
17 [v, m] = variancef(x, fre, "r");
18 assert_checkalmostequal([v m], [refV refM]);
19 refM = [0.367985066667; 0.79718087];
20 refV = [0.112075256021; 0.021077902385];
21 [v, m] = variancef(x, fre, "c");
22 assert_checkalmostequal([v m], [refV refM]);
23 // With the a priori mean
24 refM = 0.63623244375;
25 refV = 0.090053830785;
26 [v, m] = variancef(x, fre, "*", meanf(x,fre));
27 assert_checkalmostequal([v m], [refV refM]);
28 refM = [0.61986415 0.636490633333 0.64688645];
29 refV = [0.055634772930 0.202419459524 0.000342052079];
30 [v, m] = variancef(x, fre, "r", meanf(x,fre,"r"));
31 assert_checkalmostequal([v m], [refV refM]);
32 refM = [0.367985066667; 0.79718087];
33 refV = [0.093396046684; 0.018970112146];
34 [v, m] = variancef(x, fre, "c", meanf(x,fre,"c"));
35 assert_checkalmostequal([v m], [refV refM]);
36 // Biased variance
37 refM = 0.63623244375;
38 refV = 0.090053830785;
39 [v, m] = variancef(x, fre, "*", %nan);
40 assert_checkalmostequal([v m], [refV refM]);