Lines Matching refs:ds
153 struct dataset *ds;
155 ds = calloc(1, sizeof *ds);
156 assert(ds != NULL);
157 ds->lpoints = 100000;
158 ds->points = calloc(sizeof *ds->points, ds->lpoints);
159 assert(ds->points != NULL);
160 ds->syy = NAN;
161 return(ds);
165 AddPoint(struct dataset *ds, double a)
169 if (ds->n >= ds->lpoints) {
170 dp = ds->points;
171 ds->lpoints *= 4;
172 ds->points = calloc(sizeof *ds->points, ds->lpoints);
173 assert(ds->points != NULL);
174 memcpy(ds->points, dp, sizeof *dp * ds->n);
177 ds->points[ds->n++] = a;
178 ds->sy += a;
182 Min(const struct dataset *ds)
185 return (ds->points[0]);
189 Max(const struct dataset *ds)
192 return (ds->points[ds->n -1]);
196 Avg(const struct dataset *ds)
199 return(ds->sy / ds->n);
203 Median(const struct dataset *ds)
205 const size_t m = ds->n / 2;
207 if ((ds->n % 2) == 0)
208 return ((ds->points[m] + (ds->points[m - 1])) / 2);
209 return (ds->points[m]);
213 Var(struct dataset *ds)
216 const double a = Avg(ds);
218 if (isnan(ds->syy)) {
219 ds->syy = 0.0;
220 for (z = 0; z < ds->n; z++)
221 ds->syy += (ds->points[z] - a) * (ds->points[z] - a);
224 return (ds->syy / (ds->n - 1.0));
228 Stddev(struct dataset *ds)
231 return sqrt(Var(ds));
242 Vitals(struct dataset *ds, int flag)
246 ds->n, Min(ds), Max(ds), Median(ds), Avg(ds), Stddev(ds));
251 Relative(struct dataset *ds, struct dataset *rs, int confidx)
257 z = ds->n + rs->n - 2;
262 spool = (ds->n - 1) * Var(ds) + (rs->n - 1) * Var(rs);
263 spool /= ds->n + rs->n - 2;
265 s = spool * sqrt(1.0 / ds->n + 1.0 / rs->n);
266 d = Avg(ds) - Avg(rs);
269 re = (ds->n - 1) * Var(ds) + (rs->n - 1) * Var(rs) *
270 (Avg(ds) * Avg(ds)) / (Avg(rs) * Avg(rs));
271 re *= (ds->n + rs->n) / (ds->n * rs->n * (ds->n + rs->n - 2.0));
333 DimPlot(struct dataset *ds)
335 AdjPlot(Min(ds));
336 AdjPlot(Max(ds));
337 AdjPlot(Avg(ds) - Stddev(ds));
338 AdjPlot(Avg(ds) + Stddev(ds));
342 PlotSet(struct dataset *ds, int val)
375 for (n = 0; n < ds->n; n++) {
376 x = (ds->points[n] - pl->x0) / pl->dx;
395 for (n = 0; n < ds->n; n++) {
396 x = (ds->points[n] - pl->x0) / pl->dx;
405 av = Avg(ds);
406 sd = Stddev(ds);
416 x = (Median(ds) - pl->x0) / pl->dx;
557 struct dataset *ds[MAX_DS - 1];
662 ds[i] = ReadSet(setfiles[i], setfilenames[i], column, delim);
668 printf("%c %s\n", symbol[i+1], ds[i]->name);
673 DimPlot(ds[i]);
675 PlotSet(ds[i], i + 1);
679 Vitals(ds[0], 1);
681 Vitals(ds[i], i + 1);
683 Relative(ds[i], ds[0], ci);