Add multipass procedures. Add two-pass moments calculation. Rewrite

[pspp-builds.git] / src / levene.c

diff --git a/src/levene.c b/src/levene.c
index 635c68689e93af7d872b7d184756e8f3122669e2..7e0bd6a6f2ddc1547524ea914119f6b8adb9a8f6 100644 (file)
--- a/src/levene.c
+++ b/src/levene.c
@@ -27,7 +27,7 @@
 #include "var.h"
 #include "vfm.h"
 #include "alloc.h"
-#include "stats.h"
+#include "misc.h"
 
 #include <math.h>
 #include <stdlib.h>
@@ -348,7 +348,7 @@ levene2_calc (struct ccase *c, void *_l)
       if ( ! l->is_missing(v,var) )
        {
          levene_z = fabs(v->f - gs->mean); 
-         lz_denominator[i] += weight * sqr(levene_z - gs->lz_mean);
+         lz_denominator[i] += weight * pow2(levene_z - gs->lz_mean);
        }
     }
 
@@ -373,7 +373,7 @@ levene2_postcalc (void *_l)
          g = (struct group_statistics *) hsh_next(hash[v],&hi) )
        {
 
-         lz_numerator += g->n * sqr(g->lz_mean - lz[v].grand_mean );
+         lz_numerator += g->n * pow2(g->lz_mean - lz[v].grand_mean );
       
 
        }