CS()0WYTUC-PS/8 FOCAL, 1971?MB
DCS DM=MM(1)-MM(2)?MB=DDS PC=((NN(1)-1)*SS(1)+(NN(2)-1)*SS(2))/(NN(1)+NN(2)-2)?MCKDES PD=(1/NN(1)+1/NN(2))?MC%DFS PE=FSQT(PC*PD);S TT=DM/PE;S DD=NN(1)+NN(2)-2?MC7DGI (SS(1)-SS(2))2.08,2.09,2.09?MDFDHS VR=SS(1)/SS(2);G 2.10?MU\DIS VR=SS(2)/SS(1)?MD]HAS SZ=
Z+(X(I)-M)^2?MHCHBR?MD2NAI (X(I)-X(J))7.03,7.03,7.02?MEENBS TX=X(I);S X(I)=X(J);S X(J)=TX?MEWNCS J=J+1;I (J-N)7.01,7.01,7.04?ME)NDS I=I+1;I (I-K)7.05,7.05,7.06?ME3NES J=I+1;G 7.01?MS<NFR?MPCQ?MFKB$S K=N-1;S I=1;S J=I+1;D 7?MF[B%S M=S/N;S SZ=0;F I=1,N;D 4?
IDB&S SZ=SZ/(N-1);S SD=FSQT(SZ);S DF=(N-1)?MGQB,I (G-2) 8.01,2.01,8.01?MF1B+S L=L+1;I (SZ) 1.44,1.45,1.44?MKFB-T "VARIANCE OF 0"!!G 8.01?MG$BE?MM]BFC     THIS PROGRAM DETERMINES VARIOUS STATISTICS FROM A?MD!HJT "TOO MANY GROUPS";G 8.01?MQ-DMI (SS(2)) 3.01
8.01,3.01?ME;B"T !!"NO MORE INPUT!"!!?ML;B^S X=FOUT(1);F K=0,25;S X=FOUT(24)?MI%B	S SE=FSQT(SZ/N);S ST=M/SE;S R=X(N)-X(1);T !!"MAX":30,X(N)!!?MJIB(T "MIN":30,X(1)!!"RANGE":30,R,!!"MEAN":30,M,!!"VARIANCE":30,SZ,!!?MJ(B)T "STANDARD DEVIATION":30,SD,!!"STAND
RD ERROR":30,SE,!!?MF?B*T "T-STATISTIC":30,ST,!!"DEGREES OF FREEDOM":30,DF,!!?MUQDAS MM(L)=M;S NN(L)=N;S SS(L)=SZ;S SX(L)=SE?MT2DKT "TWO GROUP T-STATISTIC":30,TT,!!"VARIANCE RATIO":30,VR,!!?MLQBBC     MODIFIED BY JOHN HUISMAN?MWIBCC     OMSI SOFTWARE DEVE
OPMENT GROUP?MH"B!F C=1,N;A "SAMPLE:"X(C);S S=S+X(C)?MV	B
T !!!%10.04;A "NUMBER OF GROUPS:"G;S L=0;I (G-2) 1.32,1.32,4.1?MM<BGC MAXIMUM OF TWO SAMPLE GROUPS.  IT WILL ASK YOU FOR THE?MNZBHC NUMBER OF GROUPS, AND THEN FOR THE NUMBER OF SAMPLES?MN8BIC IN THE FIRST GROUP.  WHEN THAT HAS BEEN ENTERED, IT?MOUBJC WILL
HAVE YOU INPUT THAT MANY SAMPLES, AND FIGURE?MO5BKC CERTAIN STATISTICS FOR THAT GROUP.  THEN IT WILL ASK FOR?MPVBLC THE NUMBER OF SAMPLES IN THE OTHER GROUP IF YOU ASKED FOR?MQBMC TWO GROUPS, FIGURE ITS STATISTICS, AND THEN TWO-GROUP?MP8BO?MQIBZI (FITR(/5
-(/5)) 1.27,1.3?MVRB[T "YOU DO NOT HAVE ENOUGH FEATURES.  "?MH0B]E?MP5BNC STATISTICS.?MRAFAS AC=FSQT(SS(1)/NN(1)+SS(2)/NN(2))?MROFBS X1=(MM(1)-MM(2))/AC?MR"FCS RS=(SS(1)/NN(1)+SS(2)/NN(2))^2?MVBFDS L=(SS(1)/NN(1))^2/(NN(1)-1)+(SS(2)/NN(2))^2/(NN(2)-1)?MSF
S K=RS/L;T "T-STATISTIC":30,X1,!!"DEGREES OF FREEDOM":30,K,!!?MDQFGT "UNPOOLED STANDARD ERROR":30,AC,!!;G 8.01?ME7PAA !!"AGAIN
" ,AN;I (AN-0YES) 8.03,1.29,8.03?MK?BAC     FROM STAT1 BY M.J. MC KEOWN (DECUS: FOCAL8-16)?MHSDLT "DEGREES OF FREEDOM":30,MF,!!;I (SS(1)) 2.13,4.1,2.13?MBTDBI (L-2)1.32,2.03?MK^DJT !"DIFFERENCE BETWEEN MEANS":30,DM,!!"POOLED STANDARD ERROR":30,PE,!!?MSAFE
 !!"WELCH STATISTICS:"!!!?MQ
B\T "RESTART FOCAL WITH FEATURE 3."!;Q?ML	B S S=0;T %1.00,"HOW MANY SAMPLES IN GROUP",L+1;A %10.04,"
 ",N?MG!BDC     LAST CHANGE: 8/19/71?MT:L S STAT1?M? CHANGE: 8/19/71?MT:,E;L R STAT1?M?M%10.04,"
 ",N?TING LOCATION:3M;A3D^ "M&J3DGC2 
CO7Q2JNI3C/H*I>I3SWE2 [I* 3MN7G2J*I+MCP3C?G+B>IN7O6JBI3USJ*M2INI3HSR"+&22/WD7LSE3E6 ",7N2J
C+O3A2IGT3D& 2 
CO	Q6/WN3MWB3R> 3FN 7HKA3S& 7NF 2-WRG*M6IWP3S+H"M&J6ICUP	N2/?M+VCE7TE3TJ 3U[F+ECRU7P6JFI3S7U*I>I3C7O7PSU7EF 2-WR7GB 3O	S

Note: Lines longer than 256 characters were wrapped