8990  REM  ***  HP TIME-SHARED BASIC PROGRAM LIBRARY  ********************
8991  REM  
8992  REM          ANVAR1:  36039  (A410)  REV B -- 9/71
8993  REM          ANALYSIS OF VARIANCE FOR A ONE-WAY RANDOMIZED DESIGN
8994  REM
8995  REM  ***  CONTRIBUTED PROGRAM  **************************************
9001  REM  ANALYSIS OF VARIANCE
9002  MAT T=ZER
9003  MAT S=ZER
9004  U=R=V=A=0
9005  DIM X[50,20],N[20],T[20],S[20]
9006  READ M
9007  MAT  READ N[M]
9008  FOR J=1 TO M
9009  FOR I=1 TO N[J]
9010  READ X[I,J]
9011  NEXT I
9012  LET A=A+N[J]
9013  NEXT J
9014  FOR J=1 TO M
9015  FOR I=1 TO N[J]
9016  LET T[J]=T[J]+X[I,J]
9017  LET S[J]=S[J]+X[I,J]*X[I,J]
9018  NEXT I
9019  LET U=U+T[J]
9020  LET R=R+S[J]
9021  LET V=V+T[J]*T[J]/N[J]
9022  NEXT J
9023  LET C=U*U/A
9024  LET W=V-C
9025  LET E=R-V
9026  FOR I=1 TO 8
9027  PRINT 
9028  NEXT I
9029  PRINT "ANALYSIS OF VARIANCE TABLE"
9030  PRINT 
9031  PRINT "GRAND TOTAL="U"NO. OBS.="A"MEAN="U/A
9032  PRINT 
9033  PRINT 
9034  PRINT "SOURCE","   SS","   DF","   MS"
9035  PRINT 
9036  PRINT "TREATMENTS",W,M-1,W/(M-1)
9037  PRINT "ERROR",E,A-M,E/(A-M)
9038  PRINT 
9039  PRINT "TOTAL",R-C,A-1
9040  PRINT 
9041  PRINT 
9042  PRINT 
9043  LET F=(W/(M-1))/(E/(A-M))
9044  PRINT "F ="F"ON"M-1"AND"A-M"DEGREES OF FREEDOM."
9045  LET N=A-M
9046  LET M=M-1
9047  A=2*INT(M/2)-M+2
9048  B=2*INT(N/2)-N+2
9049  W=F*M/N
9050  Z=1/(1+W)
9051  IF A#1 THEN 9061
9052  IF B#1 THEN 9058
9053  P=SQR(W)
9054  K=.31831
9055  D=K*Z/P
9056  P=2*K*ATN(P)
9057  GOTO 9068
9058  P=SQR(W*Z)
9059  D=.5*P*Z/W
9060  GOTO 9068
9061  IF B#1 THEN 9066
9062  P=SQR(Z)
9063  D=.5*Z*P
9064  P=1-P
9065  GOTO 9068
9066  D=Z^2
9067  P=W*Z
9068  K=2*W/Z
9069  FOR J=B+2 TO N STEP 2
9070  D=(1+A/(J-2))*D*Z
9071  IF A=1 THEN 9074
9072  P=(P+W)*Z
9073  GOTO 9075
9074  P=P+D*K/(J-1)
9075  NEXT J
9076  K=W*Z
9077  Z=2/Z
9078  B=N-2
9079  FOR I=A+2 TO M STEP 2
9080  J=I+B
9081  D=K*D*J/(I-2)
9082  P=P-Z*D/J
9083  NEXT I
9084  P=(P>1) MAX (0<P AND P<1)*P
9085  P=1-P
9086  IF P=1 AND F>1 THEN 9088
9087  IF P >= .001 THEN 9089
9088  P=0
9089  PRINT "PROBABILITY OF F>=";F;"WITH";M;"AND";N;"D.F. IS";P
9900  DATA 5
9901  DATA 2,6,11,4,2
9902  DATA 83,85
9903  DATA 84,85,86,86,87,86
9904  DATA 87,87,87,88,88,88,88,88,85,88,90
9905  DATA 89,90,90,91
9906  DATA 90,92
9999  END 
