X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fmoments.c;h=517219d54fcf5d3cc159b4b64c8f6454f64a5bd4;hb=919daf4b00ff2072c0fc957ae8554d047336f984;hp=40180f7b7ecdb461c46eaa31b42c1927835c7b66;hpb=b3fcf4b1644bf4af9b5eb7b0b0f8856c51118128;p=pspp diff --git a/src/math/moments.c b/src/math/moments.c index 40180f7b7e..517219d54f 100644 --- a/src/math/moments.c +++ b/src/math/moments.c @@ -28,8 +28,6 @@ #include "gl/xalloc.h" -#include "gettext.h" -#define _(msgid) gettext (msgid) /* Calculates variance, skewness, and kurtosis into *VARIANCE, *SKEWNESS, and *KURTOSIS if they are non-null and not greater @@ -246,14 +244,14 @@ moments_destroy (struct moments *m) free (m); } -/* Calculates the requested moments on the CNT values in ARRAY. +/* Calculates the requested moments on the N values in ARRAY. Each value is given a weight of 1. The total weight is stored into *WEIGHT (trivially) and the mean, variance, skewness, and kurtosis are stored into *MEAN, *VARIANCE, *SKEWNESS, and *KURTOSIS, respectively. Any of the result pointers may be null, in which case no value is stored. */ void -moments_of_doubles (const double *array, size_t cnt, +moments_of_doubles (const double *array, size_t n, double *weight, double *mean, double *variance, double *skewness, double *kurtosis) @@ -272,21 +270,21 @@ moments_of_doubles (const double *array, size_t cnt, max_moment = MOMENT_MEAN; init_moments (&m, max_moment); - for (idx = 0; idx < cnt; idx++) + for (idx = 0; idx < n; idx++) moments_pass_one (&m, array[idx], 1.); - for (idx = 0; idx < cnt; idx++) + for (idx = 0; idx < n; idx++) moments_pass_two (&m, array[idx], 1.); moments_calculate (&m, weight, mean, variance, skewness, kurtosis); } -/* Calculates the requested moments on the CNT numeric values in +/* Calculates the requested moments on the N numeric values in ARRAY. Each value is given a weight of 1. The total weight is stored into *WEIGHT (trivially) and the mean, variance, skewness, and kurtosis are stored into *MEAN, *VARIANCE, *SKEWNESS, and *KURTOSIS, respectively. Any of the result pointers may be null, in which case no value is stored. */ void -moments_of_values (const union value *array, size_t cnt, +moments_of_values (const union value *array, size_t n, double *weight, double *mean, double *variance, double *skewness, double *kurtosis) @@ -305,9 +303,9 @@ moments_of_values (const union value *array, size_t cnt, max_moment = MOMENT_MEAN; init_moments (&m, max_moment); - for (idx = 0; idx < cnt; idx++) + for (idx = 0; idx < n; idx++) moments_pass_one (&m, array[idx].f, 1.); - for (idx = 0; idx < cnt; idx++) + for (idx = 0; idx < n; idx++) moments_pass_two (&m, array[idx].f, 1.); moments_calculate (&m, weight, mean, variance, skewness, kurtosis); } @@ -459,9 +457,9 @@ moments1_destroy (struct moments1 *m) double -calc_semean (double var, double W) +calc_semean (double variance, double W) { - return sqrt (var / W); + return sqrt (variance / W); } /* Returns the standard error of the skewness for the given total