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()=INT(18.5+5*)
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+.1 AND 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 
