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.200138037385 0.334558519179 0.125031231272];
  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.492838922640; 0.060105711640];
  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]);