projects
/
pspp
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
render: Add comment.
[pspp]
/
src
/
language
/
stats
/
logistic.c
diff --git
a/src/language/stats/logistic.c
b/src/language/stats/logistic.c
index 6d2366518b301606e68b720c218dc769a4148e69..2cf440c31f8829f34ee0416a26feea1ea12adf0f 100644
(file)
--- a/
src/language/stats/logistic.c
+++ b/
src/language/stats/logistic.c
@@
-242,7
+242,7
@@
predictor_value (const struct ccase *c,
/*
/*
- Return the probability beta_hat (that is the estimator logit(y)
)
+ Return the probability beta_hat (that is the estimator logit(y))
corresponding to the coefficient estimator for case C
*/
static double
corresponding to the coefficient estimator for case C
*/
static double
@@
-321,7
+321,7
@@
hessian (const struct lr_spec *cmd,
}
casereader_destroy (reader);
}
casereader_destroy (reader);
- if (
max_w < cmd->min_epsilon)
+ if (max_w < cmd->min_epsilon)
{
*converged = true;
msg (MN, _("All predicted values are either 1 or 0"));
{
*converged = true;
msg (MN, _("All predicted values are either 1 or 0"));
@@
-418,7
+418,7
@@
frq_update (const void *aux1 UNUSED, void *aux2 UNUSED,
}
static void
}
static void
-frq_destroy (const void *aux1 UNUSED, void *aux2 UNUSED, void *user_data
UNUSED
)
+frq_destroy (const void *aux1 UNUSED, void *aux2 UNUSED, void *user_data)
{
free (user_data);
}
{
free (user_data);
}
@@
-509,7
+509,7
@@
initial_pass (const struct lr_spec *cmd, struct lr_result *res, struct casereade
}
else if (!v1set)
{
}
else if (!v1set)
{
- if (
!value_equal (&res->y0, depval, width))
+ if (!value_equal (&res->y0, depval, width))
{
value_clone (&res->y1, depval, width);
v1set = true;
{
value_clone (&res->y1, depval, width);
v1set = true;
@@
-520,7
+520,7
@@
initial_pass (const struct lr_spec *cmd, struct lr_result *res, struct casereade
if (! value_equal (&res->y0, depval, width)
&&
! value_equal (&res->y1, depval, width)
if (! value_equal (&res->y0, depval, width)
&&
! value_equal (&res->y1, depval, width)
-
)
+ )
{
msg (ME, _("Dependent variable's values are not dichotomous."));
case_unref (c);
{
msg (ME, _("Dependent variable's values are not dichotomous."));
case_unref (c);
@@
-640,7
+640,7
@@
run_lr (const struct lr_spec *cmd, struct casereader *input,
work.hessian = gsl_matrix_calloc (work.beta_hat->size, work.beta_hat->size);
/* Start the Newton Raphson iteration process... */
work.hessian = gsl_matrix_calloc (work.beta_hat->size, work.beta_hat->size);
/* Start the Newton Raphson iteration process... */
- for(
i = 0 ; i < cmd->max_iter ; ++i)
+ for(i = 0 ; i < cmd->max_iter ; ++i)
{
double min, max;
gsl_vector *v ;
{
double min, max;
gsl_vector *v ;
@@
-669,7
+669,7
@@
run_lr (const struct lr_spec *cmd, struct casereader *input,
gsl_vector_minmax (delta, &min, &max);
gsl_vector_minmax (delta, &min, &max);
- if (
fabs (min) < cmd->bcon && fabs (max) < cmd->bcon)
+ if (fabs (min) < cmd->bcon && fabs (max) < cmd->bcon)
{
msg (MN, _("Estimation terminated at iteration number %d because parameter estimates changed by less than %g"),
i + 1, cmd->bcon);
{
msg (MN, _("Estimation terminated at iteration number %d because parameter estimates changed by less than %g"),
i + 1, cmd->bcon);
@@
-697,8
+697,8
@@
run_lr (const struct lr_spec *cmd, struct casereader *input,
- if (
! converged)
- msg (MW, _("Estimation terminated at iteration number %d because maximum iterations has been reached"), i
);
+ if (! converged)
+ msg (MW, _("Estimation terminated at iteration number %d because maximum iterations has been reached"), i);
output_model_summary (&work, initial_log_likelihood, log_likelihood);
output_model_summary (&work, initial_log_likelihood, log_likelihood);
@@
-739,8
+739,6
@@
lookup_variable (const struct hmap *map, const struct variable *var, unsigned in
{
if (vn->var == var)
break;
{
if (vn->var == var)
break;
-
- fprintf (stderr, "Warning: Hash table collision\n");
}
return vn;
}
return vn;
@@
-885,7
+883,7
@@
cmd_logistic (struct lexer *lexer, struct dataset *ds)
}
lr.confidence = lex_number (lexer);
lex_get (lexer);
}
lr.confidence = lex_number (lexer);
lex_get (lexer);
- if (
! lex_force_match (lexer, T_RPAREN))
+ if (! lex_force_match (lexer, T_RPAREN))
{
lex_error (lexer, NULL);
goto error;
{
lex_error (lexer, NULL);
goto error;
@@
-919,7
+917,7
@@
cmd_logistic (struct lexer *lexer, struct dataset *ds)
}
lr.bcon = lex_number (lexer);
lex_get (lexer);
}
lr.bcon = lex_number (lexer);
lex_get (lexer);
- if (
! lex_force_match (lexer, T_RPAREN))
+ if (! lex_force_match (lexer, T_RPAREN))
{
lex_error (lexer, NULL);
goto error;
{
lex_error (lexer, NULL);
goto error;
@@
-937,7
+935,7
@@
cmd_logistic (struct lexer *lexer, struct dataset *ds)
}
lr.max_iter = lex_integer (lexer);
lex_get (lexer);
}
lr.max_iter = lex_integer (lexer);
lex_get (lexer);
- if (
! lex_force_match (lexer, T_RPAREN))
+ if (! lex_force_match (lexer, T_RPAREN))
{
lex_error (lexer, NULL);
goto error;
{
lex_error (lexer, NULL);
goto error;
@@
-955,7
+953,7
@@
cmd_logistic (struct lexer *lexer, struct dataset *ds)
}
lr.lcon = lex_number (lexer);
lex_get (lexer);
}
lr.lcon = lex_number (lexer);
lex_get (lexer);
- if (
! lex_force_match (lexer, T_RPAREN))
+ if (! lex_force_match (lexer, T_RPAREN))
{
lex_error (lexer, NULL);
goto error;
{
lex_error (lexer, NULL);
goto error;
@@
-973,7
+971,7
@@
cmd_logistic (struct lexer *lexer, struct dataset *ds)
}
lr.min_epsilon = lex_number (lexer);
lex_get (lexer);
}
lr.min_epsilon = lex_number (lexer);
lex_get (lexer);
- if (
! lex_force_match (lexer, T_RPAREN))
+ if (! lex_force_match (lexer, T_RPAREN))
{
lex_error (lexer, NULL);
goto error;
{
lex_error (lexer, NULL);
goto error;
@@
-997,7
+995,7
@@
cmd_logistic (struct lexer *lexer, struct dataset *ds)
goto error;
}
lex_get (lexer);
goto error;
}
lex_get (lexer);
- if (
! lex_force_match (lexer, T_RPAREN))
+ if (! lex_force_match (lexer, T_RPAREN))
{
lex_error (lexer, NULL);
goto error;
{
lex_error (lexer, NULL);
goto error;
@@
-1314,7
+1312,7
@@
output_model_summary (const struct lr_result *res,
double entries[] = {
-2 * log_likelihood,
cox,
double entries[] = {
-2 * log_likelihood,
cox,
- cox / (
1.0 - exp(initial_log_likelihood * (2 / res->cc)))
+ cox / (1.0 - exp(initial_log_likelihood * (2 / res->cc)))
};
for (size_t i = 0; i < sizeof entries / sizeof *entries; i++)
pivot_table_put2 (table, i, 0, pivot_value_new_number (entries[i]));
};
for (size_t i = 0; i < sizeof entries / sizeof *entries; i++)
pivot_table_put2 (table, i, 0, pivot_value_new_number (entries[i]));