}
return n_indep_vars;
}
+
static double
fill_covariance (gsl_matrix *cov, struct covariance *all_cov,
const struct variable **vars,
size_t dep_subscript;
size_t *rows;
const gsl_matrix *ssizes;
- gsl_matrix *cm;
const gsl_matrix *mean_matrix;
const gsl_matrix *ssize_matrix;
double result = 0.0;
- cm = covariance_calculate_unnormalized (all_cov);
+ gsl_matrix *cm = covariance_calculate_unnormalized (all_cov);
+
+ if ( cm == NULL)
+ return 0;
+
rows = xnmalloc (cov->size1 - 1, sizeof (*rows));
for (i = 0; i < n_all_vars; i++)
dict_get_weight (dict), MV_ANY);
reader = casereader_clone (input);
- reader = casereader_create_filter_missing (reader, v_variables, n_variables,
+ reader = casereader_create_filter_missing (reader, all_vars, n_all_vars,
MV_ANY, NULL, NULL);
+
+
for (; (c = casereader_read (reader)) != NULL; case_unref (c))
{
covariance_accumulate (cov, c);