10  COM N[64],T[64],G[64],J[64],N2,N6,N7,N8,A1,B1,C1,T1
20  COM R[64],W[64],D[10]
30  REM ***** AVAR23, VERSION 1, AUGUST 5, 1971        R.J.M. TAYLOR
40  REM ***** TRANSLATED FROM 'FORTRAN PROGRAMMING FOR THE BEHAVIORAL
50  REM ***** SCIENCES ' BY D.J. VELDMAN
60  DIM X[64],K[64]
70  FILES SX,G
80  REM ENTER FILES STATEMENT HERE
90  N5=10
100  PRINT "NUMBER OF DEPENDENT VARIABLES";
110  INPUT N2
120  MAT J=ZER[N2]
130  MAT G=ZER[N2]
140  PRINT "NUMBER OF LEVELS FOR THE A FACTOR";
150  INPUT N6
160  PRINT "NUMBER OF LEVELS FOR THE B FACTOR";
170  INPUT N7
180  PRINT "NUMBER OF LEVELS FOR THE C FACTOR";
190  INPUT N8
200  PRINT "NUMBER OF SUBJECTS PER ABC CELL";
210  INPUT N1
220  N9=N6*N7*N8
230  PRINT "MISSING DATA CODED WITH WHAT NUMBER";
240  INPUT M1
250  FOR I=1 TO N2
260  T[I]=R[I]=W[I]=0
270  NEXT I
280  MAT N=CON
290  N1=ABS(N1)
300  J9=0
310  FOR M=1 TO N9
320  N=N1
330  IF N#9999 THEN 360
340  PRINT "NUMBER OF SUBJECTS IN GROUP";M;
350  INPUT N
360  FOR I=1 TO N2
370  J[I]=K[I]=0
380  G[I]=N
390  NEXT I
400  FOR I=1 TO N
410  F5=INT(((J9+I)*2-1)/128)+3
420  I5=(J9+I)-64*(F5-3)
430  READ #F5,I5*2-1
440  MAT  READ #F5;X[N2]
450  FOR J=1 TO N2
460  IF X[J]#M1 THEN 490
470  G[J]=G[J]-1
480  X[J]=0
490  J[J]=J[J]+X[J]
500  K[J]=K[J]+X[J]^2
510  NEXT J
520  NEXT I
530  J9=J9+N
540  FOR I=1 TO N2
550  IF G[I]>0 THEN 580
560  N[I]=2
570  GOTO 610
580  W[I]=W[I]+(K[I]-J[I]^2/G[I])
590  R[I]=R[I]+1/G[I]
600  J[I]=J[I]/G[I]
610  T[I]=T[I]+G[I]
620  NEXT I
630  R=M*2-1
640  READ #1,R
650  MAT  PRINT #1;J
660  READ #2,R
670  MAT  PRINT #2;G
680  NEXT M
690  REM WRITE(2)  SX,G
700  T1=N9
710  A1=N6
720  B1=N7
730  C1=N8
740  D[2]=T1-1
750  D[3]=A1-1
760  D[4]=B1-1
770  D[5]=C1-1
780  D[6]=D[3]*D[4]
790  D[7]=D[3]*D[5]
800  D[8]=D[4]*D[5]
810  D[9]=D[3]*D[8]
820  CHAIN "$AVR23+"
830  END 
