9000 REM **** STAT3 ****** STATISTICS PROGRAM
9001 REM ***** REVISED BY LINDA RUFF ***** (12/18/72)
9002 REM  SPEARMAN RANK CORRELATION
9003 DIM A(100),B(100),C(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 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 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 C(I)=A(I)
9028 NEXT I
9029 FOR K=1 TO N
9030 FOR I=1 TO N-1
9031 IF C(I)>C(I+1) THEN A=C(I):C(I)=C(I+1):C(I+1)=A
9032 NEXT I
9033 NEXT K
9034 FOR J=1 TO N: FOR I=1 TO N
9035 IF A(I)=C(J) THEN A(I)=J
9036 NEXT I
9037 NEXT J
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 GO TO 9009
9047 REM NOW WE ATART 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 GOTO9062
9061 LET R=(N*S3-S)/SQR(D)
9062 PRINT "SPEARMAN RANK CORRELATION COEFFICIENT"
9063 PRINT "  R = ";R
9064 GO TO 9999
9900 DATA 5,480,56,500,61,520,78,540,71,560,82
9999 END