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