Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-01 - 43,50110/stat08.sta
There are 2 other files named stat08.sta in the archive. Click here to see a list.
100'  NAME--STAT08
110'
120'  DESCRIPTION--COMPARES TWO GROUPS OF DATA USING THE MEDIAN TEST.
130'
140'  SOURCE--UNKNOWN
150'
160'  INSTRUCTIONS--PUT DATA IN LINE 900 AND FOLLOWING.
170'  THE FIRST DATUM IS THE NUMBER OF ENTRIES IN THE FIRST
180'  GROUP, THEN ENTER THE NUMBER OF ENTRIES IN THE SECOND GROUP,
190'  THEN THE FIRST GROUP ITSELF IS ENTERED,AND THEN THE SECOND GROUP.
200'  THE PROGRAM PRINTS OUT THE CHI SQUARE STATISTIC OF A 2 BY 2 TABLE
210'  ON 1 DEGREE OF FREEDOM.
220'  SAMPLE DATA IS IN LINE 900.
230'
240'
250'  *  *  *  *  *  *     MAIN PROGRAM   *  *  *  *  *  *  *  *  *  *
260'
270    READ M, N
280  LET M1 = M + N
290  LET M2 = INT( M1/2 )
300    MAT READ A(M1)
310  LET L = 1
320  LET U = M
330  GOSUB 630
340  LET L = M + 1
350  LET U = M1
360  GOSUB 630
370  LET X = Y = 0
380  LET I = 1
390  LET J = M + 1
400  FOR K = 1 TO M2
410     IF A(I) < A(J) THEN 470
420     LET Y = Y + 1
430     LET J = J + 1
440     IF J <= M1 THEN 520
450     LET X = M2 - Y
460     GO TO 530
470     LET X = X + 1
480     LET I = I + 1
490     IF I <= M THEN 520
500     LET Y = M2 - X
510     GO TO 530
520  NEXT K
530  LET U = M - X
540  LET V = N - Y
550  LET Z = X*V - Y*U
560  LET T = M1*( ABS(Z) - M1/2 ) ^ 2
570  LET C2 = T / M / N / M2 / (M1 - M2)
580    PRINT "TWO SAMPLE MEDIAN TEST."
590    PRINT "GROUP 1  "; X; U
600    PRINT "GROUP 2  "; Y; V
610    PRINT "CHI-SQUARE = "; C2
620  STOP
630  FOR I = 1 TO U - L
640     LET X = A(L)
650     LET Q = L
660     FOR J = L + 1 TO U - I + 1
670        IF X >= A(J) THEN 700
680        LET X = A(J)
690        LET Q = J
700     NEXT J
710     IF J = Q THEN 740
720     LET A(Q) = A(J)
730     LET A(J) = X
740  NEXT I
750  RETURN
900  DATA  4, 6, 160, 160, 140, 190, 117, 145, 147, 120, 150, 120
910END