10  COM X[1000]
20  COM L[6],B[6],F[6],P[6],T$[5],N4,T,W,O,N7,W9,N1,N2,L1,Q9,T1
30  DIM A$[6]
40  FILES *
45  REM ANOVA:FACTORIAL ANALYSIS OF VARIANCE, HP 36870A, 6/74
50  PRINT 
60  PRINT "A N A L Y S I S   O F   V A R I A N C E   P R O G R A M"
70  PRINT 
80  PRINT "2000 SERIES VERSION: MODIFIED ON 6/30/76"
90  PRINT 
100  PRINT 
110  PRINT "DO YOU WANT INSTRUCTIONS (1=YES, 0=NO)";
120  INPUT Q8
130  IF Q8=0 THEN 480
140  PRINT 
150  PRINT "*** INSTRUCTIONS ***"
160  PRINT "THIS PROGRAM COMPUTES UP TO A FIVE-WAY FACTORIAL ANOVA"
170  PRINT "WITH A MAX. OF 1000 SUBJECTS IN THE DESIGN. THE PROGRAM WILL"
180  PRINT "WORK FOR ANY BALANCED (EQUAL # OF SUBJ. PER CELL) DESIGN"
190  PRINT "THAT HAS AT LEAST ONE SUBJECT PER CELL. SINCE NO F-RATIOS"
200  PRINT "ARE PRINTED, ONE MUST CALCULATE THEM FROM THE SUMMARY TABLE."
210  PRINT "DATA MAY BE ENTERED IN DATA STATEMENTS BEGINNING ON LINE 5000"
220  PRINT "OR FROM DATA FILES STORED ON DISC."
240  PRINT "ENTER ";
250  PRINT "DATA SO THAT SUBJECTS ARE INCREMENTED FIRST (IF MORE THAN"
260  PRINT "ONE PER CELL) AND VARIABLE 'A' IS INCREMENTED LAST."
270  PRINT "FOR EXAMPLE, IN A 2 X 2 DESIGN WITH TWO SUBJECTS PER CELL,"
280  PRINT "THE 8 DATA POINTS SHOULD BE ENTERED IN THIS ORDER:"
290  PRINT 
300  PRINT TAB(5);"A     B     S"
310  PRINT TAB(5);"1     1     1"
320  PRINT TAB(5);"1     1     2"
330  PRINT TAB(5);"1     2     1"
340  PRINT TAB(5);"1     2     2"
350  PRINT TAB(5);"2     1     1"
360  PRINT TAB(5);"2     1     2"
370  PRINT TAB(5);"2     2     1"
380  PRINT TAB(5);"2     2     2"
390  PRINT 
400  PRINT "RUNNING THE PROGRAM DESTROYS DATA IN DATA STATEMENTS. IF YOU"
410  PRINT "WANT TO SAVE YOUR DATA, PUNCH ON PAPER TAPE BY TYPING PUN-5000"
420  PRINT "AND TURNING ON THE TAPE PUNCH BEFORE RUNNING. IF MORE THAN ONE"
430  PRINT "PROBLEM IS TO BE RUN AT A SINGLE TERMINAL SESSION, IT WILL BE"
440  PRINT " NECESSARY TO TYPE GET-ANOVA BEFORE ENTERING DATA FOR ADDITIONAL"
450  PRINT " PROBLEMS. GET-ANOVA MUST ALSO BE TYPED BEFORE RUNNING ADDITIONAL"
455  PRINT " PROBLEMS USING DATA FILES. NOW GET ANOVA, ENTER YOUR DATA,"
457  PRINT "AND RUN."
460  STOP 
470  REM:THIS DETERMINES THE SIZE OF THE ANALYSIS
480  LET T$="ABCDE"
490  LET T=W=O=N7=W9=0
500  LET N1=N2=L1=1
505  PRINT 
510  PRINT "1= DATA ON FILE, 0= DATA IN DATA STATEMENTS. WHICH";
520  INPUT Z9
530  IF Z9=0 THEN 570
540  PRINT "FILE NAME";
550  INPUT A$
560  ASSIGN A$,1,Z8
565  PRINT 
570  PRINT "NUMBER OF VARIABLES";
580  INPUT N4
590  LET N4=N4+1
600  IF N4>6 THEN 950
610  PRINT "NUMBER OF REPLICATES (# OF SUBJ. PER CELL)";
620  INPUT L[N4]
630  FOR I=1 TO N4
640  IF I=N4 THEN 670
650  PRINT "# OF LEVELS FOR VARIABLE ";T$[I,I];
660  INPUT L[I]
670  LET B[I]=F[I]=P[I]=0
680  LET N1=N1*2
690  LET N2=N2*L[I]
700  NEXT I
710  PRINT 
720  PRINT "DO YOU WANT THE MEANS & SUMS OF SQUARES PRINTED FOR"
730  PRINT "POST-HOC COMPARISONS (1=YES, 0=NO)";
740  INPUT Q9
750  PRINT 
760  PRINT 
770  PRINT 
780  FOR I=1 TO N4
790  LET L1=L1*L[I]
800  LET P[I]=N2/L1
810  NEXT I
820  REM: READS DATA AND STORES IN CELLS "X"
830  LET O=N2
840  LET T1=0
850  FOR I=1 TO N2
860  IF Z9=0 THEN 890
870  READ #1;X[I]
880  GOTO 900
890  READ X[I]
900  LET T=T+X[I]^2
910  REM: DETERMINES SUM OF ALL CELLS
920  LET T1=T1+X[I]
930  NEXT I
940  CHAIN "TANOVA"
950  PRINT 
960  PRINT "YOUR ANALYSIS CONTAINS TOO MANY VARIABLES FOR THIS"
970  PRINT "PROGRAM. SEE THE COMPUTER ASSISTANT FOR USE OF THE"
980  PRINT "BIO-MED ANOVA PROGAM."
990  STOP 
9999  END 
