struct np *np = UP_CAST (s, struct np, parent.parent);
double rank = np->prev_cc + (c + 1) / 2.0;
- double ns = gsl_cdf_ugaussian_Pinv (rank / ( np->n + 1 ));
+ double ns = gsl_cdf_ugaussian_Pinv (rank / (np->n + 1));
double z = (y - np->mean) / np->stddev;
minimize (&np->y_min, y);
cp = case_create (casewriter_get_proto (np->writer));
- case_data_rw_idx (cp, NP_IDX_Y)->f = y;
- case_data_rw_idx (cp, NP_IDX_NS)->f = ns;
- case_data_rw_idx (cp, NP_IDX_DNS)->f = dns;
+ *case_num_rw_idx (cp, NP_IDX_Y) = y;
+ *case_num_rw_idx (cp, NP_IDX_NS) = ns;
+ *case_num_rw_idx (cp, NP_IDX_DNS) = dns;
casewriter_write (np->writer, cp);
np->prev_cc = cc;
struct np *
np_create (double n, double mean, double var)
{
- struct np *np = xzalloc (sizeof (*np));
+ struct np *np = XZALLOC (struct np);
struct order_stats *os = &np->parent;
struct statistic *stat = &os->parent;
struct caseproto *proto;