bugfix on P scale

2007-01-31 13:00:36 +00:00 · 2007-01-31 13:00:36 +00:00 · 37feb1c25b
parent 2cfa3ca415
commit 37feb1c25b
1 changed files with 1 additions and 1 deletions
--- a/fluents/lib/blmfuncs.py
+++ b/fluents/lib/blmfuncs.py
@ -77,7 +77,7 @@ class PCA(Model):
        jk_segments = pca_jkP(self.model['E0'], aopt, n_sets)
        Pcal = self.model['P'][:,:aopt]
        # add the scale to P
-        tnorm = scipy.apply_along_axis(norm, 0, self.model['T'])
+        tnorm = scipy.apply_along_axis(norm, 0, self.model['T'][:,:aopt])
        Pcal = Pcal*tnorm
        tsq = hotelling(jk_segments, Pcal, p_center,
                        cov_center, alpha, crot, strict)