Add tests for mean, gmean, hmean and wmean

run_tests.jl

@@ -6,7 +6,8 @@
 
 using Stats
 
-my_tests = ["test/01.jl"]
+my_tests = ["test/01.jl"
+            "test/means.jl"]
             #"test/statquiz.jl"]
 
 println("Running tests:")

test/means.jl

@@ -0,0 +1,22 @@
+using Stats
+using Base.Test
+
+@test_approx_eq mean([1, 2, 3]) 2.0
+@test_approx_eq mean(1:3) 2.0
+
+@test_approx_eq gmean([1, 2, 3]) (6.0)^(1/3)
+@test_approx_eq gmean(1:3) (6.0)^(1/3)
+@test_approx_eq gmean([2, 8]) 4.0
+@test_approx_eq gmean([4, 1, 1/32]) 0.5
+
+@test_approx_eq hmean([1, 2, 3]) 3 / (1 + 1/2 + 1/3)
+@test_approx_eq hmean(1:3) 3 / (1 + 1/2 + 1/3)
+@test_approx_eq hmean([1, 2, 4]) 12 / 7
+
+@test_approx_eq wmean([1, 2, 3], [1/3, 1/3, 1/3]) 2.0
+@test_approx_eq wmean([1, 2, 3], [1.0, 0.0, 0.0]) 1.0
+@test_approx_eq wmean([1, 2, 3], [0.0, 1.0, 0.0]) 2.0
+@test_approx_eq wmean([1, 2, 3], [0.0, 0.0, 1.0]) 3.0
+@test_approx_eq wmean([1, 2, 3], [0.5, 0.0, 0.5]) 2.0
+@test_approx_eq wmean([1, 2, 3], [0.5, 0.5, 0.0]) 1.5
+@test_approx_eq wmean([1, 2, 3], [0.0, 0.5, 0.5]) 2.5