C-PS/8 FOCAL, 1971
 1.84 C
 2.10 A !!!!,"ANALYSIS ",LI;IF (LI-3)2.2,4.1,2.2
 2.20 I (LI-4)2.3,5.3,2.3
 2.30 L R SORT 2.15
 3.50 C
 3.51 T !,"CORRELATION MATRIX"
 3.55 T !,:3;S U=0;F I=1,KB;T %7.00,I
 3.60 F BA=1,KB;T %2.00,!,BA;F BE=1,BA;DO 3.8
 3.61 F I=1,U;DO 3.63
 3.62 I (100-Z)12.10,12.10;T  !,%2.00," N= ",N;G 3.95
 3.63 I (N-Z(I))3.65,3.64,3.65
 3.64 S Z=0
 3.65 S Z=100;R
 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 I (BA-BE)3.92,3.85,3.92
 3.85 S R=1;T %6.03,R
 3.90 DO 3.7;F I=1,N;DO 9
 3.92 DO 3.90;DO 4.14;DO 4.16;S U=U+1;S Z(U)=N1;T %6.03,R
 3.95 G 2.1
 4.10 T !,"VARIABLES:   ";A "INDEPENDENT  "BA,"DEPENDENT  "BE,!
 4.12 D 3.7;D 3.90;T %8.04
 4.14 S SQ=SQ-S1*S2/N1;S Q1=Q1-S1^2/N1;S Q2=Q2-S2^2/N1
 4.16 S R=SQ/FSQT(Q1*Q2)
 4.18 T !,"Y=A+BX",!,"A=    ";S B=SQ/Q1;T S2/N1-B*S1/N1
 4.20 T "    B= ",B,!;S SQ=SQ^2/Q1;S Q2=Q2-SQ;S S=Q2/(N1-2)
 4.22 T "SD A=",FSQT(S*(1/N1+(S1/N1)^2/Q1))," SD B= ",FSQT(S/Q1)
 4.24 T "  TOTAL VARIANCE",S,!,"T VALUE FOR B",B/(FSQT(S/Q1))
 4.26 T !,"EXPL. SS",SQ,"  UNEXP. SS",Q2;G 2.1
 5.30 D 4.1;D 8.35;D 8.43;D 8.45;S N3=N1
 5.32 T %8.03
 5.33 S DF=BA;S T=BE
 5.38 S N(1)=FITR(N3/3+1); S N(2)=N3-2*N(1)
 5.40 D 3.7; F I=1,N(1);D 9
 5.42 S IN(1)=S1;S IN(2)=S2;D 3.7;F I=N(1)+N(2)+1,N3;D 9
 5.44 S M=(S1-IN(1))/N(1)
 5.45 S BA=DF;S BE=T
 5.46 S B=((S2-IN(2))/(S1-IN(1)));D 3.9;S A=S2/N1-B*S1/N3
 5.48 T !,"Y=A+BX",!
A  
,
  B 
;I (N1-5)2.1,2.1;
 5.50 A !, "CONFIDENCE LIMITS 
",DF;I (DF-0YES)2.1,5.52,2.1
 5.52 T "T FOR  ",%3.00,N3-3,"  DEGREES OF FREEDOM  ";A T,!
 5.57 F I=1,6;S IN(I)=0
 5.60 S A=1;S K=1;S Z=N(1);D 11
 5.65 S A=N(1)+1;S K=2;S Z=N(1)+N(2);D 11
 5.70 S A=Z+1;S K=1;S Z=N3; D 11
 5.75 S Q2=IN(2)-IN(5); S Q1=IN(1)-IN(4)
 5.80 S SQ=IN(3)-IN(6)
 5.85 S DF=N(1)*M^2*(N3-3)/(2*T^2)
 5.90 S I=((B*DF-SQ))/(DF-Q1)
 5.95 S J=(DF*(B^2*Q1+Q2-2*B*SQ)-Q1*Q2-SQ^2)
 5.96 I (J)5.98,5.98;S J=FSQT(J)/(DF-Q1)
 5.97 T %8.03, !, "UPPER LIMIT ",I+J," LOWER LIMIT ",I-J; G 5.5
 5.98 T ! "IMPOSSIBLE,TOO LARGE ERRORS";G 5.5
 6.10 S X=X(100*BA+I); S Y=X(100*BE+I)
 6.15 I (X-SA)6.3,6.9,6.3
 6.30 I (Y-SA)6.4,6.9,6.4
 6.40 S N1=N1+1
 6.50 S X(N1)=X;S X(N1-100)=Y
 6.90 R
 8.20 S S1=S1+X; S Q1=Q1+X^2; S N1=N1+1;R
 8.35 S N1=0;F I=1,N; D 6
 8.43 S VT=0;S VP=-1
 8.45 F I=1,N;D 10
 9.10 S X=X(BA*100+I);S Y=X(BE*100+I)
 9.20 I (X-SA)9.3,9.9,9.3
 9.30 I (Y-SA)9.4,9.9,9.4
 9.40 S S2=S2+Y;D 8.2;S Q2=Q2+Y^2;S SQ=SQ+X*Y;S DD=DD+(X-Y)^2
 9.50 S D=D+X-Y
 9.90 R
10.05 S J=N1
10.10 S X1=X(100*VT+I);S X2=X(100*VT+J)
10.20 I (X1-X2)10.4,10.3,10.4
10.30 S X(100*VT+I)=X2;S X(100*VT+J)=X1
10.35 S X(100*VP+I)=X(100*VP+I);S X(100*VP+J)=X(100*VP+J);D 10.6
10.40 S J=J-1; I (J-1)10.9,10.9,10.1
10.60 S X1=X(100*VP+I);S X(100*VP+I)=X(100*VP+J);S X(100*VP+J)=X1
10.90 R
11.10 D 3.7;F I=A,Z;D 9
11.15 S IN(1)=Q1+IN(1);S IN(2)=IN(2)+Q2
11.20 S IN(3)=SQ+IN(3);S IN(4)=IN(4)+S1^2/N(K)
11.25 S IN(5)=IN(5)+S2^2/N(K);S IN(6)=IN(6)+S1*S2/N(K)
12.10 T !!,"  NUMBER MATRIX"
12.15 DO 3.55
12.20 F BA=1,KB;T %2.00,!,BA;F BE=1,BA;DO 12.30
12.25 GO 2.10
12.30 I (BA-BE)12.40,12.35,12.40
12.35 T "      **"
12.40 S U=U+1;T :4, %7.00,Z(U)