Mercurial > hg > medcouple

comparison slowmedcouple.m @ 15:3958c6e77442

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add slow Octave medcouple implementation
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Fri, 16 Jan 2015 09:32:09 -0500
parents
children 6339e1a06aa9
comparison
equal deleted inserted replaced
14:549846d50ddb 15:3958c6e77442
1 z = load("gamma2");
2 n = length(z);
3 n2 = ceil(n/2);
4 z = sort(z, "descend");
5 zmed = median(z);
6 z -= zmed;
7 zplus = z(z >= 0);
8 zminus = z(z <= 0)';
9 zz = (zplus + zminus)./(zplus - zminus);
10 zz(isnan(zz)) = 0;
11 zzsort = sort(zz(:));
12 medc_idx = ceil(n2^2/2);
13 mc = zzsort(medc_idx);
14 printf("%.16g\n", mc);