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:>PLTS
100  DEF FNB()=INT(27.5+5*)
150  DIM A$[12],D[563],G[53]
160  DIM G$[21]
330  Z1=Q[8]
340  F=Q4
350  L=N
360  G$="SEQUENCE PLOT "
400  GOTO Q9 OF 405,405,415
405  PRINT "WANT TO PLOT ALL ROWS OF ACTIVE DATA MATRIX ";
410  GOTO 420
415  PRINT "PLOT ALL ROWS ";
420  INPUT A$
425  IF A$[1,1]#"N" THEN 500
430  PRINT "FIRST ROW ";
435  INPUT F
436  IF F#0 THEN 440
438  F=Q4
439  GOTO 498
440  IF F >= Q4 AND F<N THEN 455
445  PRINT "GIVE NUMBER LESS THAN "N" AND MORE THAN "Q4-1
450  GOTO 435
455  PRINT "LAST ROW ";
460  INPUT L
465  IF L>F AND L <= N THEN 490
470  PRINT "GIVE NUMBER LESS THAN "N+1" AND MORE THAN "F
475  GOTO 455
490  PRINT "WANT HISTOGRAM ONLY ";
492  INPUT A$
496  IF A$[1,1]="N" THEN 500
498  G$="HISTOGRAM "
500  IF X[N9+2,Z1]=0 THEN 562
510  FOR I=F TO L
550  D[I]=(X[I,Z1]-X[N9+1,Z1])/X[N9+2,Z1]
560  NEXT I
561  GOTO 2000
562  FOR I=F TO L
563  D[I]=X[N9+1,Z1]
564  NEXT I
2000  IMAGE #,"  OBS. IN ROW ",3D," AN OUTLIER; PLOTTED AS "
2002  FOR I=F TO L
2010  IF X[I,K9+2]=0 THEN 2120
2020  D[I]=FNB(D[I])
2030  IF D[I]>.9 THEN 2070
2040  D[I]=1
2050  PRINT  USING 2000;I
2060  PRINT "-5.2"
2070  IF D[I]<53.1 THEN 2120
2080  D[I]=53
2090  PRINT  USING 2000;I
2100  PRINT "5.2"
2120  NEXT I
2130  PRINT 
6005  MAT G=ZER
6006  IF Q3#1 THEN 6017
6007  A$="OF '      '"
6008  IF Z1<11 THEN 6012
6009  Z=Z1-10
6010  A$[5,10]=N$[6*Z-5,6*Z]
6011  GOTO 6018
6012  A$[5,10]=M$[6*Z1-5,6*Z1]
6014  GOTO 6018
6017  A$=" "
6018  PRINT G$" OF STANDARDIZED VALUES "A$
6019  PRINT 
6020  IF G$[1,1]="S" THEN 6032
6021  FOR I=F TO L
6022  IF X[I,K9+2]=0 THEN 6028
6024  Z=D[I]
6026  G[Z]=G[Z]+1
6028  NEXT I
6030  GOTO 6102
6032  PRINT "  ROW"
6034  FOR I=F TO L
6036  IF INT(I/5)#I/5 THEN 6045
6038  PRINT  USING "#,3DX";I
6040  GOTO 6050
6045  PRINT "   -";
6050  IF X[I,K9+2]#0 THEN 6073
6055  PRINT " DELETED"
6060  GOTO 6089
6073  Z=D[I]
6074  G[Z]=G[Z]+1
6075  IF Z >= 27 THEN 6081
6077  PRINT TAB(Z+3)"*";
6079  PRINT TAB(30)"'"
6080  GOTO 6089
6081  IF Z>27 THEN 6087
6083  PRINT TAB(30)"*"
6085  GOTO 6089
6087  PRINT TAB(30)"'";
6088  PRINT TAB(Z+3)"*"
6089  NEXT I
6090  PRINT "    -'----'----'----'----'----'----'----'----'----'----'-"
6095  PRINT "    -5   -4   -3   -2   -1    0    1    2    3    4    5"
6096  PRINT 
6100  PRINT "HISTOGRAM"
6101  PRINT 
6102  PRINT "ABS. FREQ."
6103  PRINT 
6104  Z=G[1]
6105  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 "    -'----'----'----'----'----'----'----'----'----'----'-"
6194  PRINT "    -5   -4   -3   -2   -1    0    1    2    3    4    5"
6195  PRINT 
7000  N1=0
7010  FOR I=Q4 TO N
7020  IF X[I,K9+2]=0 THEN 7040
7030  N1=N1+1
7040  NEXT I
7050  PRINT "MEAN = "X[N9+1,Z1]
7060  PRINT "STD. DEV. = "X[N9+2,Z1]
7070  PRINT "SAMPLE SIZE = "N1
7100  IF F=Q4 AND L=N THEN 9998
7105  I5=0
7110  I6=X[F,Z1]
7115  N0=X[F,K9+2]
7120  FOR I=F+1 TO L
7125  I5=I5+X[I,K9+2]*(X[I,Z1]-I6)
7130  N0=N0+X[I,K9+2]
7135  NEXT I
7137  IF N0#0 THEN 7140
7138  PRINT '10'13"NO ACTIVE OBSERVATIONS IN ROWS "F" TO "L
7139  GOTO 9998
7140  J5=I6+I5/N0
7141  PRINT '10"BUT FOR THE "N0" OBSERVATIONS PLOTTED ABOVE:"
7142  PRINT "MEAN = "J5
7143  IF N0>1 THEN 7146
7144  J7=0
7145  GOTO 7185
7146  J6=0
7150  FOR I=F TO L
7155  IF X[I,K9+2]=0 THEN 7165
7160  J6=J6+(X[I,Z1]-J5)^2
7165  NEXT I
7170  J7=SQR(J6/(N0-1))
7185  PRINT "STD. DEV. ="J7
9998  CHAIN "$IDA",150
9999  END 
