File: STAT1.FC of Tape: Sources/Focal/fc3
(.FC,.FO Focal format converted to 'WRITE' listing)
C-PS/8 FOCAL, 1971 01.01 C FROM STAT1 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 DETERMINES VARIOUS STATISTICS FROM A 01.07 C MAXIMUM OF TWO SAMPLE GROUPS. IT WILL ASK YOU FOR THE 01.08 C NUMBER OF GROUPS, AND THEN FOR THE NUMBER OF SAMPLES 01.09 C IN THE FIRST GROUP. WHEN THAT HAS BEEN ENTERED, IT 01.10 C WILL HAVE YOU INPUT THAT MANY SAMPLES, AND FIGURE 01.11 C CERTAIN STATISTICS FOR THAT GROUP. THEN IT WILL ASK FOR 01.12 C THE NUMBER OF SAMPLES IN THE OTHER GROUP IF YOU ASKED FOR 01.13 C TWO GROUPS, FIGURE ITS STATISTICS, AND THEN TWO-GROUP 01.14 C STATISTICS. 01.15 E 01.26 I (FITR(#/5)-(#/5)) 1.27,1.3 01.27 T "YOU DO NOT HAVE ENOUGH FEATURES. " 01.28 T "RESTART FOCAL WITH FEATURE 3."!;Q 01.29 01.30 S X=FOUT(1);F K=0,25;S X=FOUT(24) 01.31 T !!!%10.04;A "NUMBER OF GROUPS:"G;S L=0;I (G-2) 1.32,1.32,4.1 01.32 S S=0;T %1.00,"HOW MANY SAMPLES IN GROUP",L+1;A %10.04,"? ",N 01.33 F C=1,N;A "SAMPLE:"X(C);S S=S+X(C) 01.34 T !!"NO MORE INPUT!"!! 01.36 S K=N-1;S I=1;S J=I+1;D 7 01.37 S M=S/N;S SZ=0;F I=1,N;D 4 01.38 S SZ=SZ/(N-1);S SD=FSQT(SZ);S DF=(N-1) 01.39 S SE=FSQT(SZ/N);S ST=M/SE;S R=X(N)-X(1);T !!"MAX":30,X(N)!! 01.40 T "MIN":30,X(1)!!"RANGE":30,R,!!"MEAN":30,M,!!"VARIANCE":30,SZ,!! 01.41 T "STANDARD DEVIATION":30,SD,!!"STANDARD ERROR":30,SE,!! 01.42 T "T-STATISTIC":30,ST,!!"DEGREES OF FREEDOM":30,DF,!! 01.43 S L=L+1;I (SZ) 1.44,1.45,1.44 01.44 I (G-2) 8.01,2.01,8.01 01.45 T "VARIANCE OF 0"!!G 8.01 02.01 S MM(L)=M;S NN(L)=N;S SS(L)=SZ;S SX(L)=SE 02.02 I (L-2)1.32,2.03 02.03 S DM=MM(1)-MM(2) 02.04 S PC=((NN(1)-1)*SS(1)+(NN(2)-1)*SS(2))/(NN(1)+NN(2)-2) 02.05 S PD=(1/NN(1)+1/NN(2)) 02.06 S PE=FSQT(PC*PD);S TT=DM/PE;S DD=NN(1)+NN(2)-2 02.07 I (SS(1)-SS(2))2.08,2.09,2.09 02.08 S VR=SS(1)/SS(2);G 2.10 02.09 S VR=SS(2)/SS(1) 02.10 T !"DIFFERENCE BETWEEN MEANS":30,DM,!!"POOLED STANDARD ERROR":30,PE,!! 02.11 T "TWO GROUP T-STATISTIC":30,TT,!!"VARIANCE RATIO":30,VR,!! 02.12 T "DEGREES OF FREEDOM":30,MF,!!;I (SS(1)) 2.13,4.1,2.13 02.13 I (SS(2)) 3.01,8.01,3.01 03.01 S AC=FSQT(SS(1)/NN(1)+SS(2)/NN(2)) 03.02 S X1=(MM(1)-MM(2))/AC 03.03 S RS=(SS(1)/NN(1)+SS(2)/NN(2))^2 03.04 S L=(SS(1)/NN(1))^2/(NN(1)-1)+(SS(2)/NN(2))^2/(NN(2)-1) 03.05 T !!"WELCH STATISTICS:"!!! 03.06 S K=RS/L;T "T-STATISTIC":30,X1,!!"DEGREES OF FREEDOM":30,K,!! 03.07 T "UNPOOLED STANDARD ERROR":30,AC,!!;G 8.01 04.01 S SZ=SZ+(X(I)-M)^2 04.02 R 04.10 T "TOO MANY GROUPS";G 8.01 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 !!"AGAIN?" ,AN;I (AN-0YES) 8.03,1.29,8.03 08.03 Q