1  REM  ****  HP BASIC PROGRAM LIBRARY  *******************************
2  REM
3  REM        STAT18:   COMPUTES AN ANALYSIS OF VARIANCE TABLE AND 
4  REM                  F-RATIO FOR TREATMENTS FOR A YOUDEN SQUARE 
5  REM                  DESIGN
6  REM        36727 REV  A   10/73
7  REM
8  REM  ****  CONTRIBUTED PROGRAM  ************************************
9000  REM  STAT18--ANALYSIS OF VARIANCE, YOUDEN SQUARE
9002  PRINT "ANALYSIS OF VARIANCE TABLE"
9004  PRINT 
9006  PRINT "ITEM","SUM-SQR","DEG. FREE.","MEAN-SQR"
9008  PRINT 
9010  READ N,K
9012  LET L=K*(K-1)/(N-1)
9014  FOR I=1 TO N
9016  FOR H=1 TO K
9018  READ M[I,H]
9020  LET C[H]=0
9022  NEXT H
9024  LET R[I]=0
9026  LET T[I]=0
9028  LET P[I]=0
9030  LET Q[I]=0
9032  NEXT I
9034  FOR I=1 TO N
9036  FOR J=1 TO N
9038  READ N[I,J]
9040  NEXT J
9042  NEXT I
9044  LET S=0
9046  LET S2=0
9048  FOR I=1 TO N
9050  FOR H=1 TO K
9052  READ X[I,H]
9054  LET R[I]=R[I]+X[I,H]
9056  LET C[H]=C[H]+X[I,H]
9058  LET T[M[I,H]]=T[M[I,H]]+X[I,H]
9060  LET S=S+X[I,H]
9062  LET S2=S2+X[I,H]^2
9064  NEXT H
9066  NEXT I
9068  FOR J=1 TO N
9070  FOR I=1 TO N
9072  LET P[J]=P[J]+N[I,J]*R[I]
9074  NEXT I
9076  NEXT J
9078  FOR J=1 TO N
9080  LET Q[J]=K*T[J]-P[J]
9082  NEXT J
9084  LET C=S*S/N/K
9086  LET R1=0
9088  LET C1=0
9090  LET T1=0
9092  FOR I=1 TO N
9094  LET R1=R1+R[I]^2
9096  LET T1=T1+Q[I]^2
9098  NEXT I
9100  FOR H=1 TO K
9102  LET C1=C1+C[H]^2
9104  NEXT H
9106  LET R2=R1/K-C
9108  LET C2=C1/N-C
9110  LET T2=T1/N/K/L
9112  LET D=(N-1)*(K-2)
9114  LET D1=N*(K-2)+1
9116  LET E2=S2-R2-C2-T2-C
9118  PRINT "GRAND TOTAL",S2,N*K
9120  PRINT "GRAND MEAN",C," 1"
9122  PRINT "TREATMENTS",T2,N-1,T2/(N-1)
9124  PRINT "ROWS",R2,N-1"     ...SS NOT ADJUSTED..."
9126  PRINT "COLUMNS",C2,K-1,C2/(K-1)
9128  PRINT "ERROR",E2,D,E2/D
9130  PRINT 
9132  PRINT 
9134  LET F=T2/(N-1)/(E2/D)
9136  PRINT "TREATMENT F-RATIO ="F", ON"N-1"AND"D"DEGREES OF FREEDOM."
9138  PRINT 
9140  PRINT "IF MSC/MSE ="C2/(K-1)/(E2/D)"IS NOT SIGNIFICANT, IT MAY BE"
9142  PRINT "DESIRABLE TO POOL COLUMN AND ERROR SS TO OBTAIN AS AN"
9144  PRINT "ERROR MS ESTIMATE"(C2+E2)/D1"WITH"D1"DEGREES OF FREEDOM."
9900  DATA 4,3
9901  DATA 1,2,3
9902  DATA 4,1,2
9903  DATA 2,3,4
9904  DATA 3,4,1
9905  DATA 1,1,1,0
9906  DATA 1,1,0,1
9907  DATA 0,1,1,1
9908  DATA 1,0,1,1
9909  DATA 2,1,0
9910  DATA -2,2,2
9911  DATA -1,-1,-3
9912  DATA 0,-4,2
9999  END 
