home *** CD-ROM | disk | FTP | other *** search
- stchi2 (a,b,c,d,chi2)
-
- /*this subroutine computes chi square for 2 by 2 tables*/
- /*with yates correction for continuity.*/
-
- float a,b,c,d,*chi2;
-
- {
- float r1,r2,c1,c2,s,e,f,g;
- extern double fabs();
-
- r1 = a + b;
- r2 = c + d;
- c1 = a + c;
- c2 = b + d;
- s = r1 + r2;
- e = fabs(a*d-b*c);
- f = e-0.5*s;
- g = s*f*f;
- *chi2 = g/(r1*r2*c1*c2);
- }