File STAT6.FC (FOCAL source file)

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

C-PS/8 FOCAL, 1971
 2.10 A !!!!,"ANALYSIS ",LI;IF (LI-7)2.2,4.46,2.2
 2.20 I (LI-36)2.3,2.4,2.3
 2.30 L R SORT 2.15
 2.40 I (N)20.2,20.2;A "  REALLY ",LI;I (LI-0YES)2.1,20.2,2.1
 3.40 S SD=FSQT((Q1-S1^2/N1)/(N1-1))
 3.70 S N1=0;S S1=0;S S2=0;S SQ=0;S Q1=0;S DD=0;S D=0;S Q2=0
 3.80 D 3.7; F I=1,N; D 8
 4.30 T !;S N3=0;S Q3=0;S S3=0;S Q4=0
 4.32 F BA=1,KB;D 3.8;S N3=N3+N1;S Q3=Q3+(S1)^2/N1;S S3=S3+S1;S Q4=Q4+Q1
 4.34 S S3=S3^2/N3;S Q4=Q4-S3;S Q3=Q3-S3;S S3=Q4-Q3
 4.36 T !!,"ANOVA",!,"SOURCE    SS        DF    MS         F
 4.38 T !,"AMONG ",%8.04,Q3,%3.0,KB-1,%8.04,Q3/(KB-1)
 4.40 T Q3*(N3-KB)/(S3*(KB-1)),!,"WITHIN",S3,%3.0,N3-KB,%8.04
 4.42 T S3/(N3-KB),!,"TOTAL ",Q4,%3.0,N3-1;S S=S3/(N3-KB)
 4.44 G 2.1
 4.45 2.1
 4.46 T !;S VT=0;S VP=0;S A=0;S M=0;S S=0;S N3=0;S DF=1;S N2=0
 4.48 F BA=1,KB; D 5
 4.50 S X=S/(N2-KB); D 15
 4.52 S M=(N2-KB)*X-M
 4.54 S A=(1/(3*(KB-1)))*(A-1/(N2-KB))
 4.56 S S=KB-1; S N2=(KB+1)/A^2
 4.58 S A=N2/(1-A+2/N2)
 4.60 T !,%8.02, "F",N2*M/(S*(A-M)),"  DF",S,"/",N2;I (DF)4.9,4.9
 4.62 T !,%8.04, "SMAX/SSUM",VT/VP,%3.00,"   V",N1,"   K",KB
 4.90 G 2.1
 5.10 D 3.8;D 3.4;S SD=SD^2
 5.12 S VP=VP+SD;S N2=N2+N1
 5.14 S X=SD;D 15;S M=M+(N1-1)*X
 5.16 S S=S+(N1-1)*SD;S A=A+1/(N1-1)
 5.18 I (VT-SD)5.2,5.2,5.22
 5.20 S VT=SD
 5.22 I (N3) 5.28,5.28,5.24
 5.24 I (N3-N1)5.26,5.28,5.26
 5.26 S DF=0
 5.28 S N3=N1;R
 8.05 S X=X(BA*100+I);I (X-SA)8.2,8.15,8.2
 8.15 R
 8.20 S S1=S1+X;S Q1=Q1+X^2; S N1=N1+1;R
 8.30 D 4.32
 8.35 S N1=0;F I=1,N; D 6
 8.40 F I=1,N1;S X=X(I);S Y=X(I-100);S X(I)=FABS(X-Y);S X(I-100)=FSGN(X-Y)
 8.43 S VT=0;S VP=-1
 8.45 F I=1,N1;D 10
 8.47 D 15;D 14
 8.50 S Y(1)=0; S Y(-1)=0;
 8.52 F I=N0,N1;S Y(X(I-100))=Y(X(I-100))+X(I)
 8.55 S N1=N1-N0+1;I (Y(1)-Y(-1)) 8.6;S Y(1)=Y(-1)
 8.60 T %6.02;S T=Y(1);I (N1-25) 8.65,8.65,8.7
 8.65 T !," T ",Y(1);G 2.1
 8.70 T !,  ," Z ",(T-N1*(N1+1)/4)/FSQT((N1*(N1+1)*(2*N1+1))/24)
 8.80 G 2.1
15.10 I (X^2-2.04*X+1)15.2; S X=FSQT(X); DO 15; S X=2*X;R
15.20 S X=(X-1)/(X+1); S X=2*(X+X^3/3+X^5/5+X^7/7)
20.10 A !!,"ANALYSIS",LI;I (LI-36)2.3,20.2,2.3
20.20 I (S)20.27,20.27;A "USE SAVED MEANS,"LI;I (LI-0YES)20.27,20.45,20.27
20.27 E
20.30 A "NUMBER OF COLUMNS ",KB;F J=1,KB;D 28
20.35 A !, "ERROR M.S."S
20.40 F I=1,J;S C(I)=0
20.45 S Y=0;S N=1;A!, "NUMBER OF F VALUES ",NF
20.50 T !,"STATE",NF,"  F-VALUES,SMALLEST FIRST  ";
20.55 S I=NF+1
20.60 S I=I-1;I (I-1)21.1; A XF(I);G 20.6
20.65 G 21.1
21.10 T !!
21.20 T "GR. ";F Z=2,1,KB; T %3.00 " ",Z
21.30 T !!
23.10 S J=1;S P=J+1
23.40 T %3,J,%5;S Q=1
23.50 I (J-Q)23.7,23.7
23.60 T "     .";S Q=Q+1;G 23.5
23.70 S R=1
23.80 S C(J)=1;S C(P)=-1
24.10 S H=0;;S V=0;S SV=0
24.20 S Q=1
24.30 S H=H+C(Q)*X(Q); S V=V+C(Q)^2*(S/X(-100+Q))
24.40 I (Q-KB)24.5;G 24.6
24.50 S Q=Q+1;G 24.3
24.60 S SV=(KB-1)*XF(R);S XX=FSQT(SV)*FSQT(V)
24.70 S U=H+XX;S L=H-XX
25.10 I (L*U)25.2,25.2,25.4
25.20 S R=R+1;I (R-NF)23.8,23.8;;T "     N";G 25.5
25.40 D 27
25.50 I (KB-P)25.7,25.7
25.60 S C(P)=0;S C(J)=0;S P=P+1;G 23.7
25.70 I (KB-(J-1))26.1,26.1
25.80 S C(P)=0;S C(J)=0;S J=J+1;T !!;S P=J+1;I (J-KB)23.4
26.10 T !!;E
26.20 G 20.1
27.10 F IS=1,5-NF+R;T " "
27.20 F IS=5-NF+R,5;T "*"
28.10 A !,"NUMBER ",X(-100+J),,"MEAN ",X(J)



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