9000  REM  STAT16--ANALYSIS OF VARIANCE, GRAECO-LATIN SQUARES
9001  REM HP CONTRIBUTED LIBRARY, REV. B, 2/75
9002  PRINT "ITEM";TAB(10);"SUM-SQR";TAB(26);"DEG. FREE.";
9004  PRINT TAB(44);"MEAN-SQR";TAB(59);"F-RATIO"
9006  PRINT "----";TAB(10);"-------";TAB(26);"---------";
9008  PRINT TAB(44);"--------"TAB(59);"-------"
9010  READ N
9012  FOR I=1 TO N
9014  FOR J=1 TO N
9016  READ M[I,J]
9018  NEXT J
9020  NEXT I
9022  FOR I=1 TO N
9024  FOR J=1 TO N
9026  READ N[I,J]
9028  NEXT J
9030  NEXT I
9032  LET S=0
9034  LET S0=0
9036  FOR I=1 TO N
9038  LET R[I]=0
9040  LET C[I]=0
9042  LET T[I]=0
9044  LET G[I]=0
9046  NEXT I
9048  FOR I=1 TO N
9050  FOR J=1 TO N
9052  READ X
9054  LET R[I]=R[I]+X
9056  LET C[J]=C[J]+X
9058  LET T[M[I,J]]=T[M[I,J]]+X
9060  LET G[N[I,J]]=G[N[I,J]]+X
9062  LET S=S+X
9064  LET S0=S0+X*X
9066  NEXT J
9068  NEXT I
9070  LET S3=0
9072  LET S4=0
9074  LET S5=0
9076  LET S6=0
9078  FOR I=1 TO N
9080  LET S4=S4+R[I]^2
9082  LET S5=S5+C[I]^2
9084  LET S6=S6+T[I]^2
9086  LET S3=S3+G[I]^2
9088  NEXT I
9090  LET C=S*S/N/N
9092  LET S4=S4/N-C
9094  LET S5=S5/N-C
9096  LET S6=S6/N-C
9098  LET S3=S3/N-C
9100  LET D1=N-1
9102  LET D2=D1*(N-3)
9104  LET S7=S0-C-S4-S5-S6
9106  LET S7=S7-S3
9108  LET S8=S7/D2
9110  PRINT "ROWS";TAB(9);S4;TAB(30);D1;TAB(43);S4/D1;TAB(58);S4/D1/S8
9112  PRINT "COLS";TAB(9);S5;TAB(30);D1;TAB(43);S5/D1;TAB(58);S5/D1/S8
9114  PRINT "TREAT L";TAB(9);S6;TAB(30);D1;TAB(43);S6/D1;TAB(58);S6/D1/S8
9116  PRINT "TREAT G";TAB(9);S3;TAB(30);D1;TAB(43);S3/D1;TAB(58);S3/D1/S8
9118  PRINT "ERROR";TAB(9);S7;TAB(30);D2;TAB(43);S8
9120  STOP 
9900  DATA 4
9901  DATA 1,2,3,4,4,1,2,3,3,4,1,2,2,3,4,1
9902  DATA 4,3,2,1,3,2,1,4,2,1,4,3,1,4,3,2
9903  DATA 24,47,35,42
9904  DATA 47,85,23,47
9905  DATA 65,49,23,62
9906  DATA 12,14,19,23
9999  END 
