size_t j;
double result;
- assert (n_vals == c->n_coeffs);
if (vals == NULL || c == NULL)
{
return GSL_NAN;
}
+ assert (n_vals == c->n_coeffs);
if (c->coeff == NULL)
{
/* The stupid model: just guess the mean. */
gsl_vector_set (c->indep_means, j, m);
}
+#if 0
static void
linreg_fit_qr (const gsl_matrix *cov, struct linreg *l)
{
gsl_matrix_free (xtx);
gsl_vector_free (params);
}
+#endif
#define REG_LARGE_DATA 1000
linreg_fit_qr (cov, l);
}
else
-#endif
+#endif
{
gsl_matrix *params = gsl_matrix_calloc (cov->size1, cov->size2);
gsl_matrix_memcpy (params, cov);
}
double
-linreg_dfmodel ( const struct linreg *c)
+linreg_dfmodel (const struct linreg *c)
{
return c->dfm;
}
double
-linreg_dferror ( const struct linreg *c)
+linreg_dferror (const struct linreg *c)
{
return c->dfe;
}
double
-linreg_dftotal ( const struct linreg *c)
+linreg_dftotal (const struct linreg *c)
{
return c->dft;
}