X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Flanguage%2Fstats%2Fbinomial.c;h=d9209f1bf0499dbcf8e78f8bb9ba1b6e27861683;hb=e051d2cccc85e4ae2c46bddc7c674676ba7a5515;hp=f82647a405b176c5fa95f5ab90c4e864524e6ed4;hpb=9f087e7aa4cdff1d5d46d5e188c0017a9d2d0029;p=pspp-builds.git diff --git a/src/language/stats/binomial.c b/src/language/stats/binomial.c index f82647a4..d9209f1b 100644 --- a/src/language/stats/binomial.c +++ b/src/language/stats/binomial.c @@ -41,7 +41,6 @@ #include #include -#include #include @@ -79,7 +78,7 @@ calculate_binomial_internal (double n1, double n2, double p) /* SPSS Statistical Algorithms has completely different and WRONG advice here. */ - double sig1tailed = gslextras_cdf_binomial_P (n1, n1 + n2, p); + double sig1tailed = gsl_cdf_binomial_P (n1, p, n1 + n2); if ( p == 0.5 ) return sig1tailed > 0.5 ? 1.0 :sig1tailed * 2.0; @@ -185,7 +184,16 @@ binomial_execute (const struct dataset *ds, for (v = 0 ; v < ost->n_vars; ++v) { double n_total, sig; + struct string catstr1; + struct string catstr2; const struct variable *var = ost->vars[v]; + + ds_init_empty (&catstr1); + ds_init_empty (&catstr2); + + var_append_value_name (var, cat1[v].value, &catstr1); + var_append_value_name (var, cat2[v].value, &catstr2); + tab_hline (table, TAL_1, 0, tab_nc (table) -1, 1 + v * 3); /* Titles */ @@ -198,10 +206,8 @@ binomial_execute (const struct dataset *ds, tab_float (table, 5, 1 + v * 3, TAB_NONE, bst->p, 8, 3); /* Category labels */ - tab_text (table, 2, 1 + v * 3, TAB_NONE, - var_get_value_name (var, cat1[v].value)); - tab_text (table, 2, 2 + v * 3, TAB_NONE, - var_get_value_name (var, cat2[v].value)); + tab_text (table, 2, 1 + v * 3, TAB_NONE, ds_cstr (&catstr1)); + tab_text (table, 2, 2 + v * 3, TAB_NONE, ds_cstr (&catstr2)); /* Observed N */ tab_float (table, 3, 1 + v * 3, TAB_NONE, cat1[v].count, 8, 0); @@ -221,6 +227,9 @@ binomial_execute (const struct dataset *ds, /* Significance */ sig = calculate_binomial (cat1[v].count, cat2[v].count, bst->p); tab_float (table, 6, 1 + v * 3, TAB_NONE, sig, 8, 3); + + ds_destroy (&catstr1); + ds_destroy (&catstr2); } tab_text (table, 2, 0, TAB_CENTER, _("Category"));