File STAT3.FC (FOCAL source file)

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

()0Y9SKC-PS/8 FOCAL, 1971?MT:GFS N1=0;S S1=0;S S2=0;S SQ=0;S Q1=0;S DD=0;S D=0;S Q2=0?MX4HNS SQ=SQ-S1*S2/N1;S Q1=Q1-S1^2/N1;S Q2=Q2-S2^2/N1?MYHRT !,"Y=A+BX",!,"A=    ";S B=SQ/Q1;T S2/N1-B*S1/N1?MDIHVT "SD A=",FSQT(S*(1/N1+(S1/N1)^2/Q1))," SD B= ",FSQT(S/Q
)?MY^HXT "  TOTAL VARIANCE",S,!,"T VALUE FOR B",B/(FSQT(S/Q1))?MSYJ^D 4.1;D 8.35;D 8.43;D 8.45;S N3=N1?MEFJ!S DF=BA;S T=BE?ME\J&S N(1)=FITR(N3/3+1); S N(2)=N3-2*N(1)?ME*J(D 3.7; F I=1,N(1);D 9?MFFJ*S IN(1)=S1;S IN(2)=S2;D 3.7;F I=N(1)+N(2)+1,N3;D 9?MFSJ,S
M=(S1-IN(1))/N(1)?MF]J-S BA=DF;S BE=T?MF:J.S B=((S2-IN(2))/(S1-IN(1)));D 3.9;S A=S2/N1-B*S1/N3?MQ/J0T !,"Y=A+BX",!
A  
,
  B 
;I (N1-5)2.1,2.1;?MG
J9F I=1,6;S IN(I)=0?MG/J<S A=1;S K=1;S Z=N(1);D 11?MHDKAS A=N(1)+1;S K=2;S Z=N(1)+N(2);D 11?MHTKFS A=Z+1;S K=1;S Z=N3; D 11?MH(KKS Q2=IN(2)-IN(5); S Q1=IN(1)-IN(4)?MH3KPS SQ=IN(3)-IN(6)?MIDKUS DF=N(1)*M^2*(N3-3)/(2*T^2)?MISKZS I=((B*DF-SQ))/(DF-Q1)?MI)K
S J=(DF*(B^2*Q1+Q2-2*B*SQ)-Q1*Q2-SQ^2)?MI=K I (J)5.98,5.98;S J=FSQT(J)/(DF-Q1)?MJ]K!T %8.03, !, "UPPER LIMIT ",I+J," LOWER LIMIT ",I-J; G 5.5?MJ4K"T ! "IMPOSSIBLE,TOO LARGE ERRORS";G 5.5?MKGLJS X=X(100*BA+I); S Y=X(100*BE+I)?MKTLOI (X-SA)6.3,6.9,6.3?MK!L^
 (Y-SA)6.4,6.9,6.4?MK)L(S N1=N1+1?MK8L2S X(N1)=X;S X(N1-100)=Y?MQZMZR?MLIPS N1=0;F I=1,N; D 6?MLSP+S VT=0;S VP=-1?ML\P-F I=1,N;D 10?ML/RJS X=X(BA*100+I);S Y=X(BE*100+I)?ML<RTI (X-SA)9.3,9.9,9.3?MMIR^I (Y-SA)9.4,9.9,9.4?MM(R(S S2=S2+Y;D 8.2;S Q2=Q2+Y^2;S S
=SQ+X*Y;S DD=DD+(X-Y)^2?MM0R2S D=D+X-Y?MM4SZR?MM:TES J=N1?MNNTJS X1=X(100*VT+I);S X2=X(100*VT+J)?MN]TTI (X1-X2)10.4,10.3,10.4?MOGT^S X(100*VT+I)=X2;S X(100*VT+J)=X1?MO	T(S J=J-1; I (J-1)10.9,10.9,10.1?MP\UZR?MN1TS X(100*VP+I)=X(100*VP+I);S X(100*VP+J)=X(1
0*VP+J);D 10.6?MOCT<S X1=X(100*VP+I);S X(100*VP+I)=X(100*VP+J);S X(100*VP+J)=X1?MP(VOS IN(1)=Q1+IN(1);S IN(2)=IN(2)+Q2?MPHVJD 3.7;F I=A,Z;D 9?MP?VTS IN(3)=SQ+IN(3);S IN(4)=IN(4)+S1^2/N(K)?MWWVYS IN(5)=IN(5)+S2^2/N(K);S IN(6)=IN(6)+S1*S2/N(K)?MK<PTS S1=S1+
; S Q1=Q1+X^2; S N1=N1+1;R?MRMJ2A !, "CONFIDENCE LIMITS 
",DF;I (DF-0YES)2.1,5.52,2.1?MGSJ4T "T FOR  ",%3.00,N3-3,"  DEGREES OF FREEDOM  ";A T,!?MVJHJT !,"VARIABLES:   ";A "INDEPENDENT  "BA,"DEPENDENT  "BE,!?MS CTC?MY5D^L R SORT 2.15?MD<J T %8.03?MS8DJA !!!!,"ANALYSIS ",LI;IF (LI-3)2.2,4.1,2.2?MSODTI (LI-4)2.3,
.3,2.3?MV4F<F BA=1,KB;T %2.00,!,BA;F BE=1,BA;DO 3.8?MUIGUS R=1;T %6.03,R?MTEF7T !,:3;S U=0;F I=1,KB;T %7.00,I?MT\GPI (BA-BE)3.92,3.85,3.92?MU-GZDO 3.7;F I=1,N;DO 9?MW%XODO 3.55?MWNX^I (BA-BE)12.40,12.35,12.40?MX.G\DO 3.90;DO 4.14;DO 4.16;S U=U+1;S Z(U)=N1
T %6.03,R?MB2HLD 3.7;D 3.90;T %8.04?MV&F?I (N-Z(I))3.65,3.64,3.65?MV,GS Z=0?MBTGAS Z=100;R?MXDF=F I=1,U;DO 3.63?MX(S U=U+1;T :4, %7.00,Z(U)?MV?XT "      **"?MUVXJT !!,"  NUMBER MATRIX"?MW=XTF BA=1,KB;T %2.00,!,BA;F BE=1,BA;DO 12.30?MU]XYGO 2.10?MVWF>I (10
-Z)12.10,12.10;T  !,%2.00," N= ",N;G 3.95?MT	F3T !,"CORRELATION MATRIX"?MR+G
G 2.1?MCMHPS R=SQ/FSQT(Q1*Q2)?MC)HTT "    B= ",B,!;S SQ=SQ^2/Q1;S Q2=Q2-SQ;S S=Q2/(N1-2)?MD(HZT !,"EXPL. SS",SQ,"  UNEXP. SS",Q2;G 2.1?MX
F2C?MM?L S STAT3?M?

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