Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/twoav.ssp
There are 2 other files named twoav.ssp in the archive. Click here to see a list.
C                                                                       TWOA  10
C     ..................................................................TWOA  20
C                                                                       TWOA  30
C        SUBROUTINE TWOAV                                               TWOA  40
C                                                                       TWOA  50
C        PURPOSE                                                        TWOA  60
C           TEST WHETHER A NUMBER OF SAMPLES ARE FROM THE SAME          TWOA  70
C           POPULATION BY THE FRIEDMAN TWO-WAY ANALYSIS OF VARIANCE TESTTWOA  80
C                                                                       TWOA  90
C        USAGE                                                          TWOA 100
C           CALL TWOAV(A,R,N,M,W,XR,NDF,NR)                             TWOA 110
C                                                                       TWOA 120
C        DESCRIPTION OF PARAMETERS                                      TWOA 130
C           A   - INPUT MATRIX, N BY M, OF ORIGINAL DATA                TWOA 140
C           R   - OUTPUT MATRIX, N BY M, OF RANKED DATA                 TWOA 150
C           N   - NUMBER OF GROUPS                                      TWOA 160
C           M   - NUMBER OF CASES IN EACH GROUP                         TWOA 170
C           W   - WORK AREA OF LENGTH 2*M                               TWOA 180
C           XR  - FRIEDMAN STATISTIC (OUTPUT)                           TWOA 190
C           NDF - NUMBER OF DEGREES OF FREEDOM (OUTPUT)                 TWOA 200
C           NR  - CODE, 0 FOR UNRANKED DATA IN A, 1 FOR RANKED DATA     TWOA 210
C                 IN A (INPUT)                                          TWOA 220
C                                                                       TWOA 230
C        REMARKS                                                        TWOA 240
C           NONE                                                        TWOA 250
C                                                                       TWOA 260
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  TWOA 270
C           RANK                                                        TWOA 280
C                                                                       TWOA 290
C        METHOD                                                         TWOA 300
C           DESCRIBED IN S. SIEGEL, 'NONPARAMETRIC STATISTICS FOR THE   TWOA 310
C           BEHAVIORAL SCIENCES', MCGRAW-HILL, NEW YORK, 1956,          TWOA 320
C           CHAPTER 7                                                   TWOA 330
C                                                                       TWOA 340
C     ..................................................................TWOA 350
C                                                                       TWOA 360
      SUBROUTINE TWOAV (A,R,N,M,W,XR,NDF,NR)                            TWOA 370
      DIMENSION A(1),R(1),W(1)                                          TWOA 380
C                                                                       TWOA 390
C        DETERMINE WHETHER DATA IS RANKED                               TWOA 400
C                                                                       TWOA 410
      IF(NR-1) 10, 30, 10                                               TWOA 420
C                                                                       TWOA 430
C        RANK DATA IN EACH GROUP AND ASSIGN TIED OBSERVATIONS AVERAGE   TWOA 440
C        OF TIED RANK                                                   TWOA 450
C                                                                       TWOA 460
   10 DO 20 I=1,N                                                       TWOA 470
      IJ=I-N                                                            TWOA 480
      IK=IJ                                                             TWOA 490
      DO 15 J=1,M                                                       TWOA 500
      IJ=IJ+N                                                           TWOA 510
   15 W(J)=A(IJ)                                                        TWOA 520
      CALL RANK (W,W(M+1),M)                                            TWOA 530
      DO 20 J=1,M                                                       TWOA 540
      IK=IK+N                                                           TWOA 550
      IW=M+J                                                            TWOA 560
   20 R(IK)=W(IW)                                                       TWOA 570
      GO TO 35                                                          TWOA 580
   30 NM=N*M                                                            TWOA 590
      DO 32 I=1,NM                                                      TWOA 600
   32 R(I)=A(I)                                                         TWOA 610
C                                                                       TWOA 620
C        CALCULATE SUM OF SQUARES OF SUMS OF RANKS                      TWOA 630
C                                                                       TWOA 640
   35 RTSQ=0.0                                                          TWOA 650
      IR=0                                                              TWOA 660
      DO 50 J=1,M                                                       TWOA 670
      RT=0.0                                                            TWOA 680
      DO 40 I=1,N                                                       TWOA 690
      IR=IR+1                                                           TWOA 700
   40 RT=RT+R(IR)                                                       TWOA 710
   50 RTSQ=RTSQ+RT*RT                                                   TWOA 720
C                                                                       TWOA 730
C        CALCULATE FRIEDMAN TEST VALUE, XR                              TWOA 740
C                                                                       TWOA 750
      FNM=N*(M+1)                                                       TWOA 760
      FM=M                                                              TWOA 770
      XR=(12.0/(FM*FNM))*RTSQ-3.0*FNM                                   TWOA 780
C                                                                       TWOA 790
C        FIND DEGREES OF FREEDOM                                        TWOA 800
C                                                                       TWOA 810
      NDF=M-1                                                           TWOA 820
      RETURN                                                            TWOA 830
      END                                                               TWOA 840