wk = gsl_multifit_linear_alloc (design->size1, design->size2);
rc = gsl_multifit_linear (design, Y, param_estimates,
cache->cov, &(cache->sse), wk);
- for (i = 1; i < cache->n_coeffs; i++)
+ for (i = 0; i < cache->n_coeffs; i++)
{
- cache->coeff[i]->estimate = gsl_vector_get (param_estimates, i);
+ cache->coeff[i]->estimate = gsl_vector_get (param_estimates, i + 1);
}
cache->intercept = gsl_vector_get (param_estimates, 0);
if (rc == GSL_SUCCESS)