C-PS/8 FOCAL, 1971
 1.01 C     FROM STAT2 BY M.J. MC KEOWN (DECUS: FOCAL8-16)
 1.02 C     MODIFIED BY JOHN HUISMAN
 1.03 C     OMSI SOFTWARE DEVELOPMENT GROUP
 1.04 C     LAST CHANGE: 8/19/71
 1.05 
 1.06 C     THIS PROGRAM REQUIRES THAT YOU KNOW THE MEANS 
 1.07 C AND VARIANCES OF THE GROUPS YOU WISH TO WORK ON.  TO 
 1.08 C GET THIS INFORMATION, SIMPLY RUN 	STAT1	.
 1.09 
 1.10 C     THE PROGRAM USES THE FOLLOWING ABBREVIATIONS:
 1.11 
 1.12 C ** SSB - SUM OF SQUARES BETWEEN GROUPS
 1.13 C ** DFB - DEGREES OF FREEDOM BETWEEN GROUPS
 1.14 C ** MSB - MEAN SQUARE DIFFERENCE BETWEEN GROUPS
 1.15 C ** SSW - SUM OF SQUARES WITHIN GROUPS
 1.16 C ** DFW - DEGREES OF FREEDOM WITHIN GROUPS
 1.17 C ** MSW - MEAN SQUARE DIFFERENCE WITHIN GROUPS
 1.18 C          USING THE POOLED VARIANCE
 1.19 C ** SCT - SO-CALLED TOTAL SUM OF SQUARES
 1.20 C ** SCF - SO-CALLED TOTAL DEGREES OF FREEDOM
 1.21 C ** SSG - GRAND MEAN SUM OF SQUARES
 1.22 C ** ODF - OVERALL DEGREES OF FREEDOM
 1.23 C ** SST - TOTAL SUM OF DEGREES
 1.24 C ** DFT - TOTAL DEGREES OF FREEDOM
 1.25 C ** FR  - F-RATIO (F-STATISTIC)
 1.26 
 1.27 I (FITR(/5)-(/5)) 1.28,1.31,1.28
 1.28 T "YOU DO NOT HAVE ENOUGH FEATURES."
 1.29 T "RESTART FOCAL WITH FEATURE 3."!!;Q
 1.31 S X=FOUT(1);F K=1,20;S X=FOUT(24)
 1.32 T %1.0 !!!;A "NUMBER OF GROUPS:"G;S S=0
 1.33 F J=1,G;T "NUMBER IN GROUP",J,":";A N(J);S S=S+N(J)
 1.35 S X=0
 1.36 F J=1,G;T "MEAN OF GROUP",J,":";A M(J);S X=X+(N(J)*M(J))/S
 1.39 S SG=S*(X)^2
 2.01 S SB=0;S J=1
 2.02 F J=1,G;S SB=SB+N(J)*(M(J)-X)^2
 2.03 
 2.04 
 2.05 S SW=0
 2.06 F J=1,G;T "VARIANCE OF GROUP",J,":";A V(J);S SW=SW+(N(J)-1)*V(J)
 2.07 
 2.08 
 2.09 S CT=SB+SW;S ST=CT+SG;S XB=G-1;S MB=SB/XB
 3.01 S XW=0;S J=1
 3.02 S XW=XW+(N(J)-1)
 3.03 I (G-J)3.05,3.05,3.04
 3.04 S J=J+1;G 3.02
 3.05 S MW=SW/XW;S XR=MB/MW;S BW=XB+XW;S TT=BW+1;T ,!!!!!!
 4.01 T %6.04 "SSB":6,SB,!"DFB":6,XB,!"MSB":6,MB,!
 4.02 T "SSW":6,SW,!"DFW":6,XW,!"MSW":6,MW,!
 4.03 T "SCT":6,CT,!"SCF":6,BW,!
 4.04 T "SSG":6,SG,!"ODF":6,1,!"SST":6,ST,!
 4.05 T "DFT":6,TT,!;T "FR":6,XR,!!!
 4.06 A "MORE
"AN;I (AN-0YES)4.07,1.31,4.07
 4.07 T "GOODBYE"!!;Q
 7.01 I (X(I)-X(J))7.03,7.03,7.02
 7.02 S TX=X(I);S X(I)=X(J);S X(J)=TX
 7.03 S J=J+1;I (J-N)7.01,7.01,7.04
 7.04 S I=I+1;I (I-K)7.05,7.05,7.06
 7.05 S J=I+1;G 7.01
 7.06 R
 8.01 A "M
"E;I (E-2)8.03,1.31,1.31
 8.02 T "E2",!!;G 8.01
 8.03 Q