[Pspp-cvs] pspp/src regression.q

[Jason H Stover]

CVSROOT:        /sources/pspp
Module name:    pspp
Branch:         
Changes by:     Jason H Stover <address@hidden> 06/01/04 18:02:14

Modified files:
        src            : regression.q 

Log message:
        Fixed centering prior to sweep operation

CVSWeb URLs:
http://cvs.savannah.gnu.org/viewcvs/pspp/pspp/src/regression.q.diff?tr1=1.32&tr2=1.33&r1=text&r2=text

Patches:
Index: pspp/src/regression.q
diff -u pspp/src/regression.q:1.32 pspp/src/regression.q:1.33
--- pspp/src/regression.q:1.32  Wed Jan  4 16:33:04 2006
+++ pspp/src/regression.q       Wed Jan  4 18:02:14 2006
@@ -822,18 +822,22 @@
       X =
        design_matrix_create (n_indep, (const struct variable **) indep_vars,
                              n_data);
+      for (i = 0; i < X->m->size2; i++)
+       {
+         lopts.get_indep_mean_std[i] = 1;
+       }
       lcache = pspp_linreg_cache_alloc (X->m->size1, X->m->size2);
       lcache->indep_means = gsl_vector_alloc (X->m->size2);
       lcache->indep_std = gsl_vector_alloc (X->m->size2);
       lcache->depvar = (const struct variable *) depvar;
       /*
-         For large data sets, use QR decomposition.
-       */
+       For large data sets, use QR decomposition.
+      */
       if (n_data > sqrt (n_indep) && n_data > REG_LARGE_DATA)
        {
          lcache->method = PSPP_LINREG_SVD;
        }
-
+      
       /*
          The second pass creates the design matrix.
        */
@@ -869,7 +873,6 @@
                        {
                          design_matrix_set_numeric (X, row, v, val);
                        }
-                     lopts.get_indep_mean_std[i] = 1;
                    }
                }
              val = case_data (&c, depvar->fv);