Had to get last call to multipass_split_output() inside
[pspp] / src / levene.c
index 3d171716b9afc3af7053f9ba6f33325a24fe99bb..7e0bd6a6f2ddc1547524ea914119f6b8adb9a8f6 100644 (file)
    02111-1307, USA. */
 
 #include <config.h>
-#include <assert.h>
+#include "error.h"
 #include "levene.h"
 #include "hash.h"
+#include "str.h"
 #include "var.h"
 #include "vfm.h"
 #include "alloc.h"
-#include "stats.h"
+#include "misc.h"
 
 #include <math.h>
 #include <stdlib.h>
@@ -347,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);
        }
     }
 
@@ -372,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 );
       
 
        }