bool intercept;
double alpha;
+
+ bool dump_coding;
};
struct glm_workspace
glm.intercept = true;
glm.wv = dict_get_weight (glm.dict);
glm.alpha = 0.05;
+ glm.dump_coding = false;
if (!parse_variables_const (lexer, glm.dict,
&glm.dep_vars, &glm.n_dep_vars,
if (glm.n_interactions > 0)
design = true;
}
+ else if (lex_match_id (lexer, "SHOWCODES"))
+ /* Undocumented debug option */
+ {
+ lex_match (lexer, T_EQUALS);
+
+ glm.dump_coding = true;
+ }
else
{
lex_error (lexer, NULL);
}
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);
+ submodel_cov = gsl_matrix_alloc (cm->size1 - n_dropped_submodel, cm->size2 - n_dropped_submodel);
fill_submatrix (cm, model_cov, model_dropped);
fill_submatrix (cm, submodel_cov, submodel_dropped);
}
casereader_destroy (reader);
- for (reader = input;
+ if (cmd->dump_coding)
+ reader = casereader_clone (input);
+ else
+ reader = input;
+
+ for (;
(c = casereader_read (reader)) != NULL; case_unref (c))
{
double weight = dict_get_case_weight (dict, c, &warn_bad_weight);
}
casereader_destroy (reader);
+
+ if (cmd->dump_coding)
+ {
+ struct tab_table *t =
+ covariance_dump_enc_header (cov,
+ 1 + casereader_count_cases (input));
+ for (reader = input;
+ (c = casereader_read (reader)) != NULL; case_unref (c))
+ {
+ covariance_dump_enc (cov, c, t);
+ }
+ casereader_destroy (reader);
+ tab_submit (t);
+ }
+
{
gsl_matrix *cm = covariance_calculate_unnormalized (cov);