1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
//
// <-- CLI SHELL MODE -->
x = [0.2113249 0.0002211 0.6653811; 0.7560439 0.9546254 0.6283918];
fre = [1 2 3; 3 4 3];
refM = 0.63623244375;
refV = 0.096057419504;
[v, m] = variancef(x, fre);
assert_checkalmostequal([v m], [refV refM]);
refM = [0.61986415 0.636490633333 0.64688645];
refV = [0.074179697240 0.242903351429 0.000410462494];
[v, m] = variancef(x, fre, "r");
assert_checkalmostequal([v m], [refV refM]);
refM = [0.367985066667; 0.79718087];
refV = [0.112075256021; 0.021077902385];
[v, m] = variancef(x, fre, "c");
assert_checkalmostequal([v m], [refV refM]);
// With the a priori mean
refM = 0.63623244375;
refV = 0.090053830785;
[v, m] = variancef(x, fre, "*", meanf(x,fre));
assert_checkalmostequal([v m], [refV refM]);
refM = [0.61986415 0.636490633333 0.64688645];
refV = [0.055634772930 0.202419459524 0.000342052079];
[v, m] = variancef(x, fre, "r", meanf(x,fre,"r"));
assert_checkalmostequal([v m], [refV refM]);
refM = [0.367985066667; 0.79718087];
refV = [0.093396046684; 0.018970112146];
[v, m] = variancef(x, fre, "c", meanf(x,fre,"c"));
assert_checkalmostequal([v m], [refV refM]);
// Biased variance
refM = 0.63623244375;
refV = 0.090053830785;
[v, m] = variancef(x, fre, "*", %nan);
assert_checkalmostequal([v m], [refV refM]);
|
