- const struct variable *v = indep_vars[i];
- const union value *val = case_data (&c, v);
- if (var_is_alpha (v))
- design_matrix_set_categorical (X, row, v, val);
- else
- design_matrix_set_numeric (X, row, v, val);
+ const struct variable *v = all_vars[i];
+ const union value *val_v = case_data (&c, v);
+ for (j = i; j < n_all_vars; j++)
+ {
+ const struct variable *w = all_vars[j];
+ const union value *val_w = case_data (&c, w);
+ covariance_pass_two (X, *mom[i]->mean, *mom[j]->mean,
+ (double) 1 / n_data, (double) n_data,
+ v, w, val_v, val_w);
+ }