C-FOCAL , 8/68

01.01 W 9
01.02 W 1
01.03 Q
01.05 A "A"A
01.06 S A=A+1;A V;S X=FNEW(A,V);G 1.06

02.01 E
02.03 A ?ORDER N?,!;S R=OR+1;F M=1,R;A A(M+10)
02.08 F M=1,N;D 15
02.10 F J=1,OR;F M=J,R;S A(10*J+M)=X(10*J+M)-X(M)*X(J)/N
02.20 F Z=2,OR;F J=1,Z-1;D 5
02.30 F B=1,OR;S B(B)=A(10*B+R)/A(11*B)
02.40 I (OR-2)2.5;F B=1,OR-1;F J=1,B;D 7
02.50 S B=X(R)/N;F J=1,OR;S B=B-B(J)*X(J)/N
02.70 S X=X(R)^2/N;F J=1,OR;S X(J)=A(10*J+R)^2/A(11*J)
02.84 S Z=X(11*R);T %8.08,!!,Z;F M=0,OR;S Z=Z-X(M)
02.86 T Z;F M=0,OR;T !,"B"M,B(M),"   VAR"X(M),(N-OR)*X(M)/Z
02.88 T !!,Z/(X(11*R)-X)
02.90 Q

05.10 F M=Z,R;S A(10*Z+M)=A(10*Z+M)-A(10*J+M)*A(10*J+Z)/A(11*J)

07.10 S B(OR-B)=B(OR-B)-B(R-J)*A(10*(OR-B)+R-J)/A(11*(OR-B))

09.30 COMMENTS-- N TH ORDER REGRESSION TO 9TH WITH OUT LOGS.
09.31 TO LOAD DATA: G 1.05 ;ENTER AN ADDRESS,THEN
09.32 A COLUMN OF DATA. STOP WITH CTRL C.
09.33 
09.34 PROGRAM ASKS FOR DATA LOCATIONS. THE DATA MUST
09.36 START AT LOCATION+1. THE DEPENDENT VARIABLE MUST BE LAST.
09.38 
09.40 TYPE OUT FORMAT IS:
09.42 TOTAL VARIANCE      RESIDUAL VARIANCE
09.44 B NO.    COEFF.       VARIANCE     "F" RATIO
09.45 FRACTION OF VARIANCE NOT EXPLAINED.
09.46 
09.48 
09.50 AFTER DATA IS LOADED, ERASE 1 & 9 AND GO TO RUN REGRESSION.
09.60 
09.62 TO SOLVE RESULTING EQU. FOR ONE LINE OF DATA: DO 11
09.64 11 MUST BE ERASED FOR 9TH ORDER.

11.02 D 2.03
11.10 S J=B;F I=1,OR;S J=J+B(I)*FNEW(A(I+10)+N)
11.20 T !,J,!

15.10 F I=1,R;S B(I)=FNEW(A(I+10)+M)
15.12 F I=1,R;S X(I)=X(I)+B(I)
15.20 F J=1,OR;F I=J,R;S X(J*10+I)=X(J*10+I)+B(J)*B(I)
15.40 S X(11*R)=X(11*R)+B(R)^2
*8~