hstatistics 0.2.0.7 → 0.2.0.8
raw patch · 3 files changed
+12/−6 lines, 3 filesdep ~hmatrix-gsl-stats
Dependency ranges changed: hmatrix-gsl-stats
Files
- CHANGES +7/−1
- hstatistics.cabal +2/−2
- lib/Numeric/Statistics/ICA.hs +3/−3
CHANGES view
@@ -36,7 +36,13 @@ added PCA added ICA -0.2.0.7+0.2.0.7: removed hard-coded sigmoid from ICA algorithm internals was not using seed passed to ica minor fixes to remove warnings++0.2.0.8:++ error in decorrelation fixed++0.2.0.9:
hstatistics.cabal view
@@ -1,5 +1,5 @@ Name: hstatistics-Version: 0.2.0.7+Version: 0.2.0.8 License: GPL License-file: LICENSE Copyright: (c) A.V.H. McPhail 2010@@ -23,7 +23,7 @@ Build-Depends: base >= 3 && < 5, array, random,- hmatrix >= 0.10.0, hmatrix-gsl-stats >= 0.1.1.5+ hmatrix >= 0.10.0, hmatrix-gsl-stats >= 0.1.2.1 Extensions:
lib/Numeric/Statistics/ICA.hs view
@@ -141,9 +141,9 @@ decorrelate :: Matrix Double -> Matrix Double decorrelate w = let w' = w / (scalar $ sqrt $ pnorm PNorm2 (w <> trans w)) in decorrelate' w w'- where decorrelate' w w' - | converged 0.000001 w w' = w'- | otherwise = decorrelate' w' ((scale 1.5 w') - (scale 0.5 (w <> trans w <> w)))+ where decorrelate' m m' + | converged 0.000001 m m' = m'+ | otherwise = decorrelate' m' ((scale 1.5 m') - (scale 0.5 (m' <> trans m' <> m'))) {- don't know how to do svd of non-square matrices decorrelate m = let (u,d,v) = svd m in u <> (diag (d ** (-0.5))) <> trans v <> m