projects
/
pspp-builds.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Merge commit 'origin/stable'
[pspp-builds.git]
/
src
/
language
/
stats
/
chisquare.c
diff --git
a/src/language/stats/chisquare.c
b/src/language/stats/chisquare.c
index ed143824187593f984324bc39afeaf810c0080df..6cb7fc0b192ee2bfce5e50a620a26bd466a3ee6e 100644
(file)
--- a/
src/language/stats/chisquare.c
+++ b/
src/language/stats/chisquare.c
@@
-1,5
+1,5
@@
/* PSPP - a program for statistical analysis.
/* PSPP - a program for statistical analysis.
- Copyright (C) 2006, 2007 Free Software Foundation, Inc.
+ Copyright (C) 2006, 2007
, 2009
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
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
@@
-21,6
+21,7
@@
#include <stdlib.h>
#include <math.h>
#include <stdlib.h>
#include <math.h>
+#include <data/format.h>
#include <data/case.h>
#include <data/casereader.h>
#include <data/dictionary.h>
#include <data/case.h>
#include <data/casereader.h>
#include <data/dictionary.h>
@@
-57,7
+58,7
@@
create_freq_hash_with_range (const struct dictionary *dict,
{
bool warn = true;
float i_d;
{
bool warn = true;
float i_d;
- struct ccase c;
+ struct ccase
*
c;
struct hsh_table *freq_hash =
hsh_create (4, compare_freq, hash_freq,
struct hsh_table *freq_hash =
hsh_create (4, compare_freq, hash_freq,
@@
-78,21
+79,21
@@
create_freq_hash_with_range (const struct dictionary *dict,
hsh_insert (freq_hash, fr);
}
hsh_insert (freq_hash, fr);
}
- while (
casereader_read (input, &c)
)
+ while (
(c = casereader_read (input)) != NULL
)
{
union value obs_value;
struct freq **existing_fr;
struct freq *fr = xmalloc(sizeof (*fr));
{
union value obs_value;
struct freq **existing_fr;
struct freq *fr = xmalloc(sizeof (*fr));
- fr->value = case_data (
&
c, var);
+ fr->value = case_data (c, var);
- fr->count = dict_get_case_weight (dict,
&
c, &warn);
+ fr->count = dict_get_case_weight (dict, c, &warn);
obs_value.f = trunc (fr->value->f);
if ( obs_value.f < lo || obs_value.f > hi)
{
free (fr);
obs_value.f = trunc (fr->value->f);
if ( obs_value.f < lo || obs_value.f > hi)
{
free (fr);
- case_
destroy (&
c);
+ case_
unref (
c);
continue;
}
continue;
}
@@
-107,7
+108,7
@@
create_freq_hash_with_range (const struct dictionary *dict,
(*existing_fr)->count += fr->count;
free (fr);
(*existing_fr)->count += fr->count;
free (fr);
- case_
destroy (&
c);
+ case_
unref (
c);
}
if (casereader_destroy (input))
return freq_hash;
}
if (casereader_destroy (input))
return freq_hash;
@@
-130,20
+131,20
@@
create_freq_hash (const struct dictionary *dict,
const struct variable *var)
{
bool warn = true;
const struct variable *var)
{
bool warn = true;
- struct ccase c;
+ struct ccase
*
c;
struct hsh_table *freq_hash =
hsh_create (4, compare_freq, hash_freq,
free_freq_mutable_hash,
(void *) var);
struct hsh_table *freq_hash =
hsh_create (4, compare_freq, hash_freq,
free_freq_mutable_hash,
(void *) var);
- for (;
casereader_read (input, &c); case_destroy (&
c))
+ for (;
(c = casereader_read (input)) != NULL; case_unref (
c))
{
struct freq **existing_fr;
struct freq *fr = xmalloc(sizeof (*fr));
{
struct freq **existing_fr;
struct freq *fr = xmalloc(sizeof (*fr));
- fr->value = case_data (
&
c, var);
+ fr->value = case_data (c, var);
- fr->count = dict_get_case_weight (dict,
&
c, &warn);
+ fr->count = dict_get_case_weight (dict, c, &warn);
existing_fr = (struct freq **) hsh_probe (freq_hash, fr);
if ( *existing_fr)
existing_fr = (struct freq **) hsh_probe (freq_hash, fr);
if ( *existing_fr)
@@
-269,8
+270,8
@@
create_combo_frequency_table (const struct chisquare_test *test)
}
for ( i = test->lo ; i <= test->hi ; ++i )
}
for ( i = test->lo ; i <= test->hi ; ++i )
- tab_f
loat
(table, 0, 2 + i - test->lo,
- TAB_LEFT, 1 + i - test->lo, 8, 0);
+ tab_f
ixed
(table, 0, 2 + i - test->lo,
+
TAB_LEFT, 1 + i - test->lo, 8, 0);
tab_headers (table, 1, 0, 2, 0);
tab_headers (table, 1, 0, 2, 0);
@@
-330,6
+331,9
@@
chisquare_execute (const struct dataset *ds,
struct chisquare_test *cst = (struct chisquare_test *) test;
int n_cells = 0;
double total_expected = 0.0;
struct chisquare_test *cst = (struct chisquare_test *) test;
int n_cells = 0;
double total_expected = 0.0;
+ const struct variable *wvar = dict_get_weight (dict);
+ const struct fmt_spec *wfmt = wvar ?
+ var_get_print_format (wvar) : & F_8_0;
double *df = xzalloc (sizeof (*df) * ost->n_vars);
double *xsq = xzalloc (sizeof (*df) * ost->n_vars);
double *df = xzalloc (sizeof (*df) * ost->n_vars);
double *xsq = xzalloc (sizeof (*df) * ost->n_vars);
@@
-378,28
+382,28
@@
chisquare_execute (const struct dataset *ds,
/* The observed N */
/* The observed N */
- tab_
float
(freq_table, 1, i + 1, TAB_NONE,
- ff[i]->count,
8, 0
);
+ tab_
double
(freq_table, 1, i + 1, TAB_NONE,
+ ff[i]->count,
wfmt
);
if ( cst->n_expected > 0 )
exp = cst->expected[i] * total_obs / total_expected ;
else
exp = total_obs / (double) n_cells;
if ( cst->n_expected > 0 )
exp = cst->expected[i] * total_obs / total_expected ;
else
exp = total_obs / (double) n_cells;
- tab_
float
(freq_table, 2, i + 1, TAB_NONE,
- exp,
8, 2
);
+ tab_
double
(freq_table, 2, i + 1, TAB_NONE,
+ exp,
NULL
);
/* The residual */
/* The residual */
- tab_
float
(freq_table, 3, i + 1, TAB_NONE,
- ff[i]->count - exp,
8, 2
);
+ tab_
double
(freq_table, 3, i + 1, TAB_NONE,
+ ff[i]->count - exp,
NULL
);
xsq[v] += (ff[i]->count - exp) * (ff[i]->count - exp) / exp;
}
df[v] = n_cells - 1.0;
xsq[v] += (ff[i]->count - exp) * (ff[i]->count - exp) / exp;
}
df[v] = n_cells - 1.0;
- tab_
float
(freq_table, 1, i + 1, TAB_NONE,
- total_obs,
8, 0
);
+ tab_
double
(freq_table, 1, i + 1, TAB_NONE,
+ total_obs,
wfmt
);
tab_submit (freq_table);
tab_submit (freq_table);
@@
-450,8
+454,8
@@
chisquare_execute (const struct dataset *ds,
ds_destroy (&str);
/* The observed N */
ds_destroy (&str);
/* The observed N */
- tab_
float
(freq_table, v * 4 + 2, i + 2 , TAB_NONE,
- ff[i]->count,
8, 0
);
+ tab_
double
(freq_table, v * 4 + 2, i + 2 , TAB_NONE,
+ ff[i]->count,
wfmt
);
if ( cst->n_expected > 0 )
exp = cst->expected[i] * total_obs / total_expected ;
if ( cst->n_expected > 0 )
exp = cst->expected[i] * total_obs / total_expected ;
@@
-459,19
+463,19
@@
chisquare_execute (const struct dataset *ds,
exp = total_obs / (double) hsh_count (freq_hash);
/* The expected N */
exp = total_obs / (double) hsh_count (freq_hash);
/* The expected N */
- tab_
float
(freq_table, v * 4 + 3, i + 2 , TAB_NONE,
- exp,
8, 2
);
+ tab_
double
(freq_table, v * 4 + 3, i + 2 , TAB_NONE,
+ exp,
NULL
);
/* The residual */
/* The residual */
- tab_
float
(freq_table, v * 4 + 4, i + 2 , TAB_NONE,
- ff[i]->count - exp,
8, 2
);
+ tab_
double
(freq_table, v * 4 + 4, i + 2 , TAB_NONE,
+ ff[i]->count - exp,
NULL
);
xsq[v] += (ff[i]->count - exp) * (ff[i]->count - exp) / exp;
}
xsq[v] += (ff[i]->count - exp) * (ff[i]->count - exp) / exp;
}
- tab_
float
(freq_table, v * 4 + 2, tab_nr (freq_table) - 1, TAB_NONE,
- total_obs,
8, 0
);
+ tab_
double
(freq_table, v * 4 + 2, tab_nr (freq_table) - 1, TAB_NONE,
+ total_obs,
wfmt
);
df[v] = n_cells - 1.0;
df[v] = n_cells - 1.0;
@@
-494,11
+498,11
@@
chisquare_execute (const struct dataset *ds,
tab_text (stats_table, 1 + v, 0, TAB_CENTER, var_get_name (var));
tab_text (stats_table, 1 + v, 0, TAB_CENTER, var_get_name (var));
- tab_
float (stats_table, 1 + v, 1, TAB_NONE, xsq[v], 8,3
);
- tab_f
loat (stats_table, 1 + v, 2, TAB_NONE, df[v], 8,
0);
+ tab_
double (stats_table, 1 + v, 1, TAB_NONE, xsq[v], NULL
);
+ tab_f
ixed (stats_table, 1 + v, 2, TAB_NONE, df[v], 8,
0);
- tab_
float
(stats_table, 1 + v, 3, TAB_NONE,
- gsl_cdf_chisq_Q (xsq[v], df[v]),
8,3
);
+ tab_
double
(stats_table, 1 + v, 3, TAB_NONE,
+ gsl_cdf_chisq_Q (xsq[v], df[v]),
NULL
);
}
tab_submit (stats_table);
}
}
tab_submit (stats_table);
}