100 REM **** ANVAR2 **** STATISTICS PROGRAM ***
110 REM **** VERSION 1 ***** 1/69 ****
120 REM ANALYSIS OF VARIANCE
130 DIM C(10),T(10)
140 READ N
150 MAT READ M(N,N)
160 MAT R=ZER
170 MAT C=ZER
180 MAT T=ZER
190 LET S,S0=0
200 FOR I=1 TO N
210 FOR J=1 TON
220 READ X
230 LET R(I)=R(I)+X
240 LET C(J)=C(J)+X
250 LET T(M(I,J))=T(M(I,J))+X
260 LET S=S+X
270 LET S0=S0+X*X
280 NEXT J
290 NEXT I
300 LET S4,S5,S6=0
310 FOR I=1 TO N
320 LET S4=S4+R(I)^2
330 LET S5=S5+C(I)^2
340 LET S6=S6+T(I)^2
350 NEXT I
360 LET C=S*S/N/N
370 LET S4=S4/N-C
380 LET S5=S5/N-C
390 LET S6=S6/N-C
400 LET D1=N-1
410 LET D2=D1*(N-2)
420 LET S7=S0-C-S4-S5-S6
430 PRINT "ITEM  ","SUM-SQR  ","DEG.FREE.  ","MEAN-SQR  ","F-RATIO"
440 PRINT "----  ","------  ","-------  ","-------  ","-------"
450 LET S8=S7/D2
460 PRINT "ROWS  ",S4,D1,S4/D1,S4/D1/S8
470 PRINT "COLS  ",S5,D1,S5/D1,S5/D1/S8
480 PRINT "TREATS",S6,D1,S6/D1,S6/D1/S8
490 PRINT "ERROR",S7,D2,S8
500 LET G=S4/D1/S8
510 LET M=D1
520 LET N=D2
530 GOSUB 590
540 LET G=S5/D1/S8
550 GOSUB 590
560 LET G=S6/D1/S8
570 GOSUB 590
580 GOTO 930
590 LET P=1
600 IF G<1 THEN 650
610 LET A=M
620 LET B=N
630 LET F=G
640 GO TO 680

650 LET A=N
660 LET B=M
670 LET F=1/G
680 LET A1=2/(9*A)
690 LET B1=2/(9*B)
700 LET Z=ABS((1-B1)*F^(.333333)-1+A1)
710 LET Z=Z/SQR(B1*F^(.666667)+A1)
720 IF B<4 THEN 760
730 LET P=(1+Z*(.196854+Z*(.115194+Z*(.000344+Z*.019527))))^4
740 LET P=.5/P
750 GOTO780
760 LET Z=Z*(1+.08*Z^4/B^3)
770 GOTO730
780 IF G<1 THEN 800
790 GOTO820
800 LET P=1-P
810 GOTO820
820 PRINT
830 LET P=INT(100000.*P)/100000.
840 PRINT "PROBABILITY OF F>=";G;"WITH";M;"AND";N;"D.F. IS";P
850 PRINT
860 RETURN
870 DATA 4
880 DATA 1,2,3,4,4,1,2,3,3,4,1,2,2,3,4,1
890 DATA 81,41,44,53
900 DATA 38,97,42,49
910 DATA 31,43,67,36
920 DATA 57,33,43,81
930 END