Mercurial > hg > medcouple

changeset 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 549846d50ddb
children 9b6502f3f3bc
files slowmedcouple.m
diffstat 1 files changed, 14 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
new file mode 100644
--- /dev/null
+++ b/slowmedcouple.m
@@ -0,0 +1,14 @@
+z = load("gamma2");
+n = length(z);
+n2 = ceil(n/2);
+z = sort(z, "descend");
+zmed = median(z);
+z -= zmed;
+zplus = z(z >= 0);
+zminus = z(z <= 0)';
+zz = (zplus + zminus)./(zplus - zminus);
+zz(isnan(zz)) = 0;
+zzsort = sort(zz(:));
+medc_idx = ceil(n2^2/2);
+mc = zzsort(medc_idx);
+printf("%.16g\n", mc);