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

[Jason H Stover]

CVSROOT:        /sources/pspp
Module name:    pspp
Branch:         
Changes by:     Jason H Stover <address@hidden> 06/04/17 16:26:45

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

Log message:
        added transformation to compute residuals

CVSWeb URLs:
http://cvs.savannah.gnu.org/viewcvs/pspp/pspp/src/language/stats/ChangeLog.diff?tr1=1.6&tr2=1.7&r1=text&r2=text
http://cvs.savannah.gnu.org/viewcvs/pspp/pspp/src/language/stats/regression.q.diff?tr1=1.11&tr2=1.12&r1=text&r2=text

Patches:
Index: pspp/src/language/stats/ChangeLog
diff -u pspp/src/language/stats/ChangeLog:1.6 
pspp/src/language/stats/ChangeLog:1.7
--- pspp/src/language/stats/ChangeLog:1.6       Sun Apr 16 01:05:15 2006
+++ pspp/src/language/stats/ChangeLog   Mon Apr 17 16:26:44 2006
@@ -1,3 +1,16 @@
+2006-04-17  Jason Stover  <address@hidden>
+
+       * regression.q (regression_trns_proc): New function.
+
+       * regression.q (subcommand_save): Create variable residuals and
+       add a transformation to assign values to them. Also free the
+       linreg_cache if the SAVE command was not called. Removed the
+       casereading loop. Placed actual computation of residuals in
+       regression_trns_proc.
+
+       * regression.q (run_regression): Moved call to free
+       pspp_linreg_cache to subcommand_save.
+
 Sat Apr 15 18:01:03 2006  Ben Pfaff  <address@hidden>
 
        * examine.q (output_examine): Add casts to fix warnings.
Index: pspp/src/language/stats/regression.q
diff -u pspp/src/language/stats/regression.q:1.11 
pspp/src/language/stats/regression.q:1.12
--- pspp/src/language/stats/regression.q:1.11   Mon Apr 17 01:54:15 2006
+++ pspp/src/language/stats/regression.q        Mon Apr 17 16:26:44 2006
@@ -501,43 +501,66 @@
   statistics_keyword_output (reg_stats_tol, keywords[tol], c);
   statistics_keyword_output (reg_stats_selection, keywords[selection], c);
 }
-static void
-subcommand_save (int save, pspp_linreg_cache *lc, const struct casefile *cf, 
int *is_missing)
+static int
+regression_trns_proc (void *m, struct ccase *c, size_t case_idx)
 {
-  int i;
-  int case_num;
-  double residual;
+  size_t i;
+  size_t n_vars;
+  pspp_linreg_cache *model = m;
+  union value *output;
   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);
+  struct variable **vars = NULL;
+  struct variable *resid = NULL;
+  
+  assert (model != NULL);
+  assert (model->depvar != NULL);
+  
+  vals = xnmalloc (n_variables, sizeof (*vals));
+  dict_get_vars (default_dict, &vars, &n_vars, 1u << DC_ORDINARY);
+  resid = dict_get_var (default_dict, case_idx);
+  assert (resid != NULL);
+  output = case_data_rw (c, resid->fv);
+  assert (output);
 
-  if (save)
+  for (i = 0; i < n_vars; i++)
     {
-      vals = xnmalloc (n_variables, sizeof (*vals));
-      for (r = casefile_get_reader (cf); casereader_read (r, &c);
-          case_destroy (&c))
+      if (i != case_idx) /* Do not use the residual variable as a predictor. */
        {
-         case_num = casereader_cnum (r) - 1;
-         if (!is_missing[case_num])
+         if (vars[i]->index == model->depvar->index) 
            {
-             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->residual) ((const struct variable **) 
v_variables, 
-                                         (const union value **) vals, obs, lc, 
n_variables);
+             /* Do not use the dependent variable as a predictor. */
+             obs = case_data (c, i);
+             assert (obs != NULL);
+           }
+         else
+           {
+             vals[i] = case_data (c, i);
            }
        }
-      free (vals);
+    }
+  output->f = (*model->residual) ((const struct variable **) vars, 
+                                 vals, obs, model, i);
+  free (vals);
+  return TRNS_CONTINUE;
+}
+static void
+subcommand_save (int save, pspp_linreg_cache *lc)
+{
+  struct variable *residuals = NULL;
+
+  assert (lc != NULL);
+  assert (lc->depvar != NULL);
+
+  if (save)
+    {
+      residuals = dict_create_var_assert (default_dict, "residuals", 0);
+      assert (residuals != NULL);
+      add_transformation (regression_trns_proc, pspp_linreg_cache_free, lc);
+    }
+  else 
+    {
+      pspp_linreg_cache_free (lc);
     }
 }
 static int
@@ -1042,12 +1065,11 @@
        */
       pspp_linreg ((const gsl_vector *) Y, X->m, &lopts, lcache);
       subcommand_statistics (cmd.a_statistics, lcache);
-      subcommand_save (cmd.sbc_save, lcache, cf, is_missing_case);
       subcommand_export (cmd.sbc_export, lcache);
+      subcommand_save (cmd.sbc_save, lcache);
       gsl_vector_free (Y);
       design_matrix_destroy (X);
       free (indep_vars);
-      pspp_linreg_cache_free (lcache);
       free (lopts.get_indep_mean_std);
       casereader_destroy (r);
     }