1 /* PSPP - computes sample statistics.
2 Copyright (C) 1997-9, 2000 Free Software Foundation, Inc.
3 Written by Ben Pfaff <blp@gnu.org>.
5 This program is free software; you can redistribute it and/or
6 modify it under the terms of the GNU General Public License as
7 published by the Free Software Foundation; either version 2 of the
8 License, or (at your option) any later version.
10 This program is distributed in the hope that it will be useful, but
11 WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
23 #include <language/command.h>
24 #include <language/lexer/lexer.h>
25 #include <math/moments.h>
29 #include <libpspp/compiler.h>
31 #define _(msgid) gettext (msgid)
34 read_values (struct lexer *lexer, double **values, double **weights, size_t *cnt)
41 while (lex_is_number (lexer))
43 double value = lex_tokval (lexer);
46 if (lex_match (lexer, '*'))
48 if (!lex_is_number (lexer))
50 lex_error (lexer, _("expecting weight value"));
53 weight = lex_tokval (lexer);
60 *values = xnrealloc (*values, cap, sizeof **values);
61 *weights = xnrealloc (*weights, cap, sizeof **weights);
64 (*values)[*cnt] = value;
65 (*weights)[*cnt] = weight;
73 cmd_debug_moments (struct lexer *lexer, struct dataset *ds UNUSED)
75 int retval = CMD_FAILURE;
76 double *values = NULL;
77 double *weights = NULL;
83 if (lex_match_id (lexer, "ONEPASS"))
85 if (lex_token (lexer) != '/')
87 lex_force_match (lexer, '/');
90 fprintf (stderr, "%s => ", lex_rest_of_line (lexer, NULL));
95 struct moments *m = NULL;
97 m = moments_create (MOMENT_KURTOSIS);
98 if (!read_values (lexer, &values, &weights, &cnt))
103 for (i = 0; i < cnt; i++)
104 moments_pass_one (m, values[i], weights[i]);
105 for (i = 0; i < cnt; i++)
106 moments_pass_two (m, values[i], weights[i]);
107 moments_calculate (m, &weight, &M[0], &M[1], &M[2], &M[3]);
112 struct moments1 *m = NULL;
114 m = moments1_create (MOMENT_KURTOSIS);
115 if (!read_values (lexer, &values, &weights, &cnt))
117 moments1_destroy (m);
120 for (i = 0; i < cnt; i++)
121 moments1_add (m, values[i], weights[i]);
122 moments1_calculate (m, &weight, &M[0], &M[1], &M[2], &M[3]);
123 moments1_destroy (m);
126 fprintf (stderr, "W=%.3f", weight);
127 for (i = 0; i < 4; i++)
129 fprintf (stderr, " M%d=", i + 1);
131 fprintf (stderr, "sysmis");
132 else if (fabs (M[i]) <= 0.0005)
133 fprintf (stderr, "0.000");
135 fprintf (stderr, "%.3f", M[i]);
137 fprintf (stderr, "\n");
139 retval = lex_end_of_command (lexer);