X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;ds=sidebyside;f=src%2Flanguage%2Fstats%2Freliability.c;h=a4c4b862fdb344f98a3cf2a49cce3f2ee39265ba;hb=d5f42de822479581dded51ade0d997c16404890f;hp=144c31125b6f7fc4ca70c22b92bb79a7abfba284;hpb=04876a5077678930ea55110d1d2123c243203209;p=pspp diff --git a/src/language/stats/reliability.c b/src/language/stats/reliability.c index 144c31125b..a4c4b862fd 100644 --- a/src/language/stats/reliability.c +++ b/src/language/stats/reliability.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 2009, 2010, 2011 Free Software Foundation, Inc. + Copyright (C) 2009, 2010, 2011, 2013 Free Software Foundation, Inc. This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by @@ -113,15 +113,17 @@ reliability_destroy (struct reliability *rel) { int j; ds_destroy (&rel->scale_name); - for (j = 0; j < rel->n_sc ; ++j) - { - int x; - free (rel->sc[j].items); - moments1_destroy (rel->sc[j].total); - for (x = 0; x < rel->sc[j].n_items; ++x) - free (rel->sc[j].m[x]); - free (rel->sc[j].m); - } + if (rel->sc) + for (j = 0; j < rel->n_sc ; ++j) + { + int x; + free (rel->sc[j].items); + moments1_destroy (rel->sc[j].total); + if (rel->sc[j].m) + for (x = 0; x < rel->sc[j].n_items; ++x) + free (rel->sc[j].m[x]); + free (rel->sc[j].m); + } free (rel->sc); free (rel->variables); @@ -136,12 +138,14 @@ cmd_reliability (struct lexer *lexer, struct dataset *ds) reliability.n_variables = 0; reliability.variables = NULL; reliability.model = MODEL_ALPHA; - reliability.exclude = MV_ANY; + reliability.exclude = MV_ANY; reliability.summary = 0; - + reliability.n_sc = 0; + reliability.sc = NULL; reliability.wv = dict_get_weight (dict); - reliability.total_start = 0; + ds_init_empty (&reliability.scale_name); + lex_match (lexer, T_SLASH); @@ -168,7 +172,7 @@ cmd_reliability (struct lexer *lexer, struct dataset *ds) reliability.n_sc = 1; reliability.sc = xzalloc (sizeof (struct cronbach) * reliability.n_sc); - ds_init_cstr (&reliability.scale_name, "ANY"); + ds_assign_cstr (&reliability.scale_name, "ANY"); c = &reliability.sc[0]; c->n_items = reliability.n_variables; @@ -371,6 +375,20 @@ run_reliability (struct dataset *ds, const struct reliability *reliability) struct casereader *group; struct casegrouper *grouper = casegrouper_create_splits (proc_open (ds), dict); + int si; + + for (si = 0 ; si < reliability->n_sc; ++si) + { + struct cronbach *s = &reliability->sc[si]; + int i; + + s->m = xzalloc (sizeof *s->m * s->n_items); + s->total = moments1_create (MOMENT_VARIANCE); + + for (i = 0 ; i < s->n_items ; ++i ) + s->m[i] = moments1_create (MOMENT_VARIANCE); + } + while (casegrouper_get_next_group (grouper, &group)) { @@ -434,11 +452,10 @@ do_reliability (struct casereader *input, struct dataset *ds, { struct cronbach *s = &rel->sc[si]; - s->m = xzalloc (sizeof (s->m) * s->n_items); - s->total = moments1_create (MOMENT_VARIANCE); + moments1_clear (s->total); for (i = 0 ; i < s->n_items ; ++i ) - s->m[i] = moments1_create (MOMENT_VARIANCE); + moments1_clear (s->m[i]); } input = casereader_create_filter_missing (input,