File STAT1.FC (FOCAL source file)

Directory of image this file is from
This file as a plain text file

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



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search