projects
/
pspp
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
lexer: Fix negative 'n0' and 'n1' in lex_source_contains_macro_call().
[pspp]
/
src
/
math
/
linreg.c
diff --git
a/src/math/linreg.c
b/src/math/linreg.c
index b5ccbc599921cad6b9fd1ecf02a5d24ef4abacb2..d417599babd94f9fa9afc9030b457ca8e03bfcb4 100644
(file)
--- a/
src/math/linreg.c
+++ b/
src/math/linreg.c
@@
-275,11
+275,11
@@
linreg_predict (const struct linreg *c, const double *vals, size_t n_vals)
size_t j;
double result;
size_t j;
double result;
- assert (n_vals == c->n_coeffs);
if (vals == NULL || c == NULL)
{
return GSL_NAN;
}
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. */
if (c->coeff == NULL)
{
/* The stupid model: just guess the mean. */
@@
-322,6
+322,7
@@
linreg_set_indep_variable_mean (struct linreg *c, size_t j, double m)
gsl_vector_set (c->indep_means, j, m);
}
gsl_vector_set (c->indep_means, j, m);
}
+#if 0
static void
linreg_fit_qr (const gsl_matrix *cov, struct linreg *l)
{
static void
linreg_fit_qr (const gsl_matrix *cov, struct linreg *l)
{
@@
-415,6
+416,7
@@
linreg_fit_qr (const gsl_matrix *cov, struct linreg *l)
gsl_matrix_free (xtx);
gsl_vector_free (params);
}
gsl_matrix_free (xtx);
gsl_vector_free (params);
}
+#endif
#define REG_LARGE_DATA 1000
#define REG_LARGE_DATA 1000
@@
-443,7
+445,7
@@
linreg_fit (const gsl_matrix *cov, struct linreg *l)
linreg_fit_qr (cov, l);
}
else
linreg_fit_qr (cov, l);
}
else
-#endif
+#endif
{
gsl_matrix *params = gsl_matrix_calloc (cov->size1, cov->size2);
gsl_matrix_memcpy (params, cov);
{
gsl_matrix *params = gsl_matrix_calloc (cov->size1, cov->size2);
gsl_matrix_memcpy (params, cov);
@@
-528,19
+530,19
@@
double linreg_sst (const struct linreg *c)
}
double
}
double
-linreg_dfmodel (
const struct linreg *c)
+linreg_dfmodel (const struct linreg *c)
{
return c->dfm;
}
double
{
return c->dfm;
}
double
-linreg_dferror (
const struct linreg *c)
+linreg_dferror (const struct linreg *c)
{
return c->dfe;
}
double
{
return c->dfe;
}
double
-linreg_dftotal (
const struct linreg *c)
+linreg_dftotal (const struct linreg *c)
{
return c->dft;
}
{
return c->dft;
}