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:2OCT73
60  REM:Q1=1 FOR >CORR;=2 FOR >COVA
70  REM: L.9010,$IDA29 CHAINS TO L.8000 IF K8=0
150  DIM A[20]
200  IF Q[1]#0 THEN 250
210  CHAIN "$IDA29",210
250  IF Q[3]=0 THEN 500
260  Q[5]=1
270  Q[6]=28
280  Q[7]=500
290  CHAIN "$IDA38"
500  PRINT '10"HOW MANY VARIABLES ";
510  INPUT N0
512  IF N0 <= K THEN 520
514  PRINT '10"ONLY"K"VARIABLE(S) IN DATA MATRIX. "'10'13'7"TRY AGAIN!"
516  GOTO 500
520  GOTO Q9 OF 530,560,580
530  PRINT "GIVE COLUMN NUMBERS "'10'13"OF VARIABLES, SEPARATED BY COMMAS: ";
550  GOTO 582
560  PRINT "COL. #'S OF VARS.: ";
570  GOTO 582
580  PRINT "COL. #'S: ";
582  MAT  INPUT A[N0]
584  FOR I=1 TO N0
586  IF A[I] <= K THEN 594
588  PRINT '10"NO DATA IN COL.'S BEYOND COL."K
590  PRINT '7"TRY AGAIN!"'10'13
592  GOTO 530
594  NEXT I
595  Z1=5
596  IF Q1=2 THEN 700
600  GOTO Q9 OF 610,640,640
610  PRINT "* NUMBER OF DECIMAL PLACES WANTED IN THE PRINTOUT? ";
630  GOTO 650
640  PRINT "*  # DECIMALS = ";
650  ENTER 15,Q8,D
651  PRINT 
652  IF Q8>0 THEN 660
654  GOTO 670
660  IF (1<D) AND (D<6) THEN 695
670  PRINT "NUMBER OF DECIMAL PLACES ALLOWED"
680  PRINT "IS BETWEEN 2 AND 5 INCLUSIVE"
690  GOTO 640
695  Z1=11-D
700  PRINT 
735  Z2=1
736  Z3=Z1
740  PRINT LIN(1);TAB(9);
760  FOR I=Z2 TO Z3 MIN N0
770  GOSUB 1000
772  IF Z1>5 THEN 778
774  PRINT SPA(7);
776  GOTO 780
778  PRINT SPA(D-1);
780  NEXT I
785  PRINT 
790  FOR I=Z2 TO N0
800  GOSUB 1000
805  PRINT "  ";
810  FOR J=Z2 TO Z3 MIN I
815  IF Q1=2 THEN 894
820  GOTO D OF 830,830,850,870,890
830  PRINT  USING "#,DD.2D,2X";M[A[I],A[J]]
840  GOTO 900
850  PRINT  USING "#,DD.3D,2X";M[A[I],A[J]]
860  GOTO 900
870  PRINT  USING "#,DD.4D,2X";M[A[I],A[J]]
880  GOTO 900
890  PRINT  USING "#,DD.5D,2X";M[A[I],A[J]]
892  GOTO 900
894  PRINT  USING "#,D.5DE,2X";M[A[I],A[J]]*X[N9+2,A[I]]*X[N9+2,A[J]]
900  NEXT J
910  PRINT 
920  NEXT I
930  PRINT LIN(1)
940  IF (N0-Z3) <= 0 THEN 9998
950  PRINT TAB(9);
960  Z2=Z2+Z1
962  Z3=Z3+Z1
970  GOTO 740
1000  IF Q3=1 THEN 1040
1010  PRINT  USING 1020;A[I]
1020  IMAGE #,"COL.",DD
1030  GOTO 1090
1040  IF A[I]>10 THEN 1080
1050  PRINT  USING 1060;M$[6*A[I]-5,6*A[I]]
1060  IMAGE #,6A
1070  GOTO 1090
1080  PRINT  USING 1060;N$[6*(A[I]-10)-5,6*(A[I]-10)]
1090  RETURN 
8000  REM:IF K8=0 IN $IDA29,>REGR
8002  D=Q[6]
8005  PRINT 
8010  PRINT  USING 8060;X[N9+1,D]
8018  N0=0
8020  FOR I=Q4 TO N
8022  IF X[I,K9+2]=0 THEN 8026
8024  N0=N0+X[I,K9+2]
8026  NEXT I
8030  PRINT  USING 8070;X[N9+2,D]/SQR(N0)
8040  PRINT  USING 8080;X[N9+2,D]
8050  PRINT  USING 8090;X[N9+2,D]*SQR(1+1/N0)
8060  IMAGE "MEAN                 = ",D.5DE
8070  IMAGE "STD. ERROR OF MEAN   = ",D.5DE
8080  IMAGE "STD. DEV.            = ",D.5DE
8090  IMAGE "STD. ERROR PRED. OBS.= ",D.5DE
9998  CHAIN "$IDA",150
9999  END 
