octave_packages/ga-0.10.0/fitscalingrank.m

   1 ## Copyright (C) 2008, 2009 Luca Favatella <slackydeb@gmail.com>
   2 ##
   3 ## This program is free software; you can redistribute it and/or modify
   4 ## it under the terms of the GNU General Public License as published by
   5 ## the Free Software Foundation; either version 2 of the License, or
   6 ## (at your option) any later version.
   7 ##
   8 ## This program is distributed in the hope that it will be useful,
   9 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11 ## GNU General Public License for more details.
  12 ##
  13 ## You should have received a copy of the GNU General Public License
  14 ## along with this program; If not, see <http://www.gnu.org/licenses/>.
  15
  16 ## Author: Luca Favatella <slackydeb@gmail.com>
  17 ## Version: 1.3.3
  18
  19 function expectation = fitscalingrank (scores, nParents)
  20   nr_scores = rows (scores);
  21   r(1, 1:nr_scores) = ranks (scores(1:nr_scores, 1));
  22                                 #TODO
  23                                 #ranks ([7,2,2]) == [3.0,1.5,1.5]
  24                                 #is [3,1,2] (or [3,2,1]) useful?
  25   expectation_wo_nParents(1, 1:nr_scores) = arrayfun (@(n) 1 / sqrt (n), r);
  26   expectation(1, 1:nr_scores) = \
  27       (nParents / sum (expectation_wo_nParents)) * \
  28       expectation_wo_nParents;
  29 endfunction
  30
  31
  32 ## number of input arguments
  33 # TODO
  34
  35 ## number of output arguments
  36 # TODO
  37
  38 ## type of arguments
  39 # TODO
  40
  41 # TODO
  42
  43 %!shared scores, nParents, expectation
  44 %! scores = rand (20, 1);
  45 %! nParents = 32;
  46 %! expectation = fitscalingrank (scores, nParents);
  47 %!assert (sum (expectation), nParents, 1e-9);
  48 %!test
  49 %! [trash index_min_scores] = min (scores);
  50 %! [trash index_max_expectation] = max (expectation);
  51 %! assert (index_min_scores, index_max_expectation);
  52 %!test
  53 %! [trash index_max_scores] = max (scores);
  54 %! [trash index_min_expectation] = min (expectation);
  55 %! assert (index_max_scores, index_min_expectation);