X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;ds=sidebyside;f=src%2Flanguage%2Fstats%2Ffactor.c;h=a53e7334240c63aacfa7adad92f7a3ebdf6d2031;hb=0c120e1c2f74ace50987771f1102fa609fd271cd;hp=2dbf3b69143f88ad53a1515670000f42f09e0dc3;hpb=204a1ee35aebcc2cf955017070c1a3638cdaee22;p=pspp diff --git a/src/language/stats/factor.c b/src/language/stats/factor.c index 2dbf3b6914..a53e733424 100644 --- a/src/language/stats/factor.c +++ b/src/language/stats/factor.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 2009 Free Software Foundation, Inc. + Copyright (C) 2009, 2010 Free Software Foundation, Inc. This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by @@ -44,10 +44,10 @@ #include #include -#include +#include #include -#include +#include #include "gettext.h" #define _(msgid) gettext (msgid) @@ -524,6 +524,8 @@ cmd_factor (struct lexer *lexer, struct dataset *ds) const struct dictionary *dict = dataset_dict (ds); struct cmd_factor factor; + factor.n_vars = 0; + factor.vars = NULL; factor.method = METHOD_CORR; factor.missing_type = MISS_LISTWISE; factor.exclude = MV_ANY; @@ -552,6 +554,9 @@ cmd_factor (struct lexer *lexer, struct dataset *ds) PV_NO_DUPLICATE | PV_NUMERIC)) goto error; + if (factor.n_vars < 2) + msg (MW, _("Factor analysis on a single variable is not useful.")); + while (lex_token (lexer) != '.') { lex_match (lexer, '/'); @@ -935,7 +940,7 @@ show_scree (const struct cmd_factor *f, struct idata *idata) s = scree_create (idata->eval, label); - chart_submit (scree_get_chart (s)); + scree_submit (s); } static void @@ -964,8 +969,6 @@ show_communalities (const struct cmd_factor * factor, tab_title (t, _("Communalities")); - tab_dim (t, tab_natural_dimensions, NULL, NULL); - tab_headers (t, heading_columns, 0, heading_rows, 0); c = 1; @@ -1012,7 +1015,7 @@ static void show_factor_matrix (const struct cmd_factor *factor, struct idata *idata, const gsl_matrix *fm) { int i; - const int n_factors = n_extracted_factors (factor, idata); + const int n_factors = idata->n_extractions; const int heading_columns = 1; const int heading_rows = 2; @@ -1027,8 +1030,6 @@ show_factor_matrix (const struct cmd_factor *factor, struct idata *idata, const else tab_title (t, _("Factor Matrix")); - tab_dim (t, tab_natural_dimensions, NULL, NULL); - tab_headers (t, heading_columns, 0, heading_rows, 0); if ( factor->extraction == EXTRACTION_PC ) @@ -1136,8 +1137,6 @@ show_explained_variance (const struct cmd_factor * factor, struct idata *idata, tab_title (t, _("Total Variance Explained")); - tab_dim (t, tab_natural_dimensions, NULL, NULL); - tab_headers (t, heading_columns, 0, heading_rows, 0); /* Outline the box */ @@ -1187,6 +1186,7 @@ show_explained_variance (const struct cmd_factor * factor, struct idata *idata, for (i = 0; i < (nc - heading_columns) / 3 ; ++i) { tab_text (t, i * 3 + 1, 1, TAB_CENTER | TAT_TITLE, _("Total")); + /* xgettext:no-c-format */ tab_text (t, i * 3 + 2, 1, TAB_CENTER | TAT_TITLE, _("% of Variance")); tab_text (t, i * 3 + 3, 1, TAB_CENTER | TAT_TITLE, _("Cumulative %")); @@ -1231,7 +1231,7 @@ show_explained_variance (const struct cmd_factor * factor, struct idata *idata, if (factor->print & PRINT_EXTRACTION) { - if ( i < n_extracted_factors (factor, idata)) + if (i < idata->n_extractions) { /* Sums of squared loadings */ tab_double (t, c++, i + heading_rows, 0, e_lambda, NULL); @@ -1288,8 +1288,6 @@ show_correlation_matrix (const struct cmd_factor *factor, const struct idata *id tab_title (t, _("Correlation Matrix")); - tab_dim (t, tab_natural_dimensions, NULL, NULL); - tab_hline (t, TAL_1, 0, nc - 1, heading_rows); if (nr > heading_rows) @@ -1431,7 +1429,6 @@ do_factor (const struct cmd_factor *factor, struct casereader *r) struct tab_table *t = tab_create (nc, nr); tab_title (t, _("Descriptive Statistics")); - tab_dim (t, tab_natural_dimensions, NULL, NULL); tab_headers (t, heading_columns, 0, heading_rows, 0); @@ -1483,12 +1480,26 @@ do_factor (const struct cmd_factor *factor, struct casereader *r) gsl_eigen_symmv_sort (idata->eval, idata->evec, GSL_EIGEN_SORT_ABS_DESC); #endif + idata->n_extractions = n_extracted_factors (factor, idata); + + if (idata->n_extractions == 0) + { + msg (MW, _("The FACTOR criteria result in zero factors extracted. Therefore no analysis will be performed.")); + goto finish; + } + + if (idata->n_extractions > factor->n_vars) + { + msg (MW, _("The FACTOR criteria result in more factors than variables, which is not meaningful. No analysis will be performed.")); + goto finish; + } + { const gsl_vector *extracted_eigenvalues = NULL; gsl_vector *initial_communalities = gsl_vector_alloc (factor->n_vars); gsl_vector *extracted_communalities = gsl_vector_alloc (factor->n_vars); size_t i; - struct factor_matrix_workspace *fmw = factor_matrix_workspace_alloc (idata->msr->size, n_extracted_factors (factor, idata)); + struct factor_matrix_workspace *fmw = factor_matrix_workspace_alloc (idata->msr->size, idata->n_extractions); gsl_matrix *factor_matrix = gsl_matrix_calloc (factor->n_vars, fmw->n_factors); if ( factor->extraction == EXTRACTION_PAF) @@ -1527,10 +1538,9 @@ do_factor (const struct cmd_factor *factor, struct casereader *r) } else if (factor->extraction == EXTRACTION_PC) { - for (i = 0 ; i < factor->n_vars; ++i) - { - gsl_vector_set (initial_communalities, i, communality (idata, i, factor->n_vars)); - } + for (i = 0; i < factor->n_vars; ++i) + gsl_vector_set (initial_communalities, i, communality (idata, i, factor->n_vars)); + gsl_vector_memcpy (extracted_communalities, initial_communalities); iterate_factor_matrix (analysis_matrix, extracted_communalities, factor_matrix, fmw); @@ -1551,6 +1561,8 @@ do_factor (const struct cmd_factor *factor, struct casereader *r) gsl_vector_free (extracted_communalities); } + finish: + idata_free (idata); casereader_destroy (r);