File STAT8.FC (FOCAL source file)

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

C-PS/8 FOCAL, 1971
 1.10 A !!,"ANALYSIS ",LI;I (LI-36)1.2,1.3,1.2
 1.20 L R SORT 2.15
 1.30 S S=0;L R STAT6 2.4;WE CANNOT USE SAVED MEANS
 2.10 S S=0;S Q=0
 2.20 F I=1,KB;S QC(I)=0
 2.30 F I=1,N;S QR(I)=0
 3.05 T !!,"ANOVA",!
 3.10 F I=1,N;F J=1,KB;DO 10
 3.20 F I=1,N;F J=1,KB;DO 11
 3.30 I (Q) 1.1,1.1,3.4
 3.40 S S1=0;F I=1,KB;S S1=S1+QC(I)*QC(I)/N
 3.50 S S2=0;F I=1,N;S S2=S2+QR(I)*QR(I)/KB
 3.60 S SQ=S*S/(KB*N)
 3.70 T !,%8.04,,,"SOURCE   SS       DF      MS         F
 3.80 S S1=S1-SQ;S S2=S2-SQ
 3.90 S SQ=Q-SQ; S Q=SQ-S1-S2
 4.05 S DF=N*KB-(N+KB)+1
 4.07 T !,"COL ",%8.04,S1,%3.00,KB-1,%8.04,S1/(KB-1)
 4.10 T S1*DF/(Q*(KB-1)),!,"ROWS",S2,%3.00,N-1,%8.04,S2/(N-1)
 4.20 T S2*DF/(Q*(N-1)),!,"ERR "
 4.30 T Q,%3.00,DF,%8.04,Q/DF,!,"TOT ",SQ,%3.00,N*KB-1
 4.40 G
10.05 S X=X(100*J+I)
10.10 I (X)10.4,10.2,10.4
10.20 T !!,"I CANNOT DEAL WITH MISSING VALUES"
10.25 T %2.00,!,"COL ",J,"  ROW ",I," IS MISSING"
10.30 S I=100;S J=100;S Q=0;R
10.40 S S=S+X;S Q=Q+X*X
11.10 D 10.05
11.20 S QR(I)=QR(I)+X
11.30 S QC(J)=QC(J)+X



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