C-PS/8 FOCAL, 1971
 1.10 A !!,"ANALYSIS ",LI;I (LI-36)1.2,1.3,1.2
 1.20 L R SORT 2.15
 1.30 S S=0;L R STAT6 2.4;WE CANNOT USE SAVED MEANS
 2.10 S S=0;S Q=0
 2.20 F I=1,KB;S QC(I)=0
 2.30 F I=1,N;S QR(I)=0
 3.05 T !!,"ANOVA",!
 3.10 F I=1,N;F J=1,KB;DO 10
 3.20 F I=1,N;F J=1,KB;DO 11
 3.30 I (Q) 1.1,1.1,3.4
 3.40 S S1=0;F I=1,KB;S S1=S1+QC(I)*QC(I)/N
 3.50 S S2=0;F I=1,N;S S2=S2+QR(I)*QR(I)/KB
 3.60 S SQ=S*S/(KB*N)
 3.70 T !,%8.04,,,"SOURCE   SS       DF      MS         F
 3.80 S S1=S1-SQ;S S2=S2-SQ
 3.90 S SQ=Q-SQ; S Q=SQ-S1-S2
 4.05 S DF=N*KB-(N+KB)+1
 4.07 T !,"COL ",%8.04,S1,%3.00,KB-1,%8.04,S1/(KB-1)
 4.10 T S1*DF/(Q*(KB-1)),!,"ROWS",S2,%3.00,N-1,%8.04,S2/(N-1)
 4.20 T S2*DF/(Q*(N-1)),!,"ERR "
 4.30 T Q,%3.00,DF,%8.04,Q/DF,!,"TOT ",SQ,%3.00,N*KB-1
 4.40 G
10.05 S X=X(100*J+I)
10.10 I (X)10.4,10.2,10.4
10.20 T !!,"I CANNOT DEAL WITH MISSING VALUES"
10.25 T %2.00,!,"COL ",J,"  ROW ",I," IS MISSING"
10.30 S I=100;S J=100;S Q=0;R
10.40 S S=S+X;S Q=Q+X*X
11.10 D 10.05
11.20 S QR(I)=QR(I)+X
11.30 S QC(J)=QC(J)+X