projects
/
pspp
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Fixed a bug in the levene test, and added the levene test to the oneway cmd
[pspp]
/
src
/
levene.c
diff --git
a/src/levene.c
b/src/levene.c
index 4659f199ac018f98bd46886c2ad61d917ef24a88..e29e77629c110000b250a07667475e523522cda5 100644
(file)
--- a/
src/levene.c
+++ b/
src/levene.c
@@
-232,8
+232,8
@@
levene_postcalc (void *_l)
for (v = 0; v < l->n_dep; ++v)
{
for (v = 0; v < l->n_dep; ++v)
{
+ /* This is Z_LL */
lz[v].grand_mean = lz[v].grand_total / lz[v].total_n ;
lz[v].grand_mean = lz[v].grand_total / lz[v].total_n ;
-
}
}
}
}
@@
-344,17
+344,15
@@
levene2_postcalc (void *_l)
g != 0 ;
g = (struct group_statistics *) hsh_next(hash,&hi) )
{
g != 0 ;
g = (struct group_statistics *) hsh_next(hash,&hi) )
{
-
lz_numerator += g->n * pow2(g->lz_mean - lz[v].grand_mean );
lz_numerator += g->n * pow2(g->lz_mean - lz[v].grand_mean );
-
-
}
lz_numerator *= ( l->v_dep[v]->p.grp_data.ugs.n -
l->v_dep[v]->p.grp_data.n_groups );
}
lz_numerator *= ( l->v_dep[v]->p.grp_data.ugs.n -
l->v_dep[v]->p.grp_data.n_groups );
- lz_denominator[v]
/
= (l->v_dep[v]->p.grp_data.n_groups - 1);
+ lz_denominator[v]
*
= (l->v_dep[v]->p.grp_data.n_groups - 1);
l->v_dep[v]->p.grp_data.levene = lz_numerator/lz_denominator[v] ;
l->v_dep[v]->p.grp_data.levene = lz_numerator/lz_denominator[v] ;
+
}
/* Now clear up after ourselves */
}
/* Now clear up after ourselves */