home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / statmisc / phi.for < prev    next >
Encoding:
Text File  |  1985-12-31  |  2.9 KB  |  123 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C     SUBROUTINE PHI
  5. C
  6. C     PURPOSE
  7. C        TO COMPUTE THE PHI COEFFICIENT BETWEEN TWO VARIABLES  WHICH
  8. C        ARE DICHOTOMOUS.
  9. C
  10. C     USAGE
  11. C        CALL PHI (N,U,V,HU,HV,P,CH,XP,IE)
  12. C
  13. C     DESCRIPTION OF PARAMETERS
  14. C        N  - NUMBER OF OBSERVATIONS
  15. C        U  - INPUT VECTOR OF LENGTH N CONTAINING THE FIRST DICHOTO-
  16. C         MOUS VARIABLE
  17. C        V  - INPUT VECTOR OF LENGTH N CONTAINING THE SECOND DICHOTO-
  18. C         MOUS VARIABLE
  19. C        HU - INPUT NUMERICAL CODE WHICH INDICATES THE HIGHER
  20. C         CATEGORY OF THE FIRST VARIABLE.  ANY OBSERVATION IN
  21. C         VECTOR U WHICH HAS A VALUE EQUAL TO OR GREATER THAN HU
  22. C         WILL BE CLASSIFIED IN THE HIGHER CATEGORY.
  23. C        HV - INPUT NUMERICAL CODE FOR VECTOR V, SIMILAR TO HU
  24. C        P  - PHI COEFFICIENT COMPUTED
  25. C        CH - CHI-SQUARE COMPUTED AS A FUNCTION OF PHI COEFFICIENT
  26. C         (DEGREES OF FREEDOM FOR CHI-SQUARE = 1)
  27. C        XP - COMPUTED VALUE OF THE MAXIMAL PHI COEFFICIENT THAT
  28. C         CAN BE ATTAINED IN THE PROBLEM
  29. C        IE - IF IE IS NON-ZERO, SOME CELL IN THE 2 BY 2 TABLE IS
  30. C         NULL.    IF SO, P, CH, AND XP ARE SET TO 10**75.
  31. C
  32. C     REMARKS
  33. C        VARIABLES U AND V MUST BE SPECIFIED NUMERIC.
  34. C        THE PHI COEFFICIENT IS A SPECIAL CASE OF THE
  35. C        PEARSON PRODUCT-MOMENT CORRELATION WHEN BOTH VARIABLES ARE
  36. C        BINARY.
  37. C
  38. C     SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  39. C        NONE
  40. C
  41. C     METHOD
  42. C        REFER TO P. HORST, 'PYSCHOLOGICAL MEASUREMENT AND
  43. C        PREDICTION', P. 94 (WADSWORTH, 1966).
  44. C
  45. C     ..................................................................
  46. C
  47.       SUBROUTINE PHI (N,U,V,HU,HV,P,CH,XP,IE)
  48. C
  49.       DIMENSION U(1),V(1)
  50. C
  51. C     CONSTRUCT A 2X2 CONTINGENCY TABLE
  52. C
  53.       IE=0
  54.       A=0.0
  55.       B=0.0
  56.       C=0.0
  57.       D=0.0
  58. C
  59.       DO 40 I=1,N
  60.       IF(U(I)-HU) 10,25,25
  61.    10 IF(V(I)-HV) 15,20,20
  62.    15 D=D+1.0
  63.       GO TO 40
  64.    20 B=B+1.0
  65.       GO TO 40
  66.    25 IF(V(I)-HV) 30,35,35
  67.    30 C=C+1.0
  68.       GO TO 40
  69.    35 A=A+1.0
  70.    40 CONTINUE
  71.       IF(A) 100,100,41
  72.    41 IF(B) 100,100,42
  73.    42 IF(C) 100,100,43
  74.    43 IF(D) 100,100,44
  75. C
  76. C     COMPUTE THE PHI COEFFICIENT
  77. C
  78.    44 P=(A*D-B*C)/ SQRT((A+B)*(C+D)*(A+C)*(B+D))
  79. C
  80. C     COMPUTE CHI-SQURE
  81. C
  82.       T=N
  83.       CH=T*P*P
  84. C
  85. C     COMPUTE THE MAXIMAL PHI COEFFICIENT
  86. C
  87.       P1=(A+C)/T
  88.       P2=(B+D)/T
  89.       P3=(A+B)/T
  90.       P4=(C+D)/T
  91.       IF(P1-P2) 75, 45, 45
  92.    45 IF(P3-P4) 65, 50, 50
  93.    50 IF(P1-P3) 60, 55, 55
  94.    55 XP=SQRT((P3/P4)*(P2/P1))
  95.       GO TO 95
  96.    60 XP=SQRT((P1/P2)*(P4/P3))
  97.       GO TO 95
  98.    65 IF(P1-P4) 70, 55, 55
  99.    70 XP=SQRT((P2/P1)*(P3/P4))
  100.       GO TO 95
  101.    75 IF(P3-P4) 90, 80, 80
  102.    80 IF(P2-P3) 60, 85, 85
  103.    85 XP=SQRT((P4/P3)*(P1/P2))
  104.       GO TO 95
  105.    90 IF(P2-P4) 70, 85, 85
  106. C
  107.    95 RETURN
  108.   100 IE=1
  109.       P=1.E38
  110.       CH=1.E38
  111.       XP=1.E38
  112.       GO TO 95
  113.       END
  114. 
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121. T((P2/P1)*(P3/P4))
  122.       GO TO 95
  123.    75 IF(P3-P4) 90, 80,