C-FOCAL,1969

03.01 TYPE "ANALYSIS OF VARIANCE"
03.10 ASK,!!, "NO. OF GROUPS",G
03.20 SET SN=0
03.30 FOR I=1,G;DO 4
03.40 GOTO 5.1

04.10 TYPE,!, "GROUP",I," ,SIZE"
04.12 ASK N(I)
04.20 SET SN=SN+N(I)
04.21 TYPE "GROUP DATA",I
04.30 FOR C=1,N(I);ASK SAM(C)
04.40 SET SUM(I)=0
04.50 FOR C=1,N(I);SET SUM(I)=SUM(I)+SAM(C)
04.60 SET IDS(I)=0
04.70 FOR C=1,N(I);SET IDS(I)=SAM(C)^2+IDS(I)
04.80 RETURN

05.10 SET DCF=0
05.20 FOR I=1,G;SET DCF=DCF+SUM(I)
05.30 SET CF=DCF^2/SN
05.40 SET SSX=0
05.50 FOR I=1,G;SET SSX=SSX+IDS(I)
05.60 SET TSS=SSX-CF
05.70 SET PBS=0
05.80 FOR I=1,G;SET PBS=PBS+SUM(I)^2/N(I)
05.90 SET BSS=PBS-CF

06.10 SET WSS=TSS-BSS
06.20 SET DB=G-1
06.30 SET MB=BSS/DB
06.40 SET DW=SN-G
06.50 SET MW=WSS/DW
06.60 TYPE %4.03,!,"F SCORE", MB/MW
06.70 T,!!, "DEGREES OF FREEDOM BETWEEN GROUPS", DB
06.71 T,!, "DEGREES OF FREEDOM WITHIN GROUPS ", DW
06.72 T,!,"WITHIN GROUP MEANS SQUARED       ",MW
06.73 QUIT
*