if ( cmd->a_plot[XMN_PLT_NPPLOT] )
{
- metric->np = np_create (metric->moments);
+ double n, mean, var;
+ moments1_calculate (metric->moments,
+ &n, &mean, &var, NULL, NULL);
+
+ metric->np = np_create (n, mean, var);
n_os ++;
}
}
struct np *
-np_create (const struct moments1 *m)
+np_create (double n, double mean, double var)
{
- double variance;
struct np *np = xzalloc (sizeof (*np));
struct order_stats *os = &np->parent;
struct statistic *stat = &os->parent;
np->prev_cc = 0;
- moments1_calculate (m, &np->n, &np->mean, &variance, NULL, NULL);
+ np->n = n;
+ np->mean = mean;
- np->stddev = sqrt (variance);
+ np->stddev = sqrt (var);
np->y_min = np->ns_min = np->dns_min = DBL_MAX;
np->y_max = np->ns_max = np->dns_max = -DBL_MAX;
#include "order-stats.h"
-struct moments1;
struct casewriter;
enum
double mean;
double stddev;
-
double prev_cc;
double ns_min;
};
-struct np * np_create (const struct moments1 *);
+struct np * np_create (double n, double mean, double var);
#endif