@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C-8K FOCAL @1969

01.01 C***** HYPERB *****   OCT 13 1971
01.02 C-----A FOCAL ADAPTATION OF THE YALE COMPUTER CENTER FORTRAN
01.03 C-----   PROGRAM 'HYPERB' USED FOR FITTING DATA TO THE
01.04 C-----      MICHAELIS-MENTEN EQUATION
01.05 C-----         ...STANLEY R. VIVIAN...
01.10 E
01.20 T !!"          MICHAELIS-MENTEN REGRESSION
01.25 A !!"N:"N
01.26 I (N-3)1.28;I (N-31)1.3
01.27 T !!"TOO MANY DATA POINTS - MAXIMUM N IS 30";G 5.5
01.28 T !!"TOO FEW DATA POINTS - MINIMUM N IS 3";G 5.5
01.30 T !"   S     V";F I=1,N;A !":"XP(I)," :"Y(I)
01.40 S D=DO;S IT=0
01.50 S X=0;S XS=0;S Y=0;S X2=0;S Y2=0;S S2=0;S XX=0;S XY=0;S SY=0
01.60 F I=1,N;S XI=XP(I)/(XP(I)+D);S XJ=XI/(XP(I)+D);D 1.7;D 1.9;D 1.92
01.65 G 2.1
01.70 S X=X+XI;S XS=XS+XJ;S Y=Y+Y(I)
01.90 S Y2=Y2+Y(I)^2;S S2=S2+XJ^2;S XX=XX+XI*XJ;S X2=X2+XI^2
01.92 S XY=XY+XI*Y(I);S SY=SY+XJ*Y(I)

