40 REM     ANVAR3:  ANALYSIS OF VARIANCE FOR TWO VARIABLES OF 
50 REM              CLASSIFICATION FACTORIAL DESIGN
60 REM
80 PRINT:PRINT
150 T,W,E,D=0
160 DIM X(20,20)
165 PRINT 'ENTER THE NUMBER OF ROWS, FOLLOWED BY THE NUMBER OF COLUMNS'
170 INPUT R,C
180 FOR I=1 TO R
190 U=0
200 FOR J=1 TO C
205 PRINT 'ENTER THE OBSERVATION IN ROW ';I;' AND COLUMN ';J
210 INPUT X(I,J)
220 U=U+X(I,J)
230 W=W+(X(I,J)^2)
240 NEXT J
250 E=E+(U^2)
260 T=T+U
270 NEXT I
280 FOR J=1 TO C
290 Q=0
300 FOR I=1 TO R
310 Q=Q+X(I,J)
320 NEXT I
330 D=D+(Q^2)
340 NEXT J
350 G=(T^2)/(R*C)
360 S1=W-G
370 S2=(E/C)-G
380 S3=(D/R)-G
390 S4=S1-S2-S3
400 M2=S2/(R-1)
410 M3=S3/(C-1)
420 A=(R*C)-1
430 B=A-R-C+2
440 M4=S4/B
450 F2=M2/M4
460 F3=M3/M4
470 PRINT
480 PRINT TAB(30);'ANOVA TABLE'
490 PRINT 'SOURCE','SUM SQ','DEG FREE','MEAN SQ','F RATIO'
500 PRINT 'ROWS',S2,R-1,M2,F2
510 PRINT 'COLS',S3,C-1,M3,F3
520 PRINT 'RESID',S4,B,M4
530 PRINT 'TOTAL',S1,A
540 PRINT
550 PRINT
560 F=F2
570 M=R-1
580 N=B
590 N1=C-1
600 GOSUB 620
602 F=F3
604 M=N1
606 GOSUB 620
610 GOTO 1120
620 A=2*INT(M/2)-M+2
630 B=2*INT(N/2)-N+2
640 W=F*M/N
650 Z=1/(1+W)
660 IF A<>1 THEN 760
670 IF B<>1 THEN 730
680 P=SQR(W*Z)
690 C=.313831
700 D=C*Z/P
710 P=2*C*ATN(P)
720 GOTO 830
730 P=SQR(W*Z)
740 D=.5*P*Z/W
750 GOTO 830
760 IF B<>1 THEN 810
770 P=SQR(Z)
780 D=.5*Z*P
790 P=1-P
800 GOTO 830
810 D=Z^2
820 P=W*Z
830 C=2*W/Z
840 FOR J=B+2 TO N STEP 2
850 D=(1+A/(J-2))*D*Z
860 IF A=1 THEN 890
870 P=(P+W)*Z
880 GOTO 900
890 P=P+D*C/(J-1)
900 NEXT J
910 C=W*Z
920 Z=2/Z
930 B=N-2
940 FOR I=A+2 TO M STEP 2
950 J=I+B
960 D=C*D*J/(I-2)
970 P=P-Z*D/J:NEXT I
980 IF 0<P AND P<1 THEN P1=P ELSE P1=0
990 IF P>1 THEN P2=1 ELSE P2=0
1000 IF P1>P2 THEN P=P1 ELSE P=P2
1010 P=1-P: PRINT
1020 IF P=1 AND F>1 THEN 1040
1030 IF P>.001 THEN 1050
1040 P=0
1050 PRINT 'PROBABILITY OF F>=';F;'WITH';M;'AND';N;'D.F. IS';P
1060 PRINT
1070 RETURN
1080 DATA 3,4
1090 DATA 7,6,8,7
1100 DATA 2,4,4,4
1110 DATA 4,6,5,3
1120 END