X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Foneway.q;h=565ba61b7f64e3d565b560061e9d2d33b6c10b83;hb=4fdeb2145d081ff1b84e3f6c99f9d1c048c0d64a;hp=bcc75c7e8c23fccbcbe64044080f2a1306484f76;hpb=f500c9c2989d63465b9a93fe6f7e1600207681af;p=pspp diff --git a/src/oneway.q b/src/oneway.q index bcc75c7e8c..565ba61b7f 100644 --- a/src/oneway.q +++ b/src/oneway.q @@ -15,8 +15,8 @@ General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software -Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA -02111-1307, USA. */ +Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +02110-1301, USA. */ #include #include @@ -27,6 +27,7 @@ Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA #include "alloc.h" #include "str.h" #include "case.h" +#include "dictionary.h" #include "command.h" #include "lexer.h" #include "error.h" @@ -42,10 +43,11 @@ Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA #include "group_proc.h" #include "group.h" #include "levene.h" +/* (headers) */ /* (specification) "ONEWAY" (oneway_): - *variables=custom; + *^variables=custom; +missing=miss:!analysis/listwise, incl:include/!exclude; contrast= double list; @@ -138,6 +140,8 @@ cmd_oneway(void) multipass_procedure_with_splits (run_oneway, &cmd); + free (vars); + free_oneway (&cmd); return CMD_SUCCESS; } @@ -195,7 +199,7 @@ output_oneway(void) /* Clean up */ for (i = 0 ; i < n_vars ; ++i ) { - struct hsh_table *group_hash = vars[i]->p.grp_data.group_hash; + struct hsh_table *group_hash = group_proc_get (vars[i])->group_hash; hsh_destroy(group_hash); } @@ -281,8 +285,8 @@ show_anova_table(void) for ( i=0 ; i < n_vars ; ++i ) { - struct group_statistics *totals = &vars[i]->p.grp_data.ugs; - struct hsh_table *group_hash = vars[i]->p.grp_data.group_hash; + struct group_statistics *totals = &group_proc_get (vars[i])->ugs; + struct hsh_table *group_hash = group_proc_get (vars[i])->group_hash; struct hsh_iterator g; struct group_statistics *gs; double ssa=0; @@ -308,12 +312,13 @@ show_anova_table(void) tab_hline(t, TAL_1, 0, n_cols - 1 , i * 3 + 1); { + struct group_proc *gp = group_proc_get (vars[i]); const double sst = totals->ssq - ( totals->sum * totals->sum) / totals->n ; - const double df1 = vars[i]->p.grp_data.n_groups - 1; - const double df2 = totals->n - vars[i]->p.grp_data.n_groups ; + const double df1 = gp->n_groups - 1; + const double df2 = totals->n - gp->n_groups ; const double msa = ssa / df1; - vars[i]->p.grp_data.mse = (sst - ssa) / df2; + gp->mse = (sst - ssa) / df2; /* Sums of Squares */ @@ -329,11 +334,11 @@ show_anova_table(void) /* Mean Squares */ tab_float (t, 4, i * 3 + 1, TAB_RIGHT, msa, 8, 3); - tab_float (t, 4, i * 3 + 2, TAB_RIGHT, vars[i]->p.grp_data.mse, 8, 3); + tab_float (t, 4, i * 3 + 2, TAB_RIGHT, gp->mse, 8, 3); { - const double F = msa/vars[i]->p.grp_data.mse ; + const double F = msa/gp->mse ; /* The F value */ tab_float (t, 5, i * 3 + 1, 0, F, 8, 3); @@ -371,7 +376,7 @@ show_descriptives(void) for ( v = 0 ; v < n_vars ; ++v ) - n_rows += vars[v]->p.grp_data.n_groups + 1; + n_rows += group_proc_get (vars[v])->n_groups + 1; t = tab_create (n_cols,n_rows,0); tab_headers (t, 2, 0, 2, 0); @@ -414,16 +419,17 @@ show_descriptives(void) { double T; double std_error; - + + struct group_proc *gp = group_proc_get (vars[v]); struct hsh_iterator g; struct group_statistics *gs; - struct group_statistics *totals = &vars[v]->p.grp_data.ugs; + struct group_statistics *totals = &gp->ugs; int count = 0 ; const char *s = var_to_string(vars[v]); - struct hsh_table *group_hash = vars[v]->p.grp_data.group_hash; + struct hsh_table *group_hash = gp->group_hash; tab_text (t, 0, row, TAB_LEFT | TAT_TITLE, s); @@ -435,18 +441,8 @@ show_descriptives(void) gs != 0; gs = hsh_next(group_hash,&g)) { - const char *s = val_labs_find(indep_var->val_labs, gs->id ); - - if ( s ) - tab_text (t, 1, row + count, - TAB_LEFT | TAT_TITLE ,s); - else if ( indep_var->width != 0 ) - tab_text (t, 1, row + count, - TAB_LEFT | TAT_TITLE, gs->id.s); - else - tab_text (t, 1, row + count, - TAB_LEFT | TAT_TITLE | TAT_PRINTF, "%g", gs->id.f); - + tab_text (t, 1, row + count, + TAB_LEFT | TAT_TITLE ,value_to_string(&gs->id,indep_var)); /* Now fill in the numbers ... */ @@ -506,7 +502,7 @@ show_descriptives(void) tab_float(t, 8, row + count, 0, totals->minimum, 8, 2); tab_float(t, 9, row + count, 0, totals->maximum, 8, 2); - row += vars[v]->p.grp_data.n_groups + 1; + row += gp->n_groups + 1; } @@ -554,15 +550,16 @@ show_homogeneity(void) { double F; const struct variable *var = vars[v]; + const struct group_proc *gp = group_proc_get (vars[v]); const char *s = var_to_string(var); - const struct group_statistics *totals = &var->p.grp_data.ugs; + const struct group_statistics *totals = &gp->ugs; - const double df1 = var->p.grp_data.n_groups - 1; - const double df2 = totals->n - var->p.grp_data.n_groups ; + const double df1 = gp->n_groups - 1; + const double df2 = totals->n - gp->n_groups ; tab_text (t, 0, v + 1, TAB_LEFT | TAT_TITLE, s); - F = var->p.grp_data.levene; + F = gp->levene; tab_float (t, 1, v + 1, TAB_RIGHT, F, 8,3); tab_float (t, 2, v + 1, TAB_RIGHT, df1 ,8,0); tab_float (t, 3, v + 1, TAB_RIGHT, df2 ,8,0); @@ -716,7 +713,7 @@ show_contrast_tests(short *bad_contrast) int ci; double contrast_value = 0.0; double coef_msq = 0.0; - struct group_proc *grp_data = &vars[v]->p.grp_data ; + struct group_proc *grp_data = group_proc_get (vars[v]); struct hsh_table *group_hash = grp_data->group_hash; struct hsh_iterator g; struct group_statistics *gs; @@ -796,7 +793,7 @@ show_contrast_tests(short *bad_contrast) cmd.sbc_contrast, TAB_RIGHT, contrast_value, 8,2); - std_error_contrast = sqrt(vars[v]->p.grp_data.mse * coef_msq); + std_error_contrast = sqrt(grp_data->mse * coef_msq); /* Std. Error */ tab_float (t, 4, (v * lines_per_variable) + i + 1, @@ -882,13 +879,14 @@ precalc ( struct cmd_oneway *cmd UNUSED ) for(i=0; i< n_vars ; ++i) { - struct group_statistics *totals = &vars[i]->p.grp_data.ugs; + struct group_proc *gp = group_proc_get (vars[i]); + struct group_statistics *totals = &gp->ugs; /* Create a hash for each of the dependent variables. The hash contains a group_statistics structure, and is keyed by value of the independent variable */ - vars[i]->p.grp_data.group_hash = + gp->group_hash = hsh_create(4, (hsh_compare_func *) compare_group, (hsh_hash_func *) hash_group, @@ -961,7 +959,8 @@ run_oneway(const struct casefile *cf, void *cmd_) const union value *val = case_data (&c, v->fv); - struct hsh_table *group_hash = vars[i]->p.grp_data.group_hash; + struct group_proc *gp = group_proc_get (vars[i]); + struct hsh_table *group_hash = gp->group_hash; struct group_statistics *gs; @@ -985,7 +984,7 @@ run_oneway(const struct casefile *cf, void *cmd_) if (! value_is_missing(val,v) ) { - struct group_statistics *totals = &vars[i]->p.grp_data.ugs; + struct group_statistics *totals = &gp->ugs; totals->n+=weight; totals->sum+=weight * val->f; @@ -1008,7 +1007,7 @@ run_oneway(const struct casefile *cf, void *cmd_) gs->maximum = val->f * weight; } - vars[i]->p.grp_data.n_groups = hsh_count ( group_hash ); + gp->n_groups = hsh_count ( group_hash ); } } @@ -1040,8 +1039,9 @@ postcalc ( struct cmd_oneway *cmd UNUSED ) for(i = 0; i < n_vars ; ++i) { - struct hsh_table *group_hash = vars[i]->p.grp_data.group_hash; - struct group_statistics *totals = &vars[i]->p.grp_data.ugs; + struct group_proc *gp = group_proc_get (vars[i]); + struct hsh_table *group_hash = gp->group_hash; + struct group_statistics *totals = &gp->ugs; struct hsh_iterator g; struct group_statistics *gs;