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:>HIST
100  DEF FNB()=INT(27.5+5*)
150  DIM G[53],D[563],C$[6],D$[12]
300  Z1=Q[8]
302  IF Q3=0 THEN 340
305  IF Z1<11 THEN 325
310  Z=Z1-10
315  C$=N$[6*Z-5,6*Z]
320  GOTO 340
325  C$=M$[6*Z1-5,6*Z1]
340  N1=0
350  FOR I=Q4 TO N
360  IF X[I,K9+2]=0 THEN 390
370  N1=N1+1
380  D[N1]=X[I,Z1]
390  NEXT I
400  M1=D[1]
410  M2=D[1]
420  FOR I=2 TO N1
430  IF M1 <= D[I] THEN 460
440  M1=D[I]
450  GOTO 480
460  IF M2 >= D[I] THEN 480
470  M2=D[I]
480  NEXT I
490  PRINT "MIN. OBS. =";M1;"MAX. OBS. =";M2
500  PRINT "MEAN =";X[N9+1,Z1];"STD. DEV. =";X[N9+2,Z1]
510  PRINT "SAMPLE SIZE=";N1
512  PRINT 
515  GOTO Q9 OF 520,520,540
520  PRINT "* GIVE MIDPOINT OF A CENTRAL INTERVAL OF HISTOGRAM: ";
530  GOTO 550
540  PRINT "* MIDPOINT: ";
550  ENTER 30,Q8,M1
555  PRINT 
560  IF Q8>0 THEN 575
561  PRINT 
562  PRINT "CHOOSE A CONVENIENT NUMBER NEAR THE MEAN AND/OR ONE NOT"
564  PRINT "TOO FAR FROM HALFWAY BETWEEN THE MIN.OBS. AND MAX.OBS."
566  PRINT 
570  GOTO 540
575  GOTO Q9 OF 580,580,600
580  PRINT "* GIVE WIDTH OF EACH CLASS INTERVAL OF HISTOGRAM: ";
590  GOTO 610
600  PRINT "* WIDTH OF INTERVAL: ";
610  ENTER 30,Q8,M2
615  PRINT 
620  IF Q8>0 THEN 699
621  PRINT 
622  PRINT "TRY A 'ROUND NUMBER' THAT IS A RELATIVELY SMALL FRACTION"
624  PRINT "(SAY ABOUT 1/5) OF THE STANDARD DEVIATION."
625  PRINT 
630  GOTO 600
699  M2=M2*5
700  FOR I=1 TO N1
710  D[I]=(D[I]-M1)/M2
720  NEXT I
2000  IMAGE #,"  OBS. IN ROW ",3D," AN OUTLIER; PLOTTED IN EXTREME "
2010  FOR I=1 TO N1
2020  D[I]=FNB(D[I])
2030  IF D[I]>.9 THEN 2070
2040  D[I]=1
2045  GOSUB 8000
2050  PRINT  USING 2000;J
2060  PRINT "LEFT INTERVAL"
2070  IF D[I]<53.1 THEN 2120
2080  D[I]=53
2085  GOSUB 8000
2090  PRINT  USING 2000;J
2100  PRINT "RIGHT INTERVAL"
2120  NEXT I
2130  PRINT 
6005  MAT G=ZER
6060  FOR I=1 TO N1
6073  Z=D[I]
6074  G[Z]=G[Z]+1
6080  NEXT I
6100  PRINT "HISTOGRAM"
6101  PRINT 
6103  PRINT "ABS. FREQ."
6104  PRINT 
6105  Z=G[1]
6106  FOR I=2 TO 53
6110  IF G[I] <= G[I-1] THEN 6120
6111  IF Z>G[I] THEN 6120
6115  Z=G[I]
6120  NEXT I
6125  FOR I=Z TO 1 STEP -1
6130  IF INT(I/5) <> I/5 THEN 6145
6135  PRINT  USING "#,DDD";I
6140  GOTO 6150
6145  PRINT TAB(3)"-";
6150  FOR J=1 TO 53
6155  IF I-G[J]>.99 THEN 6170
6160  PRINT TAB(3+J)"*";
6170  NEXT J
6180  PRINT 
6190  NEXT I
6192  PRINT "    -'----'----'----'----'----'----'----'----'----'----'-"
6200  PRINT "    ";
6210  FOR I=-4 TO 4 STEP 2
6220  PRINT  USING "#,2X4DE";M1+I*M2
6230  NEXT I
6233  PRINT LIN(1)
6234  IF Q3#1 THEN 6240
6235  PRINT TAB(27);C$
6240  PRINT 
6250  PRINT "MEAN =";X[N9+1,Z1]
6260  PRINT "STD. DEV. =";X[N9+2,Z1]
6270  PRINT "SAMPLE SIZE =";N1
6280  PRINT 
6290  PRINT "WANT CDF ALSO";
6295  INPUT D$
6300  IF D$[1,1]="Y" THEN 6320
6310  GOTO 9998
6320  FOR I=2 TO 53
6330  G[I]=G[I]+G[I-1]
6340  NEXT I
6350  PRINT 
6360  PRINT "CUMULATIVE"
6370  PRINT "ABS. FREQ."
6375  PRINT 
6380  FOR I=N1 TO 0 STEP -1
6383  IF INT(I/2)#I/2 THEN 6400
6385  D$="-"
6390  PRINT  USING "#,3D,A";I,D$
6400  FOR J=1 TO 53
6410  IF G[J]<I THEN 6440
6420  IF G[J]>I THEN 6445
6425  IF INT(I/2)#I/2 THEN 6435
6430  PRINT TAB(J+3)".";
6432  GOTO 6440
6435  PRINT TAB(J+3)"'";
6440  NEXT J
6445  IF INT(I/2)#I/2 THEN 6455
6450  PRINT 
6452  GOTO 6460
6455  PRINT  USING 6456
6456  IMAGE +,""
6460  NEXT I
6470  PRINT "    -'----'----'----'----'----'----'----'----'----'----'-"
6480  PRINT "    ";
6490  FOR I=-4 TO 4 STEP 2
6500  PRINT  USING "#,4X2DE";M1+I*M2
6510  NEXT I
6520  PRINT LIN(1)
6525  IF Q3#1 THEN 9998
6530  PRINT TAB(27);C$
7999  GOTO 9998
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 
