10 COM X[103,22],M[19,19],U[19,19],Q[20],V[20],P[78] 30 COM M$[60],N$[72] 40 COM N,K,N8,K8,N9,K9,Q9,Q7,Q5,Q4,Q3,Q2,Q1 45 COM I3,I4,U9,X$[20] 50 REM:12MAY74 60 REM:>NORM 100 DEF FNA(Q)=INT(18.5+5*Q) 155 DIM C[5,35],H[35],D[375],E[375] 330 Z1=Q[8] 400 X5=Q[9] 500 N1=0 508 X3=X4=0 510 FOR I=Q4 TO N 520 IF X[I,K9+2]=0 THEN 560 530 N1=N1+1 535 D[N1]=X[I,Z1]-X[N9+1,Z1] 540 X3=X3+D[N1]^3 545 X4=X4+D[N1]^4 550 D[N1]=D[N1]/X[N9+2,Z1] 560 NEXT I 570 F=0 2000 IMAGE #," OBS. IN ROW ",3D," AN OUTLIER; PLOTTED AS " 2010 MAT H=ZER 2013 PRINT 2015 FOR I=1 TO N1 2020 D[I]=FNA(D[I]) 2025 IF D[I]>.9 THEN 2045 2030 D[I]=1 2032 GOSUB 8000 2035 PRINT USING 2000;J 2040 PRINT "-3.4" 2045 IF D[I]<35.1 THEN 2060 2050 D[I]=35 2052 GOSUB 8000 2055 PRINT USING 2000;J 2056 PRINT "3.4" 2060 H[D[I]]=H[D[I]]+1 2065 NEXT I 2067 Z5=1 2070 FOR I=1 TO 35 2075 IF H[I] <= 0 THEN 2095 2080 FOR J=1 TO H[I] 2085 D[Z5]=(I-18)/5 2087 Z5=Z5+1 2090 NEXT J 2095 NEXT I 2097 E[INT(N1/2+1)]=0 2100 FOR I=INT((N1+3)/2) TO N1 2110 P=(N1-I+1)/(N1+1) 2120 T=SQR(LOG(1/(P*P))) 2130 T1=T-(2.30753+.27061*T)/(1+.99229*T+.0448*T*T) 2140 E[I]=T1 2150 E[N1-I+1]=-T1 2160 NEXT I 7005 MAT C=ZER 7010 FOR I=1 TO N1 7015 D[I]=FNA(D[I]) 7020 E[I]=FNA(E[I]) 7083 NEXT I 7084 PRINT 7085 FOR I1=7 TO 1 STEP -1 7086 FOR I=1 TO N1 7087 IF E[I]I1*5-5 THEN 7090 7088 GOTO 7092 7090 C[I1*5-E[I]+1,D[I]]=C[I1*5-E[I]+1,D[I]]+1 7092 NEXT I 7094 FOR I2=1 TO 5 7095 I=I1*5-I2+1 7100 IF ABS(INT((I-18)/5)-(I-18)/5)>.01 THEN 7120 7105 PRINT USING "#,DD";(I-18)/5 7115 GOTO 7125 7120 PRINT " -"; 7125 FOR J=1 TO 35 7130 IF C[I2,J]<.01 THEN 7157 7135 PRINT TAB(1+2*J)""; 7140 IF C[I2,J]>9.1 THEN 7146 7142 PRINT USING "#,D";C[I2,J] 7144 GOTO 7160 7146 PRINT TAB(1+2*J)"*"; 7156 GOTO 7160 7157 IF J <> I THEN 7160 7158 PRINT TAB(1+2*J)"+"; 7160 NEXT J 7165 PRINT 7170 NEXT I2 7171 MAT C=ZER 7172 NEXT I1 7175 GOSUB 7400 7180 GOSUB 7500 7190 GOSUB 7400 7191 IF F=0 THEN 7300 7192 PRINT '10'10"NOTE: FREQUENCIES OVER 9 INDICATED BY '*'" 7300 PRINT '10'10"MEAN = "X[N9+1,Z1] 7320 PRINT "STD.DEV. = "X[N9+2,Z1] 7330 X3=X3/(N1*X[N9+2,Z1]^3) 7340 X4=X4/(N1*X[N9+2,Z1]^4)-3 7350 PRINT "SKEWNESS = "X3 7360 PRINT "KURTOSIS = "X4 7365 PRINT "STUDENTIZED RANGE = "X5 7370 PRINT "SAMPLE SIZE ="N1 7390 GOTO 9998 7400 PRINT " - - '"; 7405 FOR I=1 TO 6 7410 PRINT " - - - - '"; 7415 NEXT I 7420 PRINT " - -" 7425 PRINT " -3"; 7430 FOR I=-2 TO 3 7435 PRINT " "; 7440 PRINT USING "#,2D";I 7445 NEXT I 7450 RETURN 7500 PRINT LIN(3);TAB(25);"FREQUENCY DISTRIBUTION"'13'10'10" "; 7510 FOR I=1 TO 35 7520 IF H[I]>0 AND H[I]<9.1 THEN 7580 7530 IF H[I]>9 THEN 7560 7540 PRINT " "; 7550 GOTO 7590 7560 PRINT " *"; 7562 F=1 7570 GOTO 7590 7580 PRINT USING "#,XD";H[I] 7590 NEXT I 7600 RETURN 8000 N3=0 8010 FOR J=Q4 TO N 8020 IF X[J,K9+2]=0 THEN 8050 8030 N3=N3+1 8040 IF N3=I THEN 8060 8050 NEXT J 8060 RETURN 9998 CHAIN "$IDA",150 9999 END