static gboolean
tool_random_engine_run_bernoulli (data_analysis_output_t *dao,
tools_data_random_t *info,
bernoulli_random_tool_t *param)
{
int i, n;
for (i = 0; i < info->n_vars; i++) {
for (n = 0; n < info->count; n++) {
gnm_float tmp = random_bernoulli (param->p);
dao_set_cell_int (dao, i, n, (int)tmp);
}
}
return FALSE;
}
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;
}
