1 ## Copyright (C) 2008, 2009 Luca Favatella <slackydeb@gmail.com>
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.
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.
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/>.
16 ## Author: Luca Favatella <slackydeb@gmail.com>
19 function expectation = fitscalingrank (scores, nParents)
20 nr_scores = rows (scores);
21 r(1, 1:nr_scores) = ranks (scores(1:nr_scores, 1));
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;
32 ## number of input arguments
35 ## number of output arguments
43 %!shared scores, nParents, expectation
44 %! scores = rand (20, 1);
46 %! expectation = fitscalingrank (scores, nParents);
47 %!assert (sum (expectation), nParents, 1e-9);
49 %! [trash index_min_scores] = min (scores);
50 %! [trash index_max_expectation] = max (expectation);
51 %! assert (index_min_scores, index_max_expectation);
53 %! [trash index_max_scores] = max (scores);
54 %! [trash index_min_expectation] = min (expectation);
55 %! assert (index_max_scores, index_min_expectation);