File STAT6.FC (FOCAL source file)

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

()0Z\Y&C-PS/8 FOCAL, 1971?MB!DTI (LI-36)2.3,2.4,2.3?MX<D^L R SORT 2.15?MB=F(S SD=FSQT((Q1-S1^2/N1)/(N1-1))?MC[GFS N1=0;S S1=0;S S2=0;S SQ=0;S Q1=0;S DD=0;S D=0;S Q2=0?MC(GPD 3.7; F I=1,N; D 8?MC;H^T !;S N3=0;S Q3=0;S S3=0;S Q4=0?MD
H 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?MD9H"S S3=S3^2/N3;S Q4=Q4-S3;S Q3=Q3-S3;S S3=Q4-Q3?MEWH$T !!,"ANOVA",!,"SOURCE    SS        DF    MS         F?ME2H&T !,"AMONG ",%8.04,Q3,%3.0,KB-1,%8.04,Q3/(KB-1)?MFQH(T Q3*(N3-KB)/(S3*
KB-1)),!,"WITHIN",S3,%3.0,N3-KB,%8.04?MY H*T S3/(N3-KB),!,"TOTAL ",Q4,%3.0,N3-1;S S=S3/(N3-KB)?MGMH.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?MGWH0F BA=1,KB; D 5?MG$H2S X=S/(N2-KB); D 15?MG/H4S M=(N2-KB)*X-M?MHBH6S A=(1/(3*(KB-1)))*(A-1/(N2
KB))?MHRH8S S=KB-1; S N2=(KB+1)/A^2?MH^H:S A=N2/(1-A+2/N2)?MH?H<T !,%8.02, "F",N2*M/(S*(A-M)),"  DF",S,"/",N2;I (DF)4.9,4.9?MI]H>T !,%8.04, "SMAX/SSUM",VT/VP,%3.00,"   V",N1,"   K",KB?MIIZG 2.1?MI1JJD 3.8;D 3.4;S SD=SD^2?MI?JLS VP=VP+SD;S N2=N2+N1?MJOJNS 
=SD;D 15;S M=M+(N1-1)*X?MJ!JPS S=S+(N1-1)*SD;S A=A+1/(N1-1)?MJ/JRI (VT-SD)5.2,5.2,5.22?MJ6JTS VT=SD?MKDJVI (N3) 5.28,5.28,5.24?MKSJXI (N3-N1)5.26,5.28,5.26?MKYJZS DF=0?MK!J\S N3=N1;R?MK6PES X=X(BA*100+I);I (X-SA)8.2,8.15,8.2?MK:POR?MLNPTS S1=S1+X;S Q1=Q1+
^2; S N1=N1+1;R?MLTP^D 4.32?ML!PS N1=0;F I=1,N; D 6?MMFP(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)?MMPP+S VT=0;S VP=-1?MMZP-F I=1,N1;D 10?MM"P/D 15;D 14?MM/P2S Y(1)=0; S Y(-1)=0;?MNFP4F I=N0,N1;S Y(X(I-100))=Y(X(I-100))+X(I)?MN
P7S N1=N1-N0+1;I (Y(1)-Y(-1)) 8.6;S Y(1)=Y(-1)?MN6P<T %6.02;S T=Y(1);I (N1-25) 8.65,8.65,8.7?MOCQAT !," T ",Y(1);G 2.1?MOQFT !,  ," Z ",(T-N1*(N1+1)/4)/FSQT((N1*(N1+1)*(2*N1+1))/24)?MO)QPG 2.1?MPF^JI (X^2-2.04*X+1)15.2; S X=FSQT(X); DO 15; S X=2*X;R?MZE^T
 X=(X-1)/(X+1); S X=2*(X+X^3/3+X^5/5+X^7/7)?MP([E?MP:(^A "NUMBER OF COLUMNS ",KB;F J=1,KB;D 28?MQF(A !, "ERROR M.S."S?MQQ((F I=1,J;S C(I)=0?MQ((-S Y=0;S N=1;A!, "NUMBER OF F VALUES ",NF?MRB(2T !,"STATE",NF,"  F-VALUES,SMALLEST FIRST  ";?MRI(7S I=NF+1?MR^(
S I=I-1;I (I-1)21.1; A XF(I);G 20.6?MR$)AG 21.1?MR)*JT !!?MR>*TT "GR. ";F Z=2,1,KB; T %3.00 " ",Z?MSC*^T !!?MSM.JS J=1;S P=J+1?MSX.(T %3,J,%5;S Q=1?MS.2I (J-Q)23.7,23.7?MS3.<T "     .";S Q=Q+1;G 23.5?MS9/FS R=1?MTE/PS C(J)=1;S C(P)=-1?MTR0JS H=0;;S V=0;S 
V=0?MTX0TS Q=1?MT10^S H=H+C(Q)*X(Q); S V=V+C(Q)^2*(S/X(-100+Q))?MT>0(I (Q-KB)24.5;G 24.6?MUH02S Q=Q+1;G 24.3?MU
0<S SV=(KB-1)*XF(R);S XX=FSQT(SV)*FSQT(V)?MU+1FS U=H+XX;S L=H-XX?MU92JI (L*U)25.2,25.2,25.4?MVR2TS R=R+1;I (R-NF)23.8,23.8;;T "     N";G 25.5?MVW2(D 27?MV22I (KB-P)25.7,25.7?MV62<S C(P)=0;S C(J)=0;S P=P+1;G 23.7?MWD3FI (KB-(J-1))26.1,26.1?MW!3PS C(P)=0;S 
(J)=0;S J=J+1;T !!;S P=J+1;I (J-KB)23.4?MY?4JT !!;E?MW46JF IS=1,5-NF+R;T " "?MXA6TF IS=5-NF+R,5;T "*"?M8JA !,"NUMBER ",X(-100+J),,"MEAN ",X(J)?MP
(TI (S)20.27,20.27;A "USE SAVED MEANS,"LI;I (LI-0YES)20.27,20.45,20.27?MB+D(I (N)20.2,20.2;A "  REALLY ",LI;I (LI-0YES)2.1,20.2,2.1?MF.H-2.1?MY[H,G 2.1?MBTDJA !!!!,"ANALYSIS ",LI;IF (LI-7)2.2,4.46,2.2?MW	4TG 20.1?MXW(JA !!,"ANALYSIS",LI;I (LI-36)2.3,20.2,
.3?M15L S STAT6?M?ALYSIS",LI;I (LI-36)2.3,20.2,2.3?M152,2.1?,20.45,20.27?MX6H68 H7%JFT+QH*%UX)2IL&A=H)\&2AP&P*P6+F-7NQSX%\XTL?8 X6L?%JFV+QB\0Y+"%X0^A;<(+"%NTZY/!H2/!A;X6:A%JFT2IL&A=H6\&2AP&//BCL,	$

Note: Lines longer than 256 characters were wrapped



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