40 REM        ANALYSIS OF VARIANCE FOR A TWO-WAY EXPERIMENT
50 REM        WITH REPEATED OBSERVATIONS
60 REM		TRANSLATED TO BASIC-PLUS BY DSAA
70 REM
80 DIM X(20,40),Y(20,20),D(20),E(20)
82 PRINT 'PLEASE ENTER THE NUMBER OF ROWS, FOLLOWED BY THE NUMBER OF'
84 PRINT 'COLUMNS, THEN THE NUMBER OF REPLICATIONS PER BOX.'
90 INPUT R,C,P
100 T2,T3=0
110 FOR I=1 TO R
120 T4=0
130 FOR J=1 TO C
140 T1=0
145 PRINT ' ENTER THE OBSERVATIONS OF BLOCK ('I;',';J;')'
150 FOR K=1 TO P
160 INPUT X(J,K)
170 T1=T1+X(J,K)
180 T2=T2+X(J,K)^2
190 NEXT K
200 Y(I,J)=T1
210 T3=T3+T1
220 T4=T4+T1
230 NEXT J
240 E(I)=T4
250 NEXT I
260 T6=0
270 FOR J=1 TO C
280 T5=0
290 FOR I=1 TO R
300 T5=T5+Y(I,J)
310 T6=T6+Y(I,J)^2
320 NEXT I
330 D(J)=T5
340 NEXT J
350 H=T3^2/(R*C*P)
360 S1=T2-H
370 S2=(T6/P)-H
380 S3=S1-S2
390 N1=(R*C*P)-1
400 N2=(R*C)-1
410 N3=N1-N2
420 M2=S2/N2
430 M3=S3/N3
440 F2=M2/M3
445 PRINT:PRINT
450 PRINT TAB (30);'ANOVA TABLE I'
460 PRINT 'SOURCE','SUM','DEG FREE','MEAN SQ','F RATIO'
470 PRINT 'MEANS',S2,N2,M2,F2
480 PRINT 'WITHIN',S3,N3,M3
490 PRINT 'TOTAL',S1,N1
500 R3=P
510 F=F2
520 M=N2
530 N=N3
540 GOSUB 1170
550 T7,T8=0
560 FOR I=1 TO R
570 T7=T7+E(I)^2
580 NEXT I
590 FOR J=1 TO C
600 T8=T8+D(J)^2
610 NEXT J
620 S4=(T7/(R3*C))-H
630 S5=(T8/(R3*R))-H
640 N4=R-1
650 N5=C-1
660 S6=S2-S4-S5
670 N6=N2-N4-N5
680 M4=S4/N4
690 M5=S5/N5
700 M6=S6/N6
710 F4=M4/M3
720 F5=M5/M3
730 F6=M6/M3
740 PRINT
750 PRINT TAB(29);'ANOVA TABLE II'
760 PRINT 'SOURCE','SUM','DEG FREE','MEAN SQ','F RATIO'
770 PRINT 'ROWS',S4,N4,M4,F4
780 PRINT 'COLS',S5,N5,M5,F5
790 PRINT 'INTER',S6,N6,M6,F6
800 PRINT 'SUBTOT',S2,N2
810 PRINT 'WITHIN',S3,N3,M3
820 PRINT 'TOTAL',S1,N1
830 F=F4
840 M=N4
850 GOSUB 1170
860 F=F5
870 M=N5
880 GOSUB 1170
890 F=F6
900 M=N6
910 GOSUB 1170
920 PRINT 'IF THE INTERACTION EFFECT IS NOT SIGNIFICANT AND '
930 PRINT 'IF YOU AWISH TO POOL INTERACTION AND WITHIN SUMS OF '
940 PRINT 'SQUARES TO FORM RESIDUAL SUM OF SQUARES TYPE THE'
950 PRINT 'NUMBER 1 OTHERWISE TYPE NUMBER 0.'
960 INPUT Q
970 IF Q<.5 THEN 1660
980 S3=S3+S6
990 N3=N3+N6
1000 M3=S3/N3
1010 F4=M4/M3
1020 F5=M5/M3
1030 PRINT TAB(38);'ANOVA TABLE III'
1040 PRINT 'SOURCE','SUM','DEG FREE','MEAN SQ','F RATIO'
1050 PRINT 'ROWS',S4,N4,M4,F4
1060 PRINT 'COLS',S5,N5,M5,F5
1070 PRINT 'RESID',S3,N3,M3
1080 PRINT 'TOTAL',S1,N1
1090 F=F4
1100 M=N4
1110 N=N3
1120 GOSUB 1170
1130 F=F5
1140 M=N5
1150 GOSUB 1170
1160 GOTO 1660
1170 A=2*INT(M/2)-M+2
1180 B=2*INT(N/2)-N+2
1190 W=F*M/N
1200 Z=1/(1+W)
1210 IF A<>1 THEN 1310
1220 IF B<>1 THEN 1280
1230 P=SQR(W)
1240 K=.31831
1250 D=K*Z/P
1260 P=2*K*ATN(P)
1270 GOTO 1380
1280 P=SQR(W*Z)
1290 D=.5*P*Z/W
1300 GOTO 1380
1310 IF B<>1 THEN 1360
1320 P=SQR(Z)
1330 D=.5*Z*P
1340 P=1-P
1350 GOTO 1380
1360 D=Z^2
1370 P=W*Z
1380 K=2*W/Z
1390 FOR J=B+2 TO N STEP 2
1400 D=(1+A/(J-2))*D*Z
1410 IF A=1 THEN 1440
1420 P=(P+W)*Z
1430 GOTO 1450
1440 P=P+D*K/(J-1)
1450 NEXT J
1460 K=W*Z
1470 Z=2/Z
1480 B=N-2
1490 FOR I=A+2 TO M STEP 2
1500 J=I+B
1510 D=K*D*J/(I-2)
1520 P=P-Z*D/J:NEXT I
1530 IF 0<P AND P<1 THEN P1=P ELSE P1=0
1540 IF P>1 THEN P2=1 ELSE P2=0
1550 IF P1>P2 THEN P=P1 ELSE P=P2
1560 P=1-P:IF P=1 AND F>1 THEN 1580
1570 IF P>.001 THEN 1590
1580 P=0
1590 PRINT
1600 PRINT 'PROBABILITY OF F>=';F;'WITH';M;'AND';N;'D.F. IS';P
1610 PRINT
1620 RETURN
1660 END