3 REM GANOVA: ANALYSIS OF VARIANCE (2-WAY) 4 REM 5 REM TRANLATED TO BASIC-PLUS BY DSAA 7 REM 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 1101 PRINT'ENTER YOUR OBSERVATIONS IN THE FOLLOWING MANNER. ENTER ALL' 1102 PRINT'OBSERVATIONS IN THE FIRST CELL OF THE FIRST COLUMN, THEN' 1103PRINT'THE OBSERVATIONS OF THE SECOND CELL OF THE FIRST COLUMN,ETC.' 1110 FOR I=1 TO K1 1120 FOR J=1 TO J1 1130 FOR L=1 TO N1 1140 INPUT 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 NEXTI 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 1321 PRINT:PRINT 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 3305 DATA 23,25,18,26,18,26,20,19,22,17 9990 END