100 REM **** ANVAR1 12/70 **** STATISTICS PROGRAM **** 110 REM ANALYSIS OF VARIANCE 120 MAT T=ZER 130 MAT S=ZER 140 U,R,V,A=0 150 DIM X(50,20),N(20),T(20),S(20) 160 READ M 170 MAT READ N(M) 180 FOR J=1 TO M 190 FOR I=1 TO N(J) 200 READ X(I,J) 210 NEXT I 220 LET A=A+N(J) 230 NEXT J 240 FOR J=1 TO M 250 FOR I=1 TO N(J) 260 LET T(J)=T(J)+X(I,J) 270 LET S(J)=S(J)+X(I,J)*X(I,J) 280 NEXT I 290 LET U=U+T(J) 300 LET R=R+S(J) 310 LET V=V+T(J)*T(J)/N(J) 320 NEXT J 330 LET C=U*U/A 340 LET W=V-C 350 LET E=R-V 360 FOR I=1 TO 8 370 PRINT 380 NEXT I 390 PRINT "ANALYSIS OF VARIANCE TABLE" 400 PRINT 410 PRINT "GRAND TOTAL="U"N0. OBS.="A"MEAN="U/A 420 PRINT 430 PRINT 440 PRINT "SOURCE"," SS"," DF"," MS" 450 PRINT 460 PRINT"TREATMENTS",W,M-1,W/(M-1) 470 PRINT "ERROR",E,A-M,E/(A-M) 480 PRINT 490 PRINT "TOTAL",R-C,A-1 500 PRINT 510 PRINT 520 PRINT 530 LET F=(W/(M-1))/(E/(A-M)) 540 PRINT "F ="F"ON"M-1"AND"A-M"DEGREES OF FREEDOM." 550 LET G=F 560 LET N=A-M 570 LET M=M-1 580 GO SUB 630 590 FOR I=1 TO 8 600 PRINT 610 NEXT I 620 GOTO 990 630 REM SUBROUTINE TO COMPUTE F PROBABILITIES. 640 LET P=1 650 IF G<1 THEN 700 660 LET A=M 670 LET B=N 680 LET F=G 690 GO TO 730 700 LET A=N 710 LET B=M 720 LET F=1/G 730 LET A1=2/(9*A) 740 LET B1=2/(9*B) 750 LET Z=ABS((1-B1)*F^(.333333)-1+A1) 760 LET Z=Z/SQR(B1*F^(.666667)+A1) 770 IF B<4 THEN 810 780 LET P=(1+Z*(.196854+Z*(.115194+Z*(.000344+Z*.019527))))^4 790 LET P=.5/P 800 GO TO 830 810 LET Z=Z*(1+.08*Z^4/B^3) 820 GO TO 780 830 IF G<1 THEN 850 840 GOTO 870 850 LET P=1-P 860 GO TO 870 870 PRINT 880 LET P=INT(100000.*P)/100000. 890 PRINT"EXACT PROB. OF F="G;" WITH ("M;",";N") D.F. IS ";P 900 PRINT 910 RETURN 920 DATA 5 930 DATA 2,6,11,4,2 940 DATA 83,85 950 DATA 84,85,86,86,87,86 960 DATA 87,87,87,88,88,88,88,88,85,88,90 970 DATA 89,90,90,91 980 DATA 90,92 990 END