X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;ds=sidebyside;f=src%2Flanguage%2Fstats%2Fexamine.q;h=2649968b31485b105f6a80f3119ab39bdc851eea;hb=3c52a21736b3777ced3e47148fa2fc8f9eff59c8;hp=b564ba4415991aa1593f828fcdfc197635e6437d;hpb=a1efcf97ca2f75f4be6a0389ff2372c03ed2d4e1;p=pspp

diff --git a/src/language/stats/examine.q b/src/language/stats/examine.q
index b564ba4415..2649968b31 100644
--- a/src/language/stats/examine.q
+++ b/src/language/stats/examine.q
@@ -122,6 +122,9 @@ struct factor_metrics
   /* Sum of all weights, including those for missing values */
   double n;
 
+  /* Sum of weights of non_missing values */
+  double n_valid;
+
   double mean;
 
   double variance;
@@ -195,7 +198,7 @@ factor_destroy (struct xfactor *fctr)
 }
 
 static struct xfactor level0_factor;
-static struct ll_list factor_list = LL_INITIALIZER (factor_list);
+static struct ll_list factor_list;
 
 /* Parse the clause specifying the factors */
 static int examine_parse_independent_vars (struct lexer *lexer,
@@ -265,6 +268,8 @@ cmd_examine (struct lexer *lexer, struct dataset *ds)
   subc_list_double_create (&percentile_list);
   percentile_algorithm = PC_HAVERAGE;
 
+  ll_init (&factor_list);
+
   if ( !parse_examine (lexer, ds, &cmd, NULL) )
     {
       subc_list_double_destroy (&percentile_list);
@@ -934,23 +939,27 @@ examine_group (struct cmd_examine *cmd, struct casereader *reader, int level,
 	      case_data_idx (c, casereader_get_value_cnt (reader) - 1)->f;
 
 	  const double weight = wv ? case_data (c, wv)->f : 1.0;
+	  const union value *value = case_data (c, dependent_vars[v]);
 
 	  if (weight != SYSMIS)
 	    minimize (&result->metrics[v].cmin, weight);
 
 	  moments1_add (result->metrics[v].moments,
-			case_data (c, dependent_vars[v])->f,
+			value->f,
 			weight);
 
 	  result->metrics[v].n += weight;
 
+	  if ( ! var_is_value_missing (dependent_vars[v], value, MV_ANY) )
+	    result->metrics[v].n_valid += weight;
+
 	  extrema_add (result->metrics[v].maxima,
-		       case_data (c, dependent_vars[v])->f,
+		       value->f,
 		       weight,
 		       loc);
 
 	  extrema_add (result->metrics[v].minima,
-		       case_data (c, dependent_vars[v])->f,
+		       value->f,
 		       weight,
 		       loc);
 
@@ -985,7 +994,7 @@ examine_group (struct cmd_examine *cmd, struct casereader *reader, int level,
 	  for (i = 0 ; i < metric->n_ptiles; ++i)
 	    {
 	      metric->ptl[i] = (struct percentile *)
-		percentile_create (percentile_list.data[i] / 100.0, metric->n);
+		percentile_create (percentile_list.data[i] / 100.0, metric->n_valid);
 
 	      if ( percentile_list.data[i] == 25)
 		metric->quartiles[0] = metric->ptl[i];