6e2b070d4b80ac4d35d41cdb4b8ad3d4d27d8115
[scilab.git] / scilab / modules / statistics / tests / nonreg_tests / bug_7593.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 // <-- Non-regression test for bug 7593 -->
11 //
12 // <-- Bugzilla URL -->
13 // http://bugzilla.scilab.org/show_bug.cgi?id=7593
14 //
15 // <-- Short Description -->
16 // The function st_deviation was merged into stdev (they were duplicates),
17 // and msd is now accessible through the third input argument of stdev
18 x = matrix([5 2:8], 2, 4);
19 assert_checkalmostequal(stdev(x,"*",mean(x)), 1.870828693387);
20 assert_checkalmostequal(norm(stdev(x,1,mean(x,"r"))-[3 1 1 1]/2),   0);
21 assert_checkalmostequal(norm(stdev(x,"r",mean(x,"r"))-[3 1 1 1]/2), 0);
22 assert_checkalmostequal(norm(stdev(x,2,mean(x,"c"))-sqrt([2;5])),   0);
23 assert_checkalmostequal(norm(stdev(x,"c",mean(x,"c"))-sqrt([2;5])), 0);
24 // Biased deviation
25 assert_checkalmostequal(stdev(x,"*",%nan), 1.870828693387);
26 assert_checkalmostequal(norm(stdev(x,1,%nan)-[3 1 1 1]/2),   0);
27 assert_checkalmostequal(norm(stdev(x,"r",%nan)-[3 1 1 1]/2), 0);
28 assert_checkalmostequal(norm(stdev(x,2,%nan)-sqrt([2;5])),   0);
29 assert_checkalmostequal(norm(stdev(x,"c",%nan)-sqrt([2;5])), 0);
30 // With scalar means
31 rand("seed", 0);
32 x = rand(3, 10) + 0.5;
33 refR = [
34 0.364541128374
35 0.155589477734
36 0.316156843180
37 0.268553719068
38 0.219155131832
39 0.273506665569
40 0.262844435686
41 0.318334990307
42 0.149138182829
43 0.111372370623 ]';
44 refC = [
45 0.299608392744
46 0.216862008097
47 0.245620889607 ];
48 assert_checkalmostequal(stdev(x, "*", 1), 0.256335678990);
49 assert_checkalmostequal(stdev(x, 1, 1),   refR);
50 assert_checkalmostequal(stdev(x, "r", 1), refR);
51 assert_checkalmostequal(stdev(x, 2, 1),   refC);
52 assert_checkalmostequal(stdev(x, "c", 1), refC);
53 // Error checks
54 refMsg = msprintf(_("%s: Wrong size for input argument #%d.\n"),"stdev",3);
55 assert_checkerror("stdev(x, ""*"", [1 1])",   refMsg);
56 refMsg = msprintf(_("%s: Wrong type for input argument #%d: A real matrix expected.\n"),"stdev",3);
57 assert_checkerror("stdev(x, ""*"", ""msd"")", refMsg);