File CHISQR.FC (FOCAL source file)

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

E()0R(Q$C-PS/8 FOCAL, 1971?MGDAF J=1,1,R;S T=T+RS(J)?MB8DCS E(K)=(RS(I)/T)*CS(J);S J=J+1;S K=K+1?MCFDDI (C-J)2.05,2.03,2.03?MNRDES J=1;S I=I+1?MD/FAF K=1,1,L;S X2=X2+((N(K)-E(K))^2)/E(K)?MH7HBF J=1,1,C;D 6.0?MC<HEF J=1,1,C;D 7.0?MDXHFT "*";T !;D 8.0;T !;S
I=I+1;I (R-I)4.07,4.02,4.02?MM%HGT !;F J=1,1,C;D 9.0?MQLB(F J=1,C;D 1.41?MC&HAT !!;S K=1;S Y=1;S N=1;S I=1;D 10;D 8.0;T !?MGCB A "ROWS:"R,"COLUMNS:"C;S DF=(R-1)*(C-1);S L=R*C?ME6B^S X=FOUT(24)+FOUT(1);F K=1,100;E?MEHB
T "HOW MANY PLACES ACCURATE
";A AC;S X1=(100*AC+3)/100;T %(X1)?MF*BBC     MODIFIED BY JOHN HUISMAN?MRBBCC     OMSI SOFTWARE DEVELOPMENT GROUP?MI8BE?MGZB!S K=0;F I=1,R;S RS(I-1)=RP;S RP=0;D 1.4?MM>B"S RS(I-1)=RP?MB"DBS J=1;S I=1;S K=1?MTAT "E = EXPECTED VALUE     O = OBSERVED VALUE",
!?MG/RAT " CS:"CS(J)," "?MC1HDT "*  ","RS:"RS(N);S N=N+1;T !?MHJPAF M=1,(C*(AC+9)+1);T "*"?MHVHCT "*",!,"*";F J=1,C;D 5;T "*"?MITJAT :(1+J*(9+AC))?MI&LAT "*  O:"N(K),;D 5.01;S K=K+1?MH(NAT "*  E:"E(Y),;D 5.01;S Y=Y+1?MJVBFC     THE PROGRAM WILL PRINT THE 
ATA MATRIX WITH TWO ?MJ4BGC VALUES IN EACH CELL.  ONE IS THE OBSERVED VALUE AND?MMEBHC THE OTHER IS THE EXPECTED VALUE AS CALCULATED.  YOU MAY HAVE UP TO?MLIBLC OF YOUR OUTPUT DEVICE.?MKYBKC MAY BE DEPENDS ON THE PRECISION YOU CHOOSE, AND THE WIDTH?MLLBM?
L[BNC ** RS IS THE ROW SUM.?ML+BOC ** CS IS THE COLUMN SUM.?ML1BP?MEBT?MPSBQI (FITR(/5)-(/5)) 1.18,1.3,1.18?MP)BIC 10 DIGIT PRECISION, BUT IT IS VARIABLE.  YOU MAY ENTER A?MO=HHT %,"  TOTAL:"T,!!,"CHI-SQUARED:"X2?MIIHKQ?MD%BT %(X1),!"NO MORE INPUT!"!;G 2.
1?MCPDFI (R-I)3.01,2.03,2.03?MM:HJA "MORE
"Z;I (Z-0YES) 4.11,1.3,4.11?MOXB*A N(K);I (10^AC-N(K)) 1.43,1.43;S RP=RP+N(K);S CS(J)=CS(J)+N(K)?MBTB+T "TOO MANY DIGITS - INPUT AGAIN OR INCREASE ACCURACY",!,"
";G 1.42?MN HIT %(X1),"   DEGREES OF FREEDOM:"DF,!!?MRRBRT "YOU DO NOT HAVE ENOUGH FEATURES.  "?MK(BJC MATRIX NO LARGER THAN 25 ELEMENTS.  THE WIDTH THAT THE MATRIX?MN5B)S K=K+1;T !,%2 "SAMPLE"I,","J,": ";G 1.42?MFXBAC     FROM CHISQR BY M.J. MC KEOWN (DE
US: FOCAL8-40)?MGBDC     LAST CHANGE: 8/19/71?ML.BST "RESTART FOCAL WITH FEATURE 3.",!;Q?M?ML S CHISQR?M?M FOCAL WITH FEATURE 3.",!;Q? CHANGE: 8/9/71?      STARTING LOCATION: 1.3?MRIX?K"T "	";S D=FOUT(223);T "	 OR CTRL/U KILLS LINE"!!?MKG 1.4?M4?L D DTA3:
ILE.FD2J*I+MCP3X	FX616JWO3T[L2,&I+SCZ3X	F+X>I3L;I3EZ 3ESE"M&J7I+S* CI3X	F+X>I3ION3RGE+SCE3L;I3EN 3O;UT*M2I7J3P> 3U3TF&12/?N3TV +NCD+OCF3P
AE*M2IGT+D" .-Z66IJ/3EWS3T2 3IWN2 ?CU7N2JRI+CCA7I" 7FI2IZ/+OCR7NE7TB 3AWG"M&J7JCM6 WO3T[L*1>I7TCY3E2 3IWN2 WFE7D6JWO

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