X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Flanguage%2Fstats%2Fregression.q;h=4f6f0b862196b6a7365466147145714f3969bfa9;hb=1fe087cbcce6a9eb41928b7af6b58c7a2aac47f8;hp=1d31d1845e02e20064610bf3886cd8d11c71f5a0;hpb=73f67789df91a09ee91976434fb15c2ee1fb5e78;p=pspp

diff --git a/src/language/stats/regression.q b/src/language/stats/regression.q
index 1d31d1845e..4f6f0b8621 100644
--- a/src/language/stats/regression.q
+++ b/src/language/stats/regression.q
@@ -861,39 +861,6 @@ coeff_init (pspp_linreg_cache * c, struct design_matrix *dm)
   pspp_coeff_init (c->coeff, dm);
 }
 
-/*
-  Put the moments in the linreg cache.
- */
-static void
-compute_moments (pspp_linreg_cache * c, struct moments_var *mom,
-		 struct design_matrix *dm, size_t n)
-{
-  size_t i;
-  size_t j;
-  double weight;
-  double mean;
-  double variance;
-  double skewness;
-  double kurtosis;
-  /*
-     Scan the variable names in the columns of the design matrix.
-     When we find the variable we need, insert its mean in the cache.
-   */
-  for (i = 0; i < dm->m->size2; i++)
-    {
-      for (j = 0; j < n; j++)
-	{
-	  if (design_matrix_col_to_var (dm, i) == (mom + j)->v)
-	    {
-	      moments1_calculate ((mom + j)->m, &weight, &mean, &variance,
-				  &skewness, &kurtosis);
-	      pspp_linreg_set_indep_variable_mean (c, (mom + j)->v, mean);
-	      pspp_linreg_set_indep_variable_sd (c, (mom + j)->v, sqrt (variance));
-	    }
-	}
-    }
-}
-
 static bool
 run_regression (struct casereader *input, struct cmd_regression *cmd,
 		struct dataset *ds, pspp_linreg_cache **models)
@@ -956,9 +923,9 @@ run_regression (struct casereader *input, struct cmd_regression *cmd,
       n_indep = identify_indep_vars (indep_vars, dep_var);
       reader = casereader_clone (input);
       reader = casereader_create_filter_missing (reader, indep_vars, n_indep,
-						 MV_ANY, NULL);
+						 MV_ANY, NULL, NULL);
       reader = casereader_create_filter_missing (reader, &dep_var, 1,
-						 MV_ANY, NULL);
+						 MV_ANY, NULL, NULL);
       n_data = prepare_categories (casereader_clone (reader),
 				   indep_vars, n_indep, mom);
 
@@ -973,7 +940,8 @@ run_regression (struct casereader *input, struct cmd_regression *cmd,
 	    {
 	      lopts.get_indep_mean_std[i] = 1;
 	    }
-	  models[k] = pspp_linreg_cache_alloc (X->m->size1, X->m->size2);
+	  models[k] = pspp_linreg_cache_alloc (dep_var, (const struct variable **) indep_vars,
+					       X->m->size1, X->m->size2);
 	  models[k]->depvar = dep_var;
 	  /*
 	     For large data sets, use QR decomposition.