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:10MAY74 60 REM:>SCAT 100 DEF FNA(Q)=INT(18.5+5*Q) 150 DIM C[5,35],D[450],E[450] 400 Z1=Q[8] 410 Z2=Q[9] 420 F=I3 430 L=I4 500 N1=0 510 FOR I=F TO L 520 IF X[I,K9+2]=0 THEN 560 530 N1=N1+1 540 E[N1]=(X[I,Z1]-X[N9+1,Z1])/X[N9+2,Z1] 550 D[N1]=(X[I,Z2]-X[N9+1,Z2])/X[N9+2,Z2] 560 NEXT I 7000 PRINT 7001 IMAGE #," HOR. OBS. IN ROW ",3D," AN OUTLIER; PLOTTED AS " 7002 IMAGE #," VERT. OBS. IN ROW ",3D," AN OUTLIER; PLOTTED AS " 7005 MAT C=ZER 7010 FOR I=1 TO N1 7015 D[I]=FNA(D[I]) 7020 E[I]=FNA(E[I]) 7025 IF D[I]>.9 THEN 7040 7030 D[I]=1 7032 GOSUB 8000 7035 PRINT USING 7001;J 7036 PRINT "-3.4" 7040 IF D[I]<35.1 THEN 7055 7045 D[I]=35 7048 GOSUB 8000 7050 PRINT USING 7001;J 7051 PRINT "3.4" 7055 IF E[I]>.9 THEN 7070 7060 E[I]=1 7062 GOSUB 8000 7065 PRINT USING 7002;J 7066 PRINT "-3.4" 7070 IF E[I]<35.1 THEN 7083 7075 E[I]=35 7078 GOSUB 8000 7080 PRINT USING 7002;J 7081 PRINT "3.4" 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 7160 7135 PRINT TAB(1+2*J)""; 7140 IF C[I2,J]>9.1 THEN 7155 7145 PRINT USING "#,D";C[I2,J] 7150 GOTO 7160 7155 PRINT TAB(1+2*J)"*"; 7160 NEXT J 7165 PRINT 7170 NEXT I2 7171 MAT C=ZER 7172 NEXT I1 7175 PRINT " - - '"; 7180 FOR I=1 TO 6 7182 PRINT " - - - - '"; 7184 NEXT I 7185 PRINT " - -" 7186 PRINT " -3"; 7187 FOR I=-2 TO 3 7188 PRINT " "; 7189 PRINT USING "#,2D";I 7190 NEXT I 7195 PRINT LIN(2) 7200 PRINT " ","MEAN","STD. DEV." 7210 PRINT "VERT. VAR.",X[N9+1,Z1],X[N9+2,Z1] 7220 PRINT "HORIZ. VAR.",X[N9+1,Z2],X[N9+2,Z2] 7230 PRINT "SAMPLE SIZE =";Q[11] 7240 IF F=Q4 AND L=N THEN 9998 7250 PRINT 7260 PRINT "BUT FOR THE"Q[12]"OBSERVATIONS PLOTTED ABOVE:" 7270 PRINT 7280 PRINT " ","MEAN","STD. DEV." 7290 PRINT "VERT. VAR.",Q[13],Q[15] 7300 PRINT "HORIZ. VAR.",Q[14],Q[16] 7310 PRINT "SAMPLE SIZE =";Q[12] 7320 GOTO 9998 8000 N3=0 8010 FOR J=F TO L 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