home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / STCORR.C < prev    next >
Encoding:
Text File  |  1984-07-27  |  661 b   |  31 lines
  1.       stcorr (n,x,y,r)
  2.  
  3.       /*this subroutine computes a correlation coefficient for n*/
  4.       /*pairs of values.*/
  5.  
  6.       int n;
  7.       float x[],y[],*r;
  8.  
  9.      {
  10.  
  11.       extern double sqrt();
  12.       float sumx,sumy,sumxy,sumx2,sumy2,en,t;
  13.       int i;
  14.  
  15.       sumx = sumy = sumxy = sumx2 = sumy2 = 0.0;
  16.  
  17.       for (i = 0; i <= n-1; i++)
  18.       {
  19.         sumx = sumx + x[i];
  20.         sumy = sumy + y[i];
  21.         sumxy = sumxy + x[i]*y[i];
  22.         sumx2 = sumx2 + x[i]*x[i];
  23.         sumy2 = sumy2 + y[i]*y[i];
  24.       }
  25.  
  26.       en = n;
  27.       t = (en*sumx2-sumx*sumx)*(en*sumy2-sumy*sumy);
  28.       *r = (en*sumxy-sumx*sumy)/sqrt(t);
  29.  
  30.      }
  31.