static gboolean
analysis_tool_chi_squared_engine_run (data_analysis_output_t *dao,
analysis_tools_data_chi_squared_t *info)
{
GnmExpr const *expr_check;
GnmExpr const *expr_region;
GnmExpr const *expr_statistic;
GnmExpr const *expr_row_ones;
GnmExpr const *expr_col_ones;
GnmExpr const *expr_row;
GnmExpr const *expr_column;
GnmExpr const *expr_expect;
GnmFunc *fd_mmult = analysis_tool_get_function ("MMULT", dao);
GnmFunc *fd_row = analysis_tool_get_function ("ROW", dao);
GnmFunc *fd_column = analysis_tool_get_function ("COLUMN", dao);
GnmFunc *fd_transpose = analysis_tool_get_function ("TRANSPOSE", dao);
GnmFunc *fd_sum = analysis_tool_get_function ("SUM", dao);
GnmFunc *fd_min = analysis_tool_get_function ("MIN", dao);
GnmFunc *fd_offset = analysis_tool_get_function ("OFFSET", dao);
GnmFunc *fd_chiinv = analysis_tool_get_function ("CHIINV", dao);
GnmFunc *fd_chidist = analysis_tool_get_function ("CHIDIST", dao);
char const *label;
char *cc;
label = (info->independence)
/* translator info: The quotation marks in the next strings need to */
/* remain since these are Excel-style format strings */
? _("[>=5]\"Test of Independence\";[<5][Red]\"Invalid Test of Independence\"")
: _("[>=5]\"Test of Homogeneity\";[<5][Red]\"Invalid Test of Homogeneity\"");
dao_set_italic (dao, 0, 1, 0, 4);
set_cell_text_col (dao, 0, 1, _("/Test Statistic"
"/Degrees of Freedom"
"/p-Value"
"/Critical Value"));
cc = g_strdup_printf ("%s = %.2" GNM_FORMAT_f, "\xce\xb1", info->alpha);
dao_set_cell_comment (dao, 0, 4, cc);
g_free (cc);
if (info->labels)
expr_region = gnm_expr_new_funcall5
(fd_offset,
gnm_expr_new_constant (value_dup (info->input)),
gnm_expr_new_constant (value_new_int (1)),
gnm_expr_new_constant (value_new_int (1)),
gnm_expr_new_constant (value_new_int (info->n_r)),
gnm_expr_new_constant (value_new_int (info->n_c)));
else
expr_region = gnm_expr_new_constant (value_dup (info->input));
expr_row = gnm_expr_new_funcall1 (fd_row, gnm_expr_copy (expr_region));
expr_column = gnm_expr_new_funcall1 (fd_column, gnm_expr_copy (expr_region));
expr_col_ones = gnm_expr_new_funcall1 (fd_transpose,
gnm_expr_new_binary (gnm_expr_copy (expr_column),
GNM_EXPR_OP_DIV,
expr_column));
expr_row_ones = gnm_expr_new_funcall1 (fd_transpose,
gnm_expr_new_binary (gnm_expr_copy (expr_row),
GNM_EXPR_OP_DIV,
expr_row));
expr_expect = gnm_expr_new_binary (gnm_expr_new_funcall2
(fd_mmult,
gnm_expr_new_funcall2
(fd_mmult,
gnm_expr_copy (expr_region),
expr_col_ones),
gnm_expr_new_funcall2
(fd_mmult,
expr_row_ones,
gnm_expr_copy (expr_region))),
GNM_EXPR_OP_DIV,
gnm_expr_new_funcall1 (fd_sum, gnm_expr_copy (expr_region)));
expr_check = gnm_expr_new_funcall1 (fd_min, gnm_expr_copy (expr_expect));
dao_set_merge (dao, 0, 0, 1, 0);
dao_set_italic (dao, 0, 0, 0, 0);
dao_set_cell_expr (dao, 0, 0, expr_check);
dao_set_format (dao, 0, 0, 0, 0, label);
dao_set_align (dao, 0, 0, 0, 0, GNM_HALIGN_CENTER, GNM_VALIGN_BOTTOM);
expr_statistic = gnm_expr_new_funcall1 (fd_sum,
gnm_expr_new_binary
(gnm_expr_new_binary (gnm_expr_new_binary
(gnm_expr_copy (expr_region),
GNM_EXPR_OP_SUB,
gnm_expr_copy (expr_expect)),
GNM_EXPR_OP_EXP,
gnm_expr_new_constant (value_new_int (2))),
GNM_EXPR_OP_DIV,
gnm_expr_copy (expr_expect)));
dao_set_cell_array_expr (dao, 1, 1, expr_statistic);
dao_set_cell_int (dao, 1, 2, (info->n_r - 1)*(info->n_c - 1));
dao_set_cell_expr(dao, 1, 3, gnm_expr_new_funcall2
(fd_chidist, make_cellref (0,-2), make_cellref (0,-1)));
dao_set_cell_expr(dao, 1, 4, gnm_expr_new_funcall2
(fd_chiinv,
gnm_expr_new_constant (value_new_float (info->alpha)),
make_cellref (0,-2)));
gnm_func_unref (fd_mmult);
gnm_func_unref (fd_row);
gnm_func_unref (fd_column);
gnm_func_unref (fd_transpose);
gnm_func_unref (fd_sum);
gnm_func_unref (fd_min);
gnm_func_unref (fd_offset);
gnm_func_unref (fd_chiinv);
gnm_func_unref (fd_chidist);
gnm_expr_free (expr_expect);
gnm_expr_free (expr_region);
dao_redraw_respan (dao);
return FALSE;
}
static gboolean
analysis_tool_principal_components_engine_run (data_analysis_output_t *dao,
analysis_tools_data_generic_t *info)
{
int l = g_slist_length (info->input), i;
GSList *inputdata;
GnmFunc *fd_mean;
GnmFunc *fd_var;
GnmFunc *fd_eigen;
GnmFunc *fd_mmult;
GnmFunc *fd_munit;
GnmFunc *fd_sqrt;
GnmFunc *fd_count;
GnmFunc *fd_sum;
GnmFunc *fd_and;
GnmFunc *fd_if;
GnmExpr const *expr;
GnmExpr const *expr_count;
GnmExpr const *expr_munit;
GnmExpr const *expr_and;
int data_points;
GnmExprList *and_args = NULL;
if (!dao_cell_is_visible (dao, l, 9 + 3 * l)) {
dao_set_bold (dao, 0, 0, 0, 0);
dao_set_italic (dao, 0, 0, 0, 0);
dao_set_cell (dao, 0, 0,
_("Principal components analysis has "
"insufficient space."));
return 0;
}
fd_mean = gnm_func_lookup_or_add_placeholder
("AVERAGE", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_mean);
fd_var = gnm_func_lookup_or_add_placeholder
("VAR", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_var);
fd_eigen = gnm_func_lookup_or_add_placeholder
("EIGEN", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_eigen);
fd_mmult = gnm_func_lookup_or_add_placeholder
("MMULT", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_mmult);
fd_munit = gnm_func_lookup_or_add_placeholder
("MUNIT", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_munit);
fd_sqrt = gnm_func_lookup_or_add_placeholder
("SQRT", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_sqrt);
fd_count = gnm_func_lookup_or_add_placeholder
("COUNT", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_count);
fd_sum = gnm_func_lookup_or_add_placeholder
("SUM", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_sum);
fd_and = gnm_func_lookup_or_add_placeholder
("AND", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_and);
fd_if = gnm_func_lookup_or_add_placeholder
("IF", dao->sheet ? dao->sheet->workbook : NULL, FALSE);
gnm_func_ref (fd_if);
dao_set_bold (dao, 0, 0, 0, 0);
dao_set_italic (dao, 0, 0, 0, 11 + 3 * l);
dao_set_format (dao, 0, 0, 0, 0,
_("\"Principal Components Analysis\";"
"[Red]\"Principal Components Analysis is invalid.\""));
dao_set_align (dao, 0, 0, 0, 0,
HALIGN_LEFT, VALIGN_BOTTOM);
dao->offset_row++;
analysis_tool_table (dao, info, _("Covariances:"), "COVAR", TRUE);
dao->offset_row--;
for (i = 1, inputdata = info->input; inputdata != NULL; i++, inputdata = inputdata->next)
analysis_tools_write_label (inputdata->data, dao, info, 0, 9 + 2 * l + i, i);
data_points = value_area_get_width (info->input->data, NULL) *
value_area_get_height (info->input->data, NULL);
for (i = 0; i < l; i++)
and_args = gnm_expr_list_prepend
(and_args,
gnm_expr_new_binary
(gnm_expr_new_constant (value_new_int (data_points)),
GNM_EXPR_OP_EQUAL,
make_cellref (1 + i, 3 + l)));
expr_and = gnm_expr_new_funcall (fd_and, and_args);
dao_set_cell_expr (dao, 0, 0,
gnm_expr_new_funcall3
(fd_if,
expr_and,
gnm_expr_new_constant (value_new_int (1)),
gnm_expr_new_constant (value_new_int (-1))));
dao_set_merge (dao,0,0,2,0);
set_cell_text_col (dao, 0, 3 + l,
_("/Count:"
"/Mean:"
"/Variance:"
"//Eigenvalues:"
"/Eigenvectors:"));
dao_set_cell (dao, 0, 11 + 3 * l, _("Percent of Trace:"));
dao_set_italic (dao, 0, 9 + 2 * l, 1 + l, 9 + 2 * l);
dao_set_percent (dao, 1, 11 + 3 * l, 1 + l, 11 + 3 * l);
for (i = 1, inputdata = info->input; inputdata != NULL; i++, inputdata = inputdata->next) {
expr = gnm_expr_new_constant (value_dup (inputdata->data));
dao_set_cell_expr (dao, i, 3 + l,
gnm_expr_new_funcall1 (fd_count, gnm_expr_copy (expr)));
dao_set_cell_expr (dao, i, 4 + l,
gnm_expr_new_funcall1 (fd_mean, gnm_expr_copy (expr)));
dao_set_cell_expr (dao, i, 5 + l,
gnm_expr_new_funcall1 (fd_var, expr));
}
expr_count = gnm_expr_new_binary (make_cellref (0,-4), GNM_EXPR_OP_DIV,
gnm_expr_new_binary (make_cellref (0,-4), GNM_EXPR_OP_SUB,
gnm_expr_new_constant (value_new_int (1))));
expr = gnm_expr_new_funcall1
(fd_eigen, gnm_expr_new_binary
(expr_count, GNM_EXPR_OP_MULT, make_rangeref (0, - (5 + l), l - 1, - 6)));
dao_set_array_expr (dao, 1, 7 + l, l, l + 1, expr);
for (i = 1; i <= l; i++) {
dao_set_align (dao, i, 9 + 2 * l, i, 9 + 2 * l,
HALIGN_CENTER, VALIGN_BOTTOM);
dao_set_cell_printf (dao, i, 9 + 2 * l, "\xce\xbe%i", i);
dao_set_cell_expr (dao, i, 11 + 3 * l,
gnm_expr_new_binary (make_cellref (0,- 4 - 2 * l),
GNM_EXPR_OP_DIV,
gnm_expr_new_funcall1
(fd_sum,
dao_get_rangeref (dao, 1, 7 + l, l, 7 + l))));
}
expr_munit = gnm_expr_new_funcall1 (fd_munit, gnm_expr_new_constant (value_new_int (l)));
expr = gnm_expr_new_funcall2 (fd_mmult,
gnm_expr_new_binary
(gnm_expr_new_funcall1
(fd_sqrt, gnm_expr_new_binary
(gnm_expr_new_constant (value_new_int (1)),
GNM_EXPR_OP_DIV,
make_rangeref (0, - 5 - l, l - 1, - 5 - l))),
GNM_EXPR_OP_MULT,
gnm_expr_copy (expr_munit)),
make_rangeref (0, - 2 - l, l - 1, - 3));
expr = gnm_expr_new_funcall2 (fd_mmult, expr,
gnm_expr_new_binary
(gnm_expr_new_funcall1
(fd_sqrt, make_rangeref (0, - 3 - l, l - 1, - 3 - l)),
GNM_EXPR_OP_MULT,
expr_munit));
dao_set_array_expr (dao, 1, 10 + 2 * l, l, l, expr);
gnm_func_unref (fd_mean);
gnm_func_unref (fd_var);
gnm_func_unref (fd_eigen);
gnm_func_unref (fd_mmult);
gnm_func_unref (fd_munit);
gnm_func_unref (fd_sqrt);
gnm_func_unref (fd_count);
gnm_func_unref (fd_sum);
gnm_func_unref (fd_and);
gnm_func_unref (fd_if);
dao_redraw_respan (dao);
return 0;
}
static gboolean
analysis_tool_frequency_engine_run (data_analysis_output_t *dao,
analysis_tools_data_frequency_t *info)
{
gint i_limit, col;
GSList *l;
GnmFunc *fd_sum;
GnmFunc *fd_if;
GnmFunc *fd_index;
GnmFunc *fd_isblank;
GnmFunc *fd_rows = NULL;
GnmFunc *fd_columns = NULL;
GnmFunc *fd_exact = NULL;
fd_sum = gnm_func_lookup_or_add_placeholder ("SUM");
gnm_func_ref (fd_sum);
fd_if = gnm_func_lookup_or_add_placeholder ("IF");
gnm_func_ref (fd_if);
fd_index = gnm_func_lookup_or_add_placeholder ("INDEX");
gnm_func_ref (fd_index);
fd_isblank = gnm_func_lookup_or_add_placeholder ("ISBLANK");
gnm_func_ref (fd_isblank);
if (info->exact) {
fd_exact = gnm_func_lookup_or_add_placeholder ("EXACT");
gnm_func_ref (fd_exact);
}
if (info->percentage) {
fd_rows = gnm_func_lookup_or_add_placeholder ("ROWS");
gnm_func_ref (fd_rows);
fd_columns = gnm_func_lookup_or_add_placeholder ("COLUMNS");
gnm_func_ref (fd_columns);
}
/* General Info */
dao_set_italic (dao, 0, 0, 0, 1);
set_cell_text_col (dao, 0, 0, _("/Frequency Table"
"/Category"));
/* Setting up the categories */
if (info->predetermined) {
int row = 2, i, j, i_h_limit, i_w_limit;
GnmExpr const *expr_bin;
GnmRange range;
range_init_value (&range, info->bin);
i_h_limit = range_height (&range);
i_w_limit = range_width (&range);
i_limit = i_h_limit * i_w_limit;
expr_bin = gnm_expr_new_constant (info->bin);
for (i = 1; i <= i_h_limit; i++)
for (j = 1; j <= i_w_limit; j++) {
GnmExpr const *expr_index;
expr_index = gnm_expr_new_funcall3
(fd_index,
gnm_expr_copy (expr_bin),
gnm_expr_new_constant (value_new_int (i)),
gnm_expr_new_constant (value_new_int (j)));
dao_set_cell_expr (dao, 0, row++,
gnm_expr_new_funcall3
(fd_if,
gnm_expr_new_funcall1
(fd_isblank,
gnm_expr_copy (expr_index)),
gnm_expr_new_constant (value_new_string ("")),
expr_index));
}
gnm_expr_free (expr_bin);
} else {
i_limit = info->n;
}
for (l = info->base.input, col = 1; l; col++, l = l->next) {
GnmValue *val = value_dup ((GnmValue *)l->data);
GnmValue *val_c = NULL;
GnmExpr const *expr_count;
GnmExpr const *expr_data;
GnmExpr const *expr_if;
int i, row = 2;
dao_set_italic (dao, col, 1, col, 1);
if (info->base.labels) {
val_c = value_dup (val);
switch (info->base.group_by) {
case GROUPED_BY_ROW:
val->v_range.cell.a.col++;
break;
default:
val->v_range.cell.a.row++;
break;
}
dao_set_cell_expr (dao, col, 1,
gnm_expr_new_funcall1 (fd_index,
gnm_expr_new_constant (val_c)));
} else {
char const *format;
switch (info->base.group_by) {
case GROUPED_BY_ROW:
format = _("Row %d");
break;
case GROUPED_BY_COL:
format = _("Column %d");
break;
default:
format = _("Area %d");
break;
}
dao_set_cell_printf (dao, col, 1, format, col);
}
expr_data = gnm_expr_new_constant (val);
if (info->exact)
expr_if = gnm_expr_new_funcall2
(fd_exact, gnm_expr_copy (expr_data),
make_cellref (- col, 0));
else
expr_if = gnm_expr_new_binary
(gnm_expr_copy (expr_data),
GNM_EXPR_OP_EQUAL, make_cellref (- col, 0));
expr_count = gnm_expr_new_funcall1 (fd_sum,
gnm_expr_new_funcall3
(fd_if, expr_if,
gnm_expr_new_constant (value_new_int (1)),
gnm_expr_new_constant (value_new_int (0))));
if (info->percentage) {
dao_set_format (dao, col, 2, col, i_limit + 2, "0.0%");
expr_count = gnm_expr_new_binary (expr_count,
GNM_EXPR_OP_DIV,
gnm_expr_new_binary
(gnm_expr_new_funcall1
(fd_rows, gnm_expr_copy (expr_data)),
GNM_EXPR_OP_MULT,
gnm_expr_new_funcall1
(fd_columns, expr_data)));
} else
gnm_expr_free (expr_data);
for (i = 0; i < i_limit; i++, row++)
dao_set_cell_array_expr (dao, col, row, gnm_expr_copy (expr_count));
gnm_expr_free (expr_count);
}
gnm_func_unref (fd_if);
gnm_func_unref (fd_sum);
gnm_func_unref (fd_index);
gnm_func_unref (fd_isblank);
if (fd_rows != NULL)
gnm_func_unref (fd_rows);
if (fd_columns != NULL)
gnm_func_unref (fd_columns);
if (fd_exact != NULL)
gnm_func_unref (fd_exact);
/* Create Chart if requested */
if (info->chart != NO_CHART) {
SheetObject *so;
GogGraph *graph;
GogChart *chart;
GogPlot *plot;
GogSeries *series;
GOData *cats;
GOData *values;
int ct;
graph = g_object_new (GOG_TYPE_GRAPH, NULL);
chart = GOG_CHART (gog_object_add_by_name (
GOG_OBJECT (graph), "Chart", NULL));
plot = gog_plot_new_by_name ("GogBarColPlot");
if (info->chart == BAR_CHART)
go_object_toggle (plot, "horizontal");
gog_object_add_by_name (GOG_OBJECT (chart),
"Plot", GOG_OBJECT (plot));
cats = dao_go_data_vector (dao, 0, 2,
0, 2 + i_limit);
for (ct = 1; ct < col; ct ++) {
g_object_ref (cats);
values = dao_go_data_vector (dao, ct, 2,
ct, 2 + i_limit);
series = gog_plot_new_series (plot);
gog_series_set_dim (series, 0, cats, NULL);
gog_series_set_dim (series, 1, values, NULL);
}
g_object_unref (cats);
so = sheet_object_graph_new (graph);
g_object_unref (graph);
dao_set_sheet_object (dao, 0, 1, so);
}
dao_redraw_respan (dao);
return FALSE;
}
