1  REM ****  HP BASIC PROGRAM LIBRARY  ******************************
2  REM
3  REM       GANOVA: ANALYSIS OF VARIANCE (2-WAY)
4  REM
5  REM       36501 REV A
6  REM
7  REM ****  CONTRIBUTED PROGRAM  ***********************************
1000  DIM C[20,20],D[20],G[20]
1010  PRINT 
1011  MAT C=ZER
1012  MAT D=ZER
1013  MAT G=ZER
1014  A1=D1=G1=H1=T1=0
1020  PRINT "HOW MANY OBSERVATIONS PER CELL DO YOU HAVE";
1030  INPUT N1
1040  PRINT "HOW MANY COLUMNS DO YOU HAVE";
1050  INPUT K1
1060  PRINT "HOW MANY ROWS DO YOU HAVE";
1070  INPUT J1
1080  N2=N1*K1*J1
1090  PRINT 
1100  PRINT 
1110  FOR I=1 TO K1
1120  FOR J=1 TO J1
1130  FOR L=1 TO N1
1140  READ X1
1150  C[I,J]=C[I,J]+X1
1160  A1=A1+X1*X1
1170  NEXT L
1180  H1=H1+C[I,J]*C[I,J]
1190  T1=T1+C[I,J]
1200  D[J]=D[J]+C[I,J]
1210  G[I]=G[I]+C[I,J]
1220  NEXT J
1230  G1=G[I]*G[I]+G1
1240  NEXT I
1250  S1=G1/(N1*J1)-(T1*T1)/N2
1260  S2=A1-H1/N1
1270  S3=A1-(T1*T1)/N2
1280  FOR I=1 TO J1
1290  D1=D1+D[I]*D[I]
1300  NEXT I
1310  S4=D1/(N1*K1)-(T1*T1)/N2
1320  S5=S3-S1-S2-S4
1330  PRINT " "," ","***  ANOVA TABLE  ***"
1340  PRINT "SOURCE","DF","SUM OF SQ","VARIANCE","F RATIO"
1350  PRINT "ROW",J1-1,S4,S4/(J1-1),(S4/(J1-1))/(S2/((N1-1)*K1*J1))
1360  PRINT "COLUMN",K1-1,S1,S1/(K1-1),(S1/(K1-1))/(S2/((N1-1)*J1*K1))
1370  PRINT "INTERACTION",(K1-1)*(J1-1),S5,S5/((K1-1)*(J1-1)),
1380  PRINT (S5/((K1-1)*(J1-1)))/(S2/((N1-1)*K1*J1))
1390  PRINT "ERROR",(N1-1)*K1*J1,S2,S2/((N1-1)*K1*J1)
1400  PRINT "TOTAL",N2-1,S3
1410  PRINT 
1420  STOP 
9990  END 
