projects
/
pspp-builds.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (from parent 1:
07fac5f
)
Fixed centering prior to sweep operation
author
Jason Stover
<jhs@math.gcsu.edu>
Wed, 4 Jan 2006 18:02:14 +0000
(18:02 +0000)
committer
Jason Stover
<jhs@math.gcsu.edu>
Wed, 4 Jan 2006 18:02:14 +0000
(18:02 +0000)
src/regression.q
patch
|
blob
|
history
diff --git
a/src/regression.q
b/src/regression.q
index eb8ff34f76ad353740f843c20f129ce7e807b9ce..60f3c39f121075fee4d8e514b7197d9aea779eff 100644
(file)
--- a/
src/regression.q
+++ b/
src/regression.q
@@
-822,18
+822,22
@@
run_regression (const struct casefile *cf, void *cmd_ UNUSED)
X =
design_matrix_create (n_indep, (const struct variable **) indep_vars,
n_data);
X =
design_matrix_create (n_indep, (const struct variable **) indep_vars,
n_data);
+ for (i = 0; i < X->m->size2; i++)
+ {
+ lopts.get_indep_mean_std[i] = 1;
+ }
lcache = pspp_linreg_cache_alloc (X->m->size1, X->m->size2);
lcache->indep_means = gsl_vector_alloc (X->m->size2);
lcache->indep_std = gsl_vector_alloc (X->m->size2);
lcache->depvar = (const struct variable *) depvar;
/*
lcache = pspp_linreg_cache_alloc (X->m->size1, X->m->size2);
lcache->indep_means = gsl_vector_alloc (X->m->size2);
lcache->indep_std = gsl_vector_alloc (X->m->size2);
lcache->depvar = (const struct variable *) depvar;
/*
-
For large data sets, use QR decomposition.
-
*/
+ For large data sets, use QR decomposition.
+ */
if (n_data > sqrt (n_indep) && n_data > REG_LARGE_DATA)
{
lcache->method = PSPP_LINREG_SVD;
}
if (n_data > sqrt (n_indep) && n_data > REG_LARGE_DATA)
{
lcache->method = PSPP_LINREG_SVD;
}
-
+
/*
The second pass creates the design matrix.
*/
/*
The second pass creates the design matrix.
*/
@@
-869,7
+873,6
@@
run_regression (const struct casefile *cf, void *cmd_ UNUSED)
{
design_matrix_set_numeric (X, row, v, val);
}
{
design_matrix_set_numeric (X, row, v, val);
}
- lopts.get_indep_mean_std[i] = 1;
}
}
val = case_data (&c, depvar->fv);
}
}
val = case_data (&c, depvar->fv);