1  REM  ****  HP BASIC PROGRAM LIBRARY  *******************************
2  REM
3  REM        STAT08:   COMPARES TWO GROUPS OF DATA USING THE 
4  REM                  MEDIAN TEST
5  REM        36732 REV  A   10/73
6  REM
7  REM  ****  CONTRIBUTED PROGRAM  ************************************
9000  REM  STAT08--MEDIAN TEST FOR 2 GROUPS
9002  DIM A[1000]
9004  READ M,N
9006  LET M1=M+N
9008  LET M2=INT(M1/2)
9010  FOR I=1 TO M1
9012  READ A[I]
9014  NEXT I
9016  LET L=1
9018  LET U=M
9020  GOSUB 9086
9022  LET L=M+1
9024  LET U=M1
9026  GOSUB 9086
9028  LET X=0
9030  LET Y=0
9032  LET I=1
9034  LET J=M+1
9036  FOR K=1 TO M2
9038  IF A[I]<A[J] THEN 9050
9040  LET Y=Y+1
9042  LET J=J+1
9044  IF J <= M1 THEN 9060
9046  LET X=M2-Y
9048  GOTO 9062
9050  LET X=X+1
9052  LET I=I+1
9054  IF I <= M THEN 9060
9056  LET Y=M2-X
9058  GOTO 9062
9060  NEXT K
9062  LET U=M-X
9064  LET V=N-Y
9066  LET Z=X*V-Y*U
9068  LET T=M1*(ABS(Z)-M1/2)^2
9070  LET C2=T/M/N/M2/(M1-M2)
9072  PRINT "TWO SAMPLE MEDIAN TEST."
9074  PRINT 
9076  PRINT "GROUP 1     ";X;U
9078  PRINT "GROUP 2     ";Y;V
9080  PRINT 
9082  PRINT "CHI-SQUARE = ";C2
9084  STOP 
9086  FOR I=1 TO U-L
9088  LET X=A[L]
9090  LET Q=L
9092  FOR J=L+1 TO U-I+1
9094  IF X >= A[J] THEN 9100
9096  LET X=A[J]
9098  LET Q=J
9100  NEXT J
9102  LET J=J-1
9104  IF J=Q THEN 9110
9106  LET A[Q]=A[J]
9108  LET A[J]=X
9110  NEXT I
9112  RETURN 
9900  DATA 4,6
9901  DATA 160,160,140,190
9902  DATA 117,145,147,120,150,120
9999  END 
