[Pspp-cvs] pspp/src/language/stats regression.q

[Jason H Stover]

CVSROOT:        /sources/pspp
Module name:    pspp
Branch:         
Changes by:     Jason H Stover <address@hidden> 06/04/09 19:46:29

Modified files:
        src/language/stats: regression.q 

Log message:
        call residual function instead of predict to get residual

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

Patches:
Index: pspp/src/language/stats/regression.q
diff -u pspp/src/language/stats/regression.q:1.9 
pspp/src/language/stats/regression.q:1.10
--- pspp/src/language/stats/regression.q:1.9    Sun Apr  9 14:00:41 2006
+++ pspp/src/language/stats/regression.q        Sun Apr  9 19:46:29 2006
@@ -507,11 +507,13 @@
   int i;
   int case_num;
   double residual;
-  const union value **vals;
+  const union value **vals = NULL;
+  const union value *obs = NULL;
   struct casereader *r;
   struct ccase c;
 
   assert (lc != NULL);
+  assert (lc->depvar != NULL);
   assert (is_missing != NULL);
 
   if (save)
@@ -526,9 +528,13 @@
              for (i = 0; i < n_variables; ++i)
                {
                  vals[i] = case_data (&c, v_variables[i]->fv);
+                 if (v_variables[i]->index == lc->depvar->index)
+                   {
+                     obs = vals[i];
+                   }
                }
-             residual = (*lc->predict) ((const struct variable **) 
v_variables, 
-                                        (const union value **) vals, lc, 
n_variables);
+             residual = (*lc->residual) ((const struct variable **) 
v_variables, 
+                                         (const union value **) vals, obs, lc, 
n_variables);
            }
        }
       free (vals);