02.10 S B1=XY*S2-SY*XX;S B2=SY*X2-XY*XX;I (B1)2.12,2.4,2.12
02.12 S BR=B2/B1
02.15 I (SW-1)2.3,2.2,2.3
02.20 T !%2,IT,%5.1,D,X,XS,X2,XX,S2,XY,SY," "%8.06,BR
02.25 I (IT-15)2.3;G 2.4
02.30 I (FABS(BR)-.000001)2.6,2.6;I (IT-40)2.5,2.5;T "";C--- 10 BELLS
02.40 I (SW-1)2.42;S SW=SW+1;G 2.47
02.42 T !!"FAILURE AT CONVERGENCE - K AND VMAX ARE NEGATIVE
02.44 A !"DO YOU WANT AN ITERATION TABLE? "IT
02.46 I (IT-0N)2.49,2.47;I (IT-0Y)2.49,2.48,2.49
02.47 A !"ENTER A NEGATIVE K OR TYPE '/' TO TERMINATE: "DO;I (DO)1.4;D 5
02.48 S SW=SW+1;D 2.52;D 2.54;G 1.4
02.49 T !"PLEASE ANSWER 'Y' OR 'N'";G 2.44
02.50 S D=D-BR;S IT=IT+1;G 1.5
02.52 T !!"          ITERATION TABLE",!"       K    
02.54 T "SUMX   SUMX*  SUMXX SUMXX* SUMX*X* SUMXY SUMX*Y   B-RATIO
02.60 S DS=X2*S2-XX^2;S B1=B1/DS;S B2=B2/DS;S BC=XY/X2;S NF=N-2
02.70 S SS=(Y2-B1*XY-B2*SY)/NF;S S=FSQT(SS);S SC=(Y2-BC*XY)/NF
02.80 S C1=S2/DS;S C2=-XX/DS;S C"=X2/DS;S VA=SS*C1;S SE=FSQT(VA)
02.90 S RV=1/VA;S RB=1/B1;S VR=(SS/B1^2)*(C"-2*C2*BR);S SV=FSQT(VR)

03.10 S RR=1/VR;S VB=VA/B1^4;S SR=FSQT(VB);S R1=1/VB
03.20 S VD=(SS/B1^4)*(C"-4*C2*BR);S SB=FSQT(VD);S RD=1/VD;S R'=1/D
03.30 S DO=D/B1;S BO=B1/D
03.40 T !!"                    PLOT DATA
03.50 T !"      S         V        X=S/(S+K)  V(CALC)      V-VCALC
03.60 F I=1,N;S XI=XP(I)/(XP(I)+D);S YY=XI*B1;S DI=Y(I)-YY;D 3.8
03.70 G 4.05
03.80 T !%8.05,XP(I),Y(I),"   "XI,YY,"   "DI

04.05 D 7
04.10 T !!!"MAXIMUM LIKELIHOOD                 CONFIDENCE INTERVAL",!
04.20 T "   ESTIMATE              -SE+               95%                99%
04.30 T !!"VMAX  ="%6.05,B1;F J=1,3;D 6.1
04.32 T !"1/VMAX="RB;F J=1,3;T "   "1/U(J),1/L(J)
04.34 T !!"K     ="D;F J=1,3;D 6.3;D 6.35
04.35 T !"1/K   =",R';F J=1,3;T "   "1/U(J),1/L(J)
04.36 T !!"K/VMAX="DO;F J=1,3;D 6.6;D 6.7
04.38 T !"VMAX/K="BO;F J=1,3;T "   "1/U(J),1/L(J)
04.40 T !!!!"     ANALYSIS OF VARIANCE
04.50 T !"VAR OF V =",SS,"    SE ="FSQT(SS)
04.60 T !!"VARIANCE    VMAX ="VA,"  SE ="SE,"  1/VAR ="RV
04.70 T !"ASYM VAR  1/VMAX ="VB,"  SE ="SR,"  1/VAR ="R1
04.80 T !"ASYM VAR  K      ="VR,"  SE ="SV,"  1/VAR ="RR
04.90 T !"ASYM VAR  K/VMAX ="VD,"  SE ="SB,"  1/VAR ="RD

05.50 T !!!!!
05.60 QUIT

06.05 T !%3,J,%5.1
06.10 S S=T(J)*SE;S L(J)=B1-S;S U(J)=B1+S;T "   "L(J),U(J)
06.30 S C=B1^2/(B1^2-C1*SS*T(J)^2);S K=(C-1)*C2/C1;S XD=D-C*BR-K
06.35 S S=(C-1)*(C*(BR-C2/C1)^2+1/C1/S2);I (0-S)6.4;G 6.36
06.36 S L(J)=99E99;S U(J)=L(J);T "    [- INFINITY +] "
06.40 S S=FSQT(S);S L(J)=XD-S;S U(J)=XD+S;D 6.5
06.50 T "   "L(J),U(J)
06.60 S C=B1^2-4*C1*SS*T(J)^2;S W=DO+(B1*B2-2*C2*SS*T(J)^2)/B1/C
06.70 S S=SS*(C*C"-4*C2*(B1*B2-C2*SS*T(J)^2));I (0-S)6.75;D 6.36
06.75 S S=T(J)*FSQT(S)/B1/C;D 6.8
06.77 T " [ -INFINITY+ ] "
06.80 S L(J)=W-S;S U(J)=W+S;T "   "L(J),U(J)

07.01 C------STUDENT'T VALUE TABLE (1-20) DF
07.10 S V=NF;D 8;S T(1)=1;S T(2)=A;S T(3)=B

08.10 I (1-V)8.11;S A=12.706;S B=63.657;R
08.11 I (2-V)8.12;S A=4.303;S B=9.925;R
08.12 I (3-V)8.13;S A=3.182;S B=5.841;R
08.13 I (4-V)8.14;S A=2.776;S B=4.604;R
08.14 I (5-V)8.15;S A=2.571;S B=4.032;R
08.15 I (6-V)8.16;S A=2.447;S B=3.707;R
08.16 I (7-V)8.17;S A=2.365;S B=3.499;R
08.17 I (8-V)8.18;S A=2.306;S B=3.355;R
08.18 I (9-V)8.19;S A=2.262;S B=3.25;R
08.19 I (10-V)8.2;S A=2.228;S B=3.169;R
08.20 I (11-V)8.21;S A=2.201;S B=3.106;R
08.21 I (12-V)8.22;S A=2.179;S B=3.055;R
08.22 I (13-V)8.23;S A=2.16;S B=3.012;R
08.23 I (14-V)8.24;S A=2.145;S B=2.977;R
08.24 I (15-V)8.25;S A=2.131;S B=2.947;R
08.25 I (16-V)8.26;S A=2.12;S B=2.921;R
08.26 I (17-V)8.27;S A=2.11;S B=2.898;R
08.27 I (18-V)8.28;S A=2.101;S B=2.878;R
08.28 I (19-V)8.29;S A=2.093;S B=2.861;R
08.29 I (20-V)8.3;S A=2.086;S B=2.845;R
08.30 T !!!"ENTER 2-SIDED STUDENT T VALUES FOR"%3,V," DEGREES OF FREEDOM
08.32 A !"P=.05: "A,"  P=.01: "B

09.01 C-----TEST DATA
09.10 D 1.1;D 1.2;T !!"** TEST DATA **";S N=6;T !!"N:"%1,N,!"   S     V"
09.20 S XP(1)=.138;S Y(1)=.148;S XP(2)=.220;S Y(2)=.171
09.30 S XP(3)=.291;S Y(3)=.234;S XP(4)=.560;S Y(4)=.324
09.40 S XP(5)=.766;S Y(5)=.390;S XP(6)=1.46;S Y(6)=.493
09.50 F I=1,N;T !":"%4.03,XP(I)," :"Y(I)
09.60 G 1.4
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@H|~