| 1 |
/* Copyright (c) 1991 Regents of the University of California */ |
| 2 |
|
| 3 |
#ifndef lint |
| 4 |
static char SCCSid[] = "$SunId$ LBL"; |
| 5 |
#endif |
| 6 |
|
| 7 |
/* |
| 8 |
* Basic linear regression calculation. |
| 9 |
*/ |
| 10 |
|
| 11 |
#include "linregr.h" |
| 12 |
|
| 13 |
|
| 14 |
lrclear(l) /* initialize sum */ |
| 15 |
register LRSUM *l; |
| 16 |
{ |
| 17 |
l->xs = l->ys = l->xxs = l->yys = l->xys = 0.0; |
| 18 |
l->n = 0; |
| 19 |
} |
| 20 |
|
| 21 |
|
| 22 |
flrpoint(x, y, l) /* add point (x,y) to sum */ |
| 23 |
double x, y; |
| 24 |
register LRSUM *l; |
| 25 |
{ |
| 26 |
l->xs += x; |
| 27 |
l->ys += y; |
| 28 |
l->xxs += x*x; |
| 29 |
l->yys += y*y; |
| 30 |
l->xys += x*y; |
| 31 |
return(++l->n); |
| 32 |
} |
| 33 |
|
| 34 |
|
| 35 |
lrfit(r, l) /* compute linear regression */ |
| 36 |
register LRLIN *r; |
| 37 |
register LRSUM *l; |
| 38 |
{ |
| 39 |
double nxvar, nyvar; |
| 40 |
|
| 41 |
if (l->n < 2) |
| 42 |
return(-1); |
| 43 |
nxvar = l->xxs - l->xs*l->xs/l->n; |
| 44 |
nyvar = l->yys - l->ys*l->ys/l->n; |
| 45 |
if (nxvar == 0.0 || nyvar == 0.0) |
| 46 |
return(-1); |
| 47 |
r->slope = (l->xys - l->xs*l->ys/l->n) / nxvar; |
| 48 |
r->intercept = (l->ys - r->slope*l->xs) / l->n; |
| 49 |
r->correlation = r->slope*sqrt(nxvar/nyvar); |
| 50 |
return(0); |
| 51 |
} |
