File: STAT2.FC of Tape: Sources/Focal/fc3
(.FC,.FO Focal format converted to 'WRITE' listing)
C-PS/8 FOCAL, 1971 01.01 C FROM STAT2 BY M.J. MC KEOWN (DECUS: FOCAL8-16) 01.02 C MODIFIED BY JOHN HUISMAN 01.03 C OMSI SOFTWARE DEVELOPMENT GROUP 01.04 C LAST CHANGE: 9/9/71 01.05 01.06 C THIS PROGRAM REQUIRES THAT YOU KNOW THE MEANS 01.07 C AND VARIANCES OF THE GROUPS YOU WISH TO WORK ON. TO 01.08 C GET THIS INFORMATION, SIMPLY RUN 'STAT1'. 01.09 01.10 C THE PROGRAM USES THE FOLLOWING ABBREVIATIONS: 01.11 01.12 C ** SSB - SUM OF SQUARES BETWEEN GROUPS 01.13 C ** DFB - DEGREES OF FREEDOM BETWEEN GROUPS 01.14 C ** MSB - MEAN SQUARE DIFFERENCE BETWEEN GROUPS 01.15 C ** SSW - SUM OF SQUARES WITHIN GROUPS 01.16 C ** DFW - DEGREES OF FREEDOM WITHIN GROUPS 01.17 C ** MSW - MEAN SQUARE DIFFERENCE WITHIN GROUPS 01.18 C USING THE POOLED VARIANCE 01.19 C ** SCT - SO-CALLED TOTAL SUM OF SQUARES 01.20 C ** SCF - SO-CALLED TOTAL DEGREES OF FREEDOM 01.21 C ** SSG - GRAND MEAN SUM OF SQUARES 01.22 C ** ODF - OVERALL DEGREES OF FREEDOM 01.23 C ** SST - TOTAL SUM OF DEGREES 01.24 C ** DFT - TOTAL DEGREES OF FREEDOM 01.25 C ** FR - F-RATIO (F-STATISTIC) 01.26 E 01.27 I (FITR(#/5)-(#/5)) 1.28,1.31,1.28 01.28 T "YOU DO NOT HAVE ENOUGH FEATURES." 01.29 T "RESTART FOCAL WITH FEATURE 3."!!;Q 01.31 S X=FOUT(1);F K=1,25;S X=FOUT(24) 01.32 T %1.0 !!!;A "NUMBER OF GROUPS:"G;S S=0 01.33 F J=1,G;T "NUMBER IN GROUP",J,":";A N(J);S S=S+N(J) 01.35 S X=0 01.36 F J=1,G;T "MEAN OF GROUP",J,":";A M(J);S X=X+(N(J)*M(J))/S 01.39 S SG=S*(X)^2 02.01 S SB=0;S J=1 02.02 F J=1,G;S SB=SB+N(J)*(M(J)-X)^2 02.03 02.04 02.05 S SW=0 02.06 F J=1,G;T "VARIANCE OF GROUP",J,":";A V(J);S SW=SW+(N(J)-1)*V(J) 02.07 02.08 02.09 S CT=SB+SW;S ST=CT+SG;S XB=G-1;S MB=SB/XB 03.01 S XW=0;S J=1 03.02 S XW=XW+(N(J)-1) 03.03 I (G-J)3.05,3.05,3.04 03.04 S J=J+1;G 3.02 03.05 S MW=SW/XW;S XR=MB/MW;S BW=XB+XW;S TT=BW+1;T ,!!!!!! 04.01 T %6.04 "SSB":6,SB,!"DFB":6,XB,!"MSB":6,MB,! 04.02 T "SSW":6,SW,!"DFW":6,XW,!"MSW":6,MW,! 04.03 T "SCT":6,CT,!"SCF":6,BW,! 04.04 T "SSG":6,SG,!"ODF":6,1,!"SST":6,ST,! 04.05 T "DFT":6,TT,!;T "FR":6,XR,!!! 04.06 A "MORE?"AN;I (AN-0YES)4.07,1.31,4.07 04.07 T "GOODBYE"!!;Q 07.01 I (X(I)-X(J))7.03,7.03,7.02 07.02 S TX=X(I);S X(I)=X(J);S X(J)=TX 07.03 S J=J+1;I (J-N)7.01,7.01,7.04 07.04 S I=I+1;I (I-K)7.05,7.05,7.06 07.05 S J=I+1;G 7.01 07.06 R 08.01 A "M?"E;I (E-2)8.03,1.31,1.31 08.02 T "E2",!!;G 8.01 08.03 Q