Directory of image this file is from
This file as a plain text file
C-PS/8 FOCAL, 1971 1.01 C FROM CHISQR BY M.J. MC KEOWN (DECUS: FOCAL8-40) 1.02 C MODIFIED BY JOHN HUISMAN 1.03 C OMSI SOFTWARE DEVELOPMENT GROUP 1.04 C LAST CHANGE: 8/19/71 1.05 1.06 C THE PROGRAM WILL PRINT THE DATA MATRIX WITH TWO 1.07 C VALUES IN EACH CELL. ONE IS THE OBSERVED VALUE AND 1.08 C THE OTHER IS THE EXPECTED VALUE AS CALCULATED. YOU MAY HAVE UP TO 1.09 C 10 DIGIT PRECISION, BUT IT IS VARIABLE. YOU MAY ENTER A 1.10 C MATRIX NO LARGER THAN 25 ELEMENTS. THE WIDTH THAT THE MATRIX 1.11 C MAY BE DEPENDS ON THE PRECISION YOU CHOOSE, AND THE WIDTH 1.12 C OF YOUR OUTPUT DEVICE. 1.13 1.14 C ** RS IS THE ROW SUM. 1.15 C ** CS IS THE COLUMN SUM. 1.16 1.17 I (FITR(/5)-(/5)) 1.18,1.3,1.18 1.18 T "YOU DO NOT HAVE ENOUGH FEATURES. " 1.19 T "RESTART FOCAL WITH FEATURE 3.",!;Q 1.20 1.30 S X=FOUT(24)+FOUT(1);F K=1,100;E 1.31 T "HOW MANY PLACES ACCURATE ";A AC;S X1=(100*AC+3)/100;T %(X1) 1.32 A "ROWS:"R,"COLUMNS:"C;S DF=(R-1)*(C-1);S L=R*C 1.33 S K=0;F I=1,R;S RS(I-1)=RP;S RP=0;D 1.4 1.34 S RS(I-1)=RP 1.35 T %(X1),!"NO MORE INPUT!"!;G 2.01 1.40 F J=1,C;D 1.41 1.41 S K=K+1;T !,%2 "SAMPLE"I,","J,": ";G 1.42 1.42 A N(K);I (10^AC-N(K)) 1.43,1.43;S RP=RP+N(K);S CS(J)=CS(J)+N(K) 1.43 T "TOO MANY DIGITS - INPUT AGAIN OR INCREASE ACCURACY",!," ";G 1.42 2.01 F J=1,1,R;S T=T+RS(J) 2.02 S J=1;S I=1;S K=1 2.03 S E(K)=(RS(I)/T)*CS(J);S J=J+1;S K=K+1 2.04 I (C-J)2.05,2.03,2.03 2.05 S J=1;S I=I+1 2.06 I (R-I)3.01,2.03,2.03 3.01 F K=1,1,L;S X2=X2+((N(K)-E(K))^2)/E(K) 4.01 T !!;S K=1;S Y=1;S N=1;S I=1;D 10;D 8.0;T ! 4.02 F J=1,1,C;D 6.0 4.03 T "*",!,"*";F J=1,C;D 5;T "*" 4.04 T "* ","RS:"RS(N);S N=N+1;T ! 4.05 F J=1,1,C;D 7.0 4.06 T "*";T !;D 8.0;T !;S I=I+1;I (R-I)4.07,4.02,4.02 4.07 T !;F J=1,1,C;D 9.0 4.08 T %," TOTAL:"T,!!,"CHI-SQUARED:"X2 4.09 T %(X1)," DEGREES OF FREEDOM:"DF,!! 4.10 A "MORE "Z;I (Z-0YES) 4.11,1.3,4.11 4.11 Q 5.01 T :(1+J*(9+AC)) 6.01 T "* O:"N(K),;D 5.01;S K=K+1 7.01 T "* E:"E(Y),;D 5.01;S Y=Y+1 8.01 F M=1,(C*(AC+9)+1);T "*" 9.01 T " CS:"CS(J)," " 10.01 T "E = EXPECTED VALUE O = OBSERVED VALUE",!!