C-PS/8 FOCAL, 1971 1.01 C FROM STAT1 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 DETERMINES VARIOUS STATISTICS FROM A 1.07 C MAXIMUM OF TWO SAMPLE GROUPS. IT WILL ASK YOU FOR THE 1.08 C NUMBER OF GROUPS, AND THEN FOR THE NUMBER OF SAMPLES 1.09 C IN THE FIRST GROUP. WHEN THAT HAS BEEN ENTERED, IT 1.10 C WILL HAVE YOU INPUT THAT MANY SAMPLES, AND FIGURE 1.11 C CERTAIN STATISTICS FOR THAT GROUP. THEN IT WILL ASK FOR 1.12 C THE NUMBER OF SAMPLES IN THE OTHER GROUP IF YOU ASKED FOR 1.13 C TWO GROUPS, FIGURE ITS STATISTICS, AND THEN TWO-GROUP 1.14 C STATISTICS. 1.15 1.26 I (FITR(/5)-(/5)) 1.27,1.3 1.27 T "YOU DO NOT HAVE ENOUGH FEATURES. " 1.28 T "RESTART FOCAL WITH FEATURE 3."!;Q 1.29 E 1.30 S X=FOUT(1);F K=0,25;S X=FOUT(24) 1.31 T !!!%10.04;A "NUMBER OF GROUPS:"G;S L=0;I (G-2) 1.32,1.32,4.1 1.32 S S=0;T %1.00,"HOW MANY SAMPLES IN GROUP",L+1;A %10.04," ",N 1.33 F C=1,N;A "SAMPLE:"X(C);S S=S+X(C) 1.34 T !!"NO MORE INPUT!"!! 1.36 S K=N-1;S I=1;S J=I+1;D 7 1.37 S M=S/N;S SZ=0;F I=1,N;D 4 1.38 S SZ=SZ/(N-1);S SD=FSQT(SZ);S DF=(N-1) 1.39 S SE=FSQT(SZ/N);S ST=M/SE;S R=X(N)-X(1);T !!"MAX":30,X(N)!! 1.40 T "MIN":30,X(1)!!"RANGE":30,R,!!"MEAN":30,M,!!"VARIANCE":30,SZ,!! 1.41 T "STANDARD DEVIATION":30,SD,!!"STANDARD ERROR":30,SE,!! 1.42 T "T-STATISTIC":30,ST,!!"DEGREES OF FREEDOM":30,DF,!! 1.43 S L=L+1;I (SZ) 1.44,1.45,1.44 1.44 I (G-2) 8.01,2.01,8.01 1.45 T "VARIANCE OF 0"!!G 8.01 2.01 S MM(L)=M;S NN(L)=N;S SS(L)=SZ;S SX(L)=SE 2.02 I (L-2)1.32,2.03 2.03 S DM=MM(1)-MM(2) 2.04 S PC=((NN(1)-1)*SS(1)+(NN(2)-1)*SS(2))/(NN(1)+NN(2)-2) 2.05 S PD=(1/NN(1)+1/NN(2)) 2.06 S PE=FSQT(PC*PD);S TT=DM/PE;S DD=NN(1)+NN(2)-2 2.07 I (SS(1)-SS(2))2.08,2.09,2.09 2.08 S VR=SS(1)/SS(2);G 2.10 2.09 S VR=SS(2)/SS(1) 2.10 T !"DIFFERENCE BETWEEN MEANS":30,DM,!!"POOLED STANDARD ERROR":30,PE,!! 2.11 T "TWO GROUP T-STATISTIC":30,TT,!!"VARIANCE RATIO":30,VR,!! 2.12 T "DEGREES OF FREEDOM":30,MF,!!;I (SS(1)) 2.13,4.1,2.13 2.13 I (SS(2)) 3.01,8.01,3.01 3.01 S AC=FSQT(SS(1)/NN(1)+SS(2)/NN(2)) 3.02 S X1=(MM(1)-MM(2))/AC 3.03 S RS=(SS(1)/NN(1)+SS(2)/NN(2))^2 3.04 S L=(SS(1)/NN(1))^2/(NN(1)-1)+(SS(2)/NN(2))^2/(NN(2)-1) 3.05 T !!"WELCH STATISTICS:"!!! 3.06 S K=RS/L;T "T-STATISTIC":30,X1,!!"DEGREES OF FREEDOM":30,K,!! 3.07 T "UNPOOLED STANDARD ERROR":30,AC,!!;G 8.01 4.01 S SZ=SZ+(X(I)-M)^2 4.02 R 4.10 T "TOO MANY GROUPS";G 8.01 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 !!"AGAIN " ,AN;I (AN-0YES) 8.03,1.29,8.03 8.03 Q