From 87686a9a3174e7fe936d668210991932f1b0a9fd Mon Sep 17 00:00:00 2001
From: John Benediktsson <mrjbq7@gmail.com>
Date: Fri, 4 May 2012 14:43:29 -0700
Subject: [PATCH] math.statistics: adding contraharmonic-mean.

---
 basis/math/statistics/statistics-tests.factor | 1 +
 basis/math/statistics/statistics.factor       | 3 +++
 2 files changed, 4 insertions(+)

diff --git a/basis/math/statistics/statistics-tests.factor b/basis/math/statistics/statistics-tests.factor
index de5ffa84ac..8e28a463f7 100644
--- a/basis/math/statistics/statistics-tests.factor
+++ b/basis/math/statistics/statistics-tests.factor
@@ -8,6 +8,7 @@ IN: math.statistics.tests
 [ t ] [ { 2 2 2 2 } geometric-mean 2.0 .0001 ~ ] unit-test
 [ 1.0 ] [ { 1 1 1 } geometric-mean ] unit-test
 [ 1/3 ] [ { 1 1 1 } harmonic-mean ] unit-test
+[ 5+1/4 ] [ { 1 3 5 7 } contraharmonic-mean ] unit-test
 
 [ 0 ] [ { 1 } range ] unit-test
 [ 89 ] [ { 1 2 30 90 } range ] unit-test
diff --git a/basis/math/statistics/statistics.factor b/basis/math/statistics/statistics.factor
index dad1f75ce6..1ab5760fa6 100644
--- a/basis/math/statistics/statistics.factor
+++ b/basis/math/statistics/statistics.factor
@@ -14,6 +14,9 @@ IN: math.statistics
 : harmonic-mean ( seq -- x )
     [ recip ] map-sum recip ;
 
+: contraharmonic-mean ( seq -- x )
+    [ [ sq ] map-sum ] [ sum ] bi / ;
+
 <PRIVATE
 
 :: ((kth-object)) ( seq k nth-quot exchange-quot quot: ( x y -- ? ) -- elt )