02110-1301, USA. */
#include <config.h>
+
#include <gsl/gsl_cdf.h>
-#include "message.h"
+#include <math.h>
#include <stdio.h>
#include <stdlib.h>
-#include <math.h>
-#include "alloc.h"
-#include "str.h"
-#include "case.h"
-#include "dictionary.h"
-#include "command.h"
-#include "lexer.h"
-#include "message.h"
-#include "magic.h"
-#include "misc.h"
-#include "table.h"
-#include "manager.h"
-#include "value-labels.h"
-#include "variable.h"
-#include "procedure.h"
-#include "hash.h"
-#include "casefile.h"
-#include "group-proc.h"
-#include "group.h"
-#include "levene.h"
+
+#include <data/case.h>
+#include <data/casefile.h>
+#include <data/dictionary.h>
+#include <data/procedure.h>
+#include <data/value-labels.h>
+#include <data/variable.h>
+#include <data/casefilter.h>
+#include <language/command.h>
+#include <language/dictionary/split-file.h>
+#include <language/lexer/lexer.h>
+#include <libpspp/alloc.h>
+#include <libpspp/compiler.h>
+#include <libpspp/hash.h>
+#include <libpspp/magic.h>
+#include <libpspp/message.h>
+#include <libpspp/message.h>
+#include <libpspp/misc.h>
+#include <libpspp/str.h>
+#include <math/group-proc.h>
+#include <math/group.h>
+#include <math/levene.h>
+#include <output/manager.h>
+#include <output/table.h>
+#include "sort-criteria.h"
#include "gettext.h"
#define _(msgid) gettext (msgid)
/* (specification)
"ONEWAY" (oneway_):
*^variables=custom;
- +missing=miss:!analysis/listwise,
- incl:include/!exclude;
- contrast= double list;
- statistics[st_]=descriptives,homogeneity.
+ missing=miss:!analysis/listwise,
+ incl:include/!exclude;
+ +contrast= double list;
+ +statistics[st_]=descriptives,homogeneity.
*/
/* (declarations) */
/* (functions) */
-
-
-static int bad_weight_warn = 1;
+static bool bad_weight_warn = true;
static struct cmd_oneway cmd;
/* The number of distinct values of the independent variable, when all
missing values are disregarded */
-static int ostensible_number_of_groups=-1;
-
-
-/* Function to use for testing for missing values */
-static is_missing_func *value_is_missing;
+static int ostensible_number_of_groups = -1;
-static bool run_oneway(const struct casefile *cf, void *_mode);
+static bool run_oneway(const struct ccase *first,
+ const struct casefile *cf,
+ void *_mode, const struct dataset *);
/* Routines to show the output tables */
int
-cmd_oneway(void)
+cmd_oneway (struct dataset *ds)
{
int i;
bool ok;
- if ( !parse_oneway(&cmd) )
+ if ( !parse_oneway (ds, &cmd, NULL) )
return CMD_FAILURE;
- /* If /MISSING=INCLUDE is set, then user missing values are ignored */
- if (cmd.incl == ONEWAY_INCLUDE )
- value_is_missing = mv_is_value_system_missing;
- else
- value_is_missing = mv_is_value_missing;
-
/* What statistics were requested */
if ( cmd.sbc_statistics )
{
}
}
- ok = multipass_procedure_with_splits (run_oneway, &cmd);
+ ok = multipass_procedure_with_splits (ds, run_oneway, &cmd);
free (vars);
free_oneway (&cmd);
/* Parser for the variables sub command */
static int
-oneway_custom_variables(struct cmd_oneway *cmd UNUSED)
+oneway_custom_variables(struct dataset *ds, struct cmd_oneway *cmd UNUSED,
+ void *aux UNUSED)
{
+ struct dictionary *dict = dataset_dict (ds);
lex_match('=');
- if ((token != T_ID || dict_lookup_var (default_dict, tokid) == NULL)
+ if ((token != T_ID || dict_lookup_var (dict, tokid) == NULL)
&& token != T_ALL)
return 2;
-
- if (!parse_variables (default_dict, &vars, &n_vars,
+ if (!parse_variables (dict, &vars, &n_vars,
PV_DUPLICATE
| PV_NUMERIC | PV_NO_SCRATCH) )
{
if ( ! lex_match(T_BY))
return 2;
-
- indep_var = parse_variable();
+ indep_var = parse_variable (dict);
if ( !indep_var )
{
return 0;
}
-
return 1;
}
}
- tab_title (t, 0, _("ANOVA"));
+ tab_title (t, _("ANOVA"));
tab_submit (t);
-
-
}
+
/* Show the descriptives table */
static void
show_descriptives(void)
tab_text (t, 9, 1, TAB_CENTER | TAT_TITLE, _("Maximum"));
- tab_title (t, 0, _("Descriptives"));
+ tab_title (t, _("Descriptives"));
row = 2;
const char *s = var_to_string(vars[v]);
- struct group_statistics *const *gs_array = hsh_sort(gp->group_hash);
+ struct group_statistics *const *gs_array =
+ (struct group_statistics *const *) hsh_sort(gp->group_hash);
int count = 0;
tab_text (t, 0, row, TAB_LEFT | TAT_TITLE, s);
tab_text (t, 4, 0, TAB_CENTER | TAT_TITLE, _("Significance"));
- tab_title (t, 0, _("Test of Homogeneity of Variances"));
+ tab_title (t, _("Test of Homogeneity of Variances"));
for ( v=0 ; v < n_vars ; ++v )
{
}
tab_submit (t);
-
-
}
/* Show the contrast coefficients table */
static void
-show_contrast_coeffs(short *bad_contrast)
+show_contrast_coeffs (short *bad_contrast)
{
int n_cols = 2 + ostensible_number_of_groups;
int n_rows = 2 + cmd.sbc_contrast;
tab_vline(t, TAL_2, 2, 0, n_rows - 1);
- tab_title (t, 0, _("Contrast Coefficients"));
+ tab_title (t, _("Contrast Coefficients"));
tab_text (t, 0, 2, TAB_LEFT | TAT_TITLE, _("Contrast"));
tab_vline(t, TAL_2, 3, 0, n_rows - 1);
- tab_title (t, 0, _("Contrast Tests"));
+ tab_title (t, _("Contrast Tests"));
tab_text (t, 2, 0, TAB_CENTER | TAT_TITLE, _("Contrast"));
tab_text (t, 3, 0, TAB_CENTER | TAT_TITLE, _("Value of Contrast"));
static bool
-run_oneway(const struct casefile *cf, void *cmd_)
+run_oneway(const struct ccase *first, const struct casefile *cf,
+ void *cmd_, const struct dataset *ds)
{
struct casereader *r;
struct ccase c;
+ struct casefilter *filter = NULL;
struct cmd_oneway *cmd = (struct cmd_oneway *) cmd_;
+ output_split_file_values (ds, first);
+
global_group_hash = hsh_create(4,
(hsh_compare_func *) compare_values,
(hsh_hash_func *) hash_value,
0,
(void *) indep_var->width );
+
precalc(cmd);
- for(r = casefile_get_reader (cf);
+ filter = casefilter_create ( (cmd->incl != ONEWAY_INCLUDE),
+ vars, n_vars );
+
+ for(r = casefile_get_reader (cf, filter);
casereader_read (r, &c) ;
case_destroy (&c))
{
size_t i;
const double weight =
- dict_get_case_weight(default_dict,&c,&bad_weight_warn);
+ dict_get_case_weight (dataset_dict (ds), &c, &bad_weight_warn);
- const union value *indep_val = case_data (&c, indep_var->fv);
+ const union value *indep_val;
- /* Deal with missing values */
- if ( value_is_missing(&indep_var->miss, indep_val) )
+ if ( casefilter_variable_missing (filter, &c, indep_var))
continue;
- /* Skip the entire case if /MISSING=LISTWISE is set */
- if ( cmd->miss == ONEWAY_LISTWISE )
- {
- for(i = 0; i < n_vars ; ++i)
- {
- const struct variable *v = vars[i];
- const union value *val = case_data (&c, v->fv);
-
- if (value_is_missing(&v->miss, val) )
- break;
- }
- if ( i != n_vars )
- continue;
-
- }
-
+ indep_val = case_data (&c, indep_var->fv);
hsh_insert ( global_group_hash, (void *) indep_val );
hsh_insert ( group_hash, (void *) gs );
}
-
- if (! value_is_missing(&v->miss, val) )
+
+ if (! casefilter_variable_missing (filter, &c, v))
{
struct group_statistics *totals = &gp->ugs;
}
}
+
casereader_destroy (r);
postcalc(cmd);
if ( stat_tables & STAT_HOMO )
- levene(cf, indep_var, n_vars, vars,
- (cmd->miss == ONEWAY_LISTWISE) ? LEV_LISTWISE : LEV_ANALYSIS ,
- value_is_missing);
+ levene (dataset_dict (ds), cf, indep_var, n_vars, vars,
+ filter);
+
+ casefilter_destroy (filter);
ostensible_number_of_groups = hsh_count (global_group_hash);