diff --git a/CHANGES b/CHANGES
--- a/CHANGES
+++ b/CHANGES
@@ -100,3 +100,6 @@
 
 0.2.5.1:
 		added pcaReduceN to PCA
+
+0.2.5.2:
+		exposed pcaReduceN as requested by Tom Nielsen
diff --git a/hstatistics.cabal b/hstatistics.cabal
--- a/hstatistics.cabal
+++ b/hstatistics.cabal
@@ -1,5 +1,5 @@
 Name:               hstatistics
-Version:            0.2.5.1
+Version:            0.2.5.2
 License:            BSD3
 License-file:       LICENSE
 Copyright:          (c) A.V.H. McPhail 2010, 2011, 2012, 2013
diff --git a/lib/Numeric/Statistics.hs b/lib/Numeric/Statistics.hs
--- a/lib/Numeric/Statistics.hs
+++ b/lib/Numeric/Statistics.hs
@@ -31,8 +31,6 @@
 
 -----------------------------------------------------------------------------
 
---import Debug.Trace
-
 --import Numeric.Vector
 --import Numeric.Matrix
 --import Numeric.Container
@@ -132,10 +130,10 @@
 cut v c  = let c' = sort c
            in mapVector (\x -> cut_helper 0 x c') v 
     where
-      cut_helper i x c 
-          | i >= dim c                       = error "Numeric.Statistics: cut: data point not within interval"
-          | x >= (c @> i) && x <= (c @> (i+1)) = i
-          | otherwise                       = cut_helper (i + 1) x c
+      cut_helper j x d 
+          | j >= dim d                       = error "Numeric.Statistics: cut: data point not within interval"
+          | x >= (d @> j) && x <= (d @> (j+1)) = j
+          | otherwise                       = cut_helper (j + 1) x d
 
 -----------------------------------------------------------------------------
 
diff --git a/lib/Numeric/Statistics/ICA.hs b/lib/Numeric/Statistics/ICA.hs
--- a/lib/Numeric/Statistics/ICA.hs
+++ b/lib/Numeric/Statistics/ICA.hs
@@ -2,7 +2,7 @@
 -- |
 -- Module      :  Numeric.Statistics.ICA
 -- Copyright   :  (c) A. V. H. McPhail 2010
--- License     :  GPL-style
+-- License     :  BSD3
 --
 -- Maintainer  :  haskell.vivian.mcphail <at> gmail <dot> com
 -- Stability   :  provisional
@@ -176,11 +176,11 @@
     -> Int                         -- ^ sampling size (must be smaller than length of data)
     -> I.Array Int (Vector Double) -- ^ data
     -> (I.Array Int (Vector Double),Matrix Double) -- ^ transformed data, ica transform
-ica r g g' n t s a = let i = I.rangeSize $ I.bounds a
-                         w = random_vector r (i,i)
+ica r g g' n t s a = let j = I.rangeSize $ I.bounds a
+                         w = random_vector r (j,j)
                          x' = fromRows $ I.elems a
                          -- next line is BAD if distribution not stationary
-                         x = concat $ toBlocksEvery i s x'
+                         x = concat $ toBlocksEvery j s x'
                          w' = ica' g g' n t w x
                          y = w' <> x'
                      in (I.listArray (1,1) $ toRows y,w') 
diff --git a/lib/Numeric/Statistics/Information.hs b/lib/Numeric/Statistics/Information.hs
--- a/lib/Numeric/Statistics/Information.hs
+++ b/lib/Numeric/Statistics/Information.hs
@@ -3,7 +3,7 @@
 -- |
 -- Module      :  Numeric.Statistics.Information
 -- Copyright   :  (c) A. V. H. McPhail 2010
--- License     :  GPL-style
+-- License     :  BSD3
 --
 -- Maintainer  :  haskell.vivian.mcphail <at> gmail <dot> com
 -- Stability   :  provisional
diff --git a/lib/Numeric/Statistics/PCA.hs b/lib/Numeric/Statistics/PCA.hs
--- a/lib/Numeric/Statistics/PCA.hs
+++ b/lib/Numeric/Statistics/PCA.hs
@@ -2,7 +2,7 @@
 -- |
 -- Module      :  Numeric.Statistics.PCA
 -- Copyright   :  (c) A. V. H. McPhail 2010
--- License     :  GPL-style
+-- License     :  BSD3
 --
 -- Maintainer  :  haskell.vivian.mcphail <at> gmail <dot> com
 -- Stability   :  provisional
@@ -13,7 +13,7 @@
 -----------------------------------------------------------------------------
 
 module Numeric.Statistics.PCA (
-                               pca, pcaN, pcaTransform, pcaReduce
+                               pca, pcaN, pcaTransform, pcaReduce, pcaReduceN
                           ) where
 
 
diff --git a/lib/Numeric/Statistics/PDF.hs b/lib/Numeric/Statistics/PDF.hs
--- a/lib/Numeric/Statistics/PDF.hs
+++ b/lib/Numeric/Statistics/PDF.hs
@@ -4,7 +4,7 @@
 -- |
 -- Module      :  Numeric.Statistics.PDF
 -- Copyright   :  (c) A. V. H. McPhail 2010
--- License     :  GPL-style
+-- License     :  BSD3
 --
 -- Maintainer  :  haskell.vivian.mcphail <at> gmail <dot> com
 -- Stability   :  provisional
diff --git a/lib/Numeric/Statistics/Surrogate.hs b/lib/Numeric/Statistics/Surrogate.hs
--- a/lib/Numeric/Statistics/Surrogate.hs
+++ b/lib/Numeric/Statistics/Surrogate.hs
@@ -2,7 +2,7 @@
 -- |
 -- Module      :  Numeric.Statistics.Surrogate
 -- Copyright   :  (c) Alexander Vivian Hugh McPhail 2010
--- License     :  GPL-style
+-- License     :  BSD3
 --
 -- Maintainer  :  haskell.vivian.mcphail <at> gmail <dot> com
 -- Stability   :  provisional