9000  REM  *****  STAT3  *****  STATISTICS PROGRAM  *****
9001  REM  *****  VERSION 1  *****  7/31/69  *****
9002  REM  SPEARMAN RANK CORRELATION
9003  DIM A[100],B[100]
9004  LET P=0
9005  READ N
9006  FOR I=1 TO N
9007  READ A[I],B[I]
9008  NEXT I
9009  FOR I=1 TO N-1
9010  LET X=A[1]
9011  LET Y=B[1]
9012  LET Q=1
9013  FOR J=2 TO N-I+1
9014  IF X >= A[J] THEN 9018
9015  LET X=A[J]
9016  LET Y=B[J]
9017  LET Q=J
9018  NEXT J
9019  LET J=J-1
9020  LET A[Q]=A[J]
9021  LET A[J]=X
9022  LET B[Q]=B[J]
9023  LET B[J]=Y
9024  NEXT I
9025  REM NOW REPLACE THE ELEMENT BY ITS AVERAGE RANK
9026  FOR I=1 TO N
9027  LET X=A[I]
9028  FOR J=I+1 TO N
9029  IF X <> A[J] THEN 9032
9030  NEXT J
9031  LET J=N+1
9032  LET Z=(J+I-1)/2
9033  FOR K=I TO J-1
9034  LET A[K]=Z
9035  NEXT K
9036  LET I=J-1
9037  NEXT I
9038  IF P=1 THEN 9047
9039  REM NOW WE INTERCHANGE A AND B AND SORT AGAIN.
9040  FOR I=1 TO N
9041  LET T=A[I]
9042  LET A[I]=B[I]
9043  LET B[I]=T
9044  NEXT I
9045  LET P=1
9046  GOTO 9009
9047  REM NOW WE START COMPUTING THE CORRELATION COEFFICIENT.
9048  LET S1=S2=S3=0
9049  FOR I=1 TO N
9050  LET S1=S1+A[I]*A[I]
9051  LET S2=S2+B[I]*B[I]
9052  LET S3=S3+A[I]*B[I]
9053  NEXT I
9054  LET S=N*N*(N+1)*(N+1)/4
9055  LET D1=N*S1-S
9056  LET D2=N*S2-S
9057  LET D=D1*D2
9058  IF D <> 0 THEN 9061
9059  LET R=0
9060  GOTO 9062
9061  LET R=(N*S3-S)/SQR(D)
9062  PRINT "SPEARMAN RANK CORRELATION COEFFICIENT"
9063  PRINT "  R = ";R
9064  STOP 
9900  DATA 5,480,56,500,61,520,78,540,71,560,82
9999  END 
