static void regression_analysis(int n, gmx_bool bXYdy, real* x, int nset, real** val)
{
- real S, chi2, a, b, da, db, r = 0;
- int ok;
+ real S, chi2, a, b, da, db, r = 0;
+ StatisticsStatus ok;
if (bXYdy || (nset == 1))
{
printf("(use option -xydy).\n\n");
if (bXYdy)
{
- if ((ok = lsq_y_ax_b_error(n, x, val[0], val[1], &a, &b, &da, &db, &r, &S)) != estatsOK)
+ if ((ok = lsq_y_ax_b_error(n, x, val[0], val[1], &a, &b, &da, &db, &r, &S))
+ != StatisticsStatus::Ok)
{
- gmx_fatal(FARGS, "Error fitting the data: %s", gmx_stats_message(ok));
+ gmx_stats_message(ok);
}
}
else
{
- if ((ok = lsq_y_ax_b(n, x, val[0], &a, &b, &r, &S)) != estatsOK)
+ if ((ok = lsq_y_ax_b(n, x, val[0], &a, &b, &r, &S)) != StatisticsStatus::Ok)
{
- gmx_fatal(FARGS, "Error fitting the data: %s", gmx_stats_message(ok));
+ gmx_stats_message(ok);
}
}
chi2 = gmx::square((n - 2) * S);