Springer. 1998. ISBN 0-387-98542-5.
*/
-#include "linreg.h"
-#include "coefficient.h"
-#include "sweep.h"
+#include <math/linreg/linreg.h>
+#include <math/linreg/coefficient.h>
#include <gsl/gsl_errno.h>
+#include <linreg/sweep.h>
/*
Get the mean and standard deviation of a vector
of doubles via a form of the Kalman filter as
}
else
{
+ gsl_multifit_linear_workspace *wk ;
/*
Use QR decomposition via GSL.
*/
gsl_matrix_set (design, j, i + 1, tmp);
}
}
- gsl_multifit_linear_workspace *wk =
- gsl_multifit_linear_alloc (design->size1, design->size2);
+
+ wk = gsl_multifit_linear_alloc (design->size1, design->size2);
rc = gsl_multifit_linear (design, Y, param_estimates,
cache->cov, &(cache->sse), wk);
for (i = 0; i < cache->n_coeffs; i++)