+static void
+subcommand_export (int export, pspp_linreg_cache *c)
+{
+ FILE *fp;
+ size_t i;
+ struct pspp_linreg_coeff coeff;
+
+ if (export)
+ {
+ assert (c != NULL);
+ assert (model_file != NULL);
+ assert (fp != NULL);
+ fp = fopen (handle_get_filename (model_file), "w");
+ fprintf (fp, "#include <string.h>\n\n");
+ fprintf (fp, "/*\n Estimate the mean of Y, the dependent variable for\n");
+ fprintf (fp, " the linear model of the form \n\n");
+ fprintf (fp, " Y = b0 + b1 * X1 + b2 * X2 + ... + bk * X2 + error\n\n");
+ fprintf (fp, " where X1, ..., Xk are the independent variables\n");
+ fprintf (fp, " whose values are stored in var_vals and whose names, \n");
+ fprintf (fp, " as known by PSPP, are stored in var_names. The estimated \n");
+ fprintf (fp, " regression coefficients (i.e., the estimates of b0,...,bk) \n");
+ fprintf (fp, " are stored in model_coeffs.\n*/\n");
+ fprintf (fp, "double\npspp_reg_estimate (const double *var_vals, const char *var_names[])\n{\n\tchar *model_depvars[%d] = {", c->n_indeps);
+ for (i = 1; i < c->n_indeps; i++)
+ {
+ coeff = c->coeff[i];
+ fprintf (fp, "\"%s\",\n\t\t", coeff.v->name);
+ }
+ coeff = c->coeff[i];
+ fprintf (fp, "\"%s\"};\n\t", coeff.v->name);
+ fprintf (fp, "double model_coeffs[%d] = {", c->n_indeps);
+ for (i = 1; i < c->n_indeps; i++)
+ {
+ coeff = c->coeff[i];
+ fprintf (fp, "%.15e,\n\t\t", coeff.estimate);
+ }
+ coeff = c->coeff[i];
+ fprintf (fp, "%.15e};\n\t", coeff.estimate);
+ coeff = c->coeff[0];
+ fprintf (fp, "double estimate = %.15e;\n\t", coeff.estimate);
+ fprintf (fp, "int i;\n\tint j;\n\n\t");
+ fprintf (fp, "for (i = 0; i < %d; i++)\n\t", c->n_indeps);
+ fprintf (fp, "{\n\t\tfor (j = 0; j < %d; j++)\n\t\t", c->n_indeps);
+ fprintf (fp, "{\n\t\t\tif (strcmp (var_names[i], model_depvars[j]) == 0)\n");
+ fprintf (fp, "\t\t\t{\n\t\t\t\testimate += var_vals[i] * model_coeffs[j];\n");
+ fprintf (fp, "\t\t\t}\n\t\t}\n\t}\n\treturn estimate;\n}\n");
+ fclose (fp);
+ }
+}
+static int
+regression_custom_export (struct cmd_regression *cmd)
+{
+ /* 0 on failure, 1 on success, 2 on failure that should result in syntax error */
+ if (!lex_force_match ('('))
+ return 0;
+
+ if (lex_match ('*'))
+ model_file = NULL;
+ else
+ {
+ model_file = fh_parse ();
+ if (model_file == NULL)
+ return 0;
+ }
+
+ if (!lex_force_match (')'))
+ return 0;
+
+ return 1;
+}