home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / STANV1.C < prev    next >
Encoding:
Text File  |  1984-08-17  |  1.0 KB  |  47 lines
  1.    stanv1 (n,m,k,x,f,av)
  2.  
  3.       /* this function uses the data of a one-way, completely randomized */
  4.       /* experiment to compute the analysis of variance table and the    */
  5.       /* f - ratio to test the significance of the difference between    */
  6.       /* treatment means.                                                */
  7.  
  8.       int k[],m,n;
  9.       float av[],x[],*f;
  10.  
  11.    {
  12.       int i,j,j1,kk;
  13.       float t,u,v,xk;
  14.  
  15.       av[0] = 0.;
  16.       u = 0.;
  17.       v = 0.;
  18.       kk = 0;
  19.  
  20.       for(i = 0; i <= m-1; i++)
  21.       {
  22.        t = 0.;
  23.        j1 = kk + 1;
  24.        kk = kk + k[i];
  25.        xk = k[i];
  26.  
  27.        for(j = j1-1; j <= kk-1; j++)
  28.        {
  29.         av[0] = av[0] + x[j] * x[j];
  30.         t = t + x[j];
  31.        }
  32.        u = u + t;
  33.        v = v + t * t / xk;
  34.       }
  35.  
  36.       av[1] = n - 1;
  37.       av[2] = u * u/av[1];
  38.       av[3] = 1.;
  39.       av[4] = v - av[2];
  40.       av[5] = m - 1;
  41.       av[6] = av[4]/av[5];
  42.       av[7] = av[0] - v;
  43.       av[8] = n - m;
  44.       av[9] = av[7]/av[8];
  45.       *f = av[6]/av[9];
  46.    }
  47.