+ }
+}
+
+static void
+get_ssq (struct covariance *cov, gsl_vector *ssq, const struct glm_spec *cmd)
+{
+ gsl_matrix *cm = covariance_calculate_unnormalized (cov);
+ size_t i;
+ size_t k;
+ size_t *model_dropped = xcalloc (covariance_dim (cov), sizeof (*model_dropped));
+ size_t *submodel_dropped = xcalloc (covariance_dim (cov), sizeof (*submodel_dropped));
+ const struct categoricals *cats = covariance_get_categoricals (cov);
+
+ for (k = 0; k < cmd->n_interactions; k++)
+ {
+ gsl_matrix *model_cov = NULL;
+ gsl_matrix *submodel_cov = NULL;
+ size_t n_dropped_model = 0;
+ size_t n_dropped_submodel = 0;
+ for (i = cmd->n_dep_vars; i < covariance_dim (cov); i++)
+ {
+ const struct interaction * x =
+ categoricals_get_interaction_by_subscript (cats, i - cmd->n_dep_vars);
+ if (is_subset (cmd->interactions [k], x))
+ {
+ assert (n_dropped_model < covariance_dim (cov));
+ model_dropped[n_dropped_model++] = i;
+ }
+ if (drop_from_submodel (cmd->interactions [k], x))
+ {
+ assert (n_dropped_submodel < covariance_dim (cov));
+ submodel_dropped[n_dropped_submodel++] = i;
+ }
+ }
+ model_cov =
+ gsl_matrix_alloc (cm->size1 - n_dropped_model, cm->size2 - n_dropped_model);
+ gsl_matrix_set (model_cov, 0, 0, gsl_matrix_get (cm, 0, 0));
+ submodel_cov =
+ gsl_matrix_calloc (cm->size1 - n_dropped_submodel, cm->size2 - n_dropped_submodel);
+ fill_submatrix (cm, model_cov, model_dropped, n_dropped_model);
+ fill_submatrix (cm, submodel_cov, submodel_dropped, n_dropped_submodel);
+
+ reg_sweep (model_cov, 0);
+ reg_sweep (submodel_cov, 0);