/* The new variables which will be introduced by /SAVE */
const struct variable **predvars;
const struct variable **residvars;
/* The new variables which will be introduced by /SAVE */
const struct variable **predvars;
const struct variable **residvars;
struct regression_workspace *ws,
struct casereader *input);
struct regression_workspace *ws,
struct casereader *input);
struct regression_workspace workspace;
struct regression regression;
const struct dictionary *dict = dataset_dict (ds);
bool save;
struct regression_workspace workspace;
struct regression regression;
const struct dictionary *dict = dataset_dict (ds);
bool save;
- workspace.psw = xrealloc (workspace.psw, ++n_splits * sizeof (*workspace.psw));
-
- run_regression (®ression, &workspace.psw[n_splits - 1],
+ run_regression (®ression,
add_transformation (ds, save_trans_func, save_trans_free, save_trans_data);
}
add_transformation (ds, save_trans_func, save_trans_free, save_trans_data);
}
struct regression_workspace *ws,
struct casereader *input)
{
size_t i;
struct regression_workspace *ws,
struct casereader *input)
{
size_t i;
const struct variable **vars = xnmalloc (cmd->n_vars, sizeof (*vars));
const struct variable *dep_var = cmd->dep_vars[k];
int n_indep = identify_indep_vars (cmd, vars, dep_var);
gsl_matrix *this_cm = gsl_matrix_alloc (n_indep + 1, n_indep + 1);
double n_data = fill_covariance (this_cm, cov, vars, n_indep,
dep_var, all_vars, n_all_vars, means);
const struct variable **vars = xnmalloc (cmd->n_vars, sizeof (*vars));
const struct variable *dep_var = cmd->dep_vars[k];
int n_indep = identify_indep_vars (cmd, vars, dep_var);
gsl_matrix *this_cm = gsl_matrix_alloc (n_indep + 1, n_indep + 1);
double n_data = fill_covariance (this_cm, cov, vars, n_indep,
dep_var, all_vars, n_all_vars, means);
- psw->models[k] = linreg_alloc (dep_var, vars, n_data, n_indep);
- psw->models[k]->depvar = dep_var;
+ models[k] = linreg_alloc (dep_var, vars, n_data, n_indep);
+ models[k]->depvar = dep_var;
/*
For large data sets, use QR decomposition.
*/
if (n_data > sqrt (n_indep) && n_data > REG_LARGE_DATA)
{
/*
For large data sets, use QR decomposition.
*/
if (n_data > sqrt (n_indep) && n_data > REG_LARGE_DATA)
{
- subcommand_statistics (cmd, psw->models[k], this_cm, dep_var);
+ subcommand_statistics (cmd, models[k], this_cm, dep_var);
- double obs = case_data (c, psw->models[k]->depvar)->f;
- double res = linreg_residual (psw->models[k], obs, vals, n_indep);
+ double obs = case_data (c, models[k]->depvar)->f;
+ double res = linreg_residual (models[k], obs, vals, n_indep);