home *** CD-ROM | disk | FTP | other *** search
/ DP Tool Club 8 / CDASC08.ISO / NEWS / RADIANCE / SRC / COMMON / LINREGR.C < prev    next >
C/C++ Source or Header  |  1993-10-07  |  979b  |  52 lines
  1. /* Copyright (c) 1991 Regents of the University of California */
  2.  
  3. #ifndef lint
  4. static char SCCSid[] = "@(#)linregr.c 2.1 11/12/91 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. }
  52.