$EXECU IBJOB IBJOB VERSION 5 HAS CONTROL. | $IBJOB GO,MAP | $IBFTC REGT REGT000! 1 12/01/64 PAGE 1 | REGT - EFN SOURCE STATEMENT - IFN(S) - | 0| C REGRESSION ANALYSIS WITH TESTS REGT0010 | C WRITTEN BY DR. F. GEBHARDT, DEUTSCHES RECHENZENTRUM REGT0020 C F.A.GRAYBILL, AN INTRODUCTION TO LINEAR STATISTICAL MODELS, REGT0030 | C VOL. I, CHAPTER 6 REGT0040 | REGT0050 | COMMON A(100,100),C(100),CON(100),NUN(31),TX1(4), REGT0060 | 1 BEOB,FQS,IUNE,IV,LTAPE,MAX,MIN,NBEOB,NHYP,NORT,NV,SSERR,TMIN REGT0070 1 FORMAT(I2,I5,2I1,2I2,I1,4A6) REGT0080 | 2 FORMAT (13H0VARIBLE NO.,I3,20H WILL BE ELIMINATED. ) REGT0090 | CALL DATV REGT0100 2 | KTAPE = 3 REGT0110 | C REGT0120 C READ PARAMETER CARD REGT0130 | C REGT0140 | 100 READ (5, 1)NVAR,NBEOB,NORT,NFMT,MTAPE,NAUS,IBIN,TX1 REGT0150 4 | BEOB = NBEOB REGT0160 | LTAPE=0 REGT0170 IF(NORT.NE.0) LTAPE=KTAPE REGT0180 | IF (NORT.NE.0) REWIND LTAPE REGT0190 16 | 102 NV = NVAR + 1 REGT0200 | C REGT0210 | C READ DATA REGT0220 C REGT0230 | DO 103 I=1,NV REGT0240 | DO 103 J=1,NV REGT0250 | 103 A(I,J) = 0.0 REGT0260 | CALL KOVA (NVAR,NBEOB,C,A(1,NV),A,100,NFMT,MTAPE,IBIN,LTAPE,CON) REGT0270 31 DO 110 I=1,NVAR REGT0280 | A(I,NV)=A(I,NV)*BEOB REGT0290 | A(NV,I)=A(I,NV) REGT0300 | DO 110 J=1,I REGT0310 | A(I,J)=A(I,J)*(BEOB-1.) + A(I,NV)*A(J,NV)/BEOB REGT0320 110 A(J,I)=A(I,J) REGT0330 | A(NV,NV)=BEOB REGT0340 | IF (NORT.NE.0) REWIND LTAPE REGT0350 | C REGT0360 53 | DO 200 K=1,NAUS REGT0370 C REGT0380 | C UNRESTRICTED HYPOTHESIS REGT0390 | C REGT0400 | CALL REGR (1,0,0) REGT0410 58 | IF (NHYP) 170,200,120 REGT0420 C REGT0430 | C RESTRICTED HYPOTHESIS REGT0440 | C REGT0450 | 120 KHYP=NHYP REGT0460 | DO 130 I=1,KHYP REGT0470 130 CALL REGR(2,I,0) REGT0480 66 | GO TO 200 REGT0490 | C REGT0500 | C STEPWISE REGRESSION REGT0510 | C REGT0520 170 IF (NORT.EQ.0 .OR. ABS( SSERR/FQS-1. ).LE..01) GO TO 180 REGT0530 | IEL=NUN(MAX) REGT0540 | IF (MAX.EQ.IV) GO TO 175 REGT0550 | 1 12/01/64 PAGE 2 | REGT - EFN SOURCE STATEMENT - IFN(S) - 0 MAXP=MAX+1 REGT0560 | DO 173 J=MAXP,IV REGT0570 | 173 NUN(J-1)=NUN(J) REGT0580 | 175 WRITE (6,2) IEL REGT0590 85 | CALL REGR (4,0,0) REGT0600 86 IF (IUNE-1) 200,200,170 REGT0610 | 180 DO 190 J=2,IUNE REGT0620 | IF( NHYP.EQ.-2.AND.TMIN.GT.2. .OR. NHYP.EQ.-3.AND.TMIN.GT.3.) REGT0630 | 1 GO TO 200 REGT0640 | IEL=NUN(MIN) REGT0650 190 CALL REGR (3,J-1,IEL) REGT0660 98 | 200 CONTINUE REGT0670 | GO TO 100 REGT0680 | END REGT0690 | 1 12/01/64 PAGE 3 | 0| -| $IBFTC DATV7 REGT0700 | 1 12/01/64 PAGE 4 DATV7 - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE DATV REGT0710 | C *** THIS SUBROUTINE MAY BE REPLACED BY THE USER BY A SUBROUTINE REGT0720 | C *** TO MANIPULATE THE DATA ENTERING COMPUTATION. REGT0730 RETURN REGT0740 | END REGT0750 | 1 12/01/64 PAGE 5 | | 0 -| $IBFTC ELIMI REGT0760 | 1 12/01/64 PAGE 6 | ELIMI - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE ELIM(NV,MIN,D,BETA,NUN) REGT0770 | C SUBROUTINE TO ELIMINATE ONE ROW AND COLUMN FROM A MATRIX. REGT0780 | C REGT0790 | DIMENSION D(31,31),BETA(1),NUN(1) REGT0800 | DO 10 I=1,NV REGT0810 10 BETA(I)=D(I,MIN) REGT0820 | B=BETA(MIN) REGT0830 | IF (NV.EQ.MIN) RETURN REGT0840 | MINM=MIN-1 REGT0850 | MINP=MIN+1 REGT0860 DO 30 I=MINP,NV REGT0870 | NUN(I-1)=NUN(I) REGT0880 | BETA(I-1)=BETA(I) REGT0890 | DO 20 J=1,MINM REGT0900 | D(I-1,J)=D(I,J) REGT0910 20 D(J,I-1)=D(J,I) REGT0920 | DO 30 J=MINP,I REGT0930 | D(I-1,J-1)=D(I,J) REGT0940 | 30 D(J-1,I-1)=D(J,I) REGT0950 | BETA(NV)=B REGT0960 RETURN REGT0970 | END REGT0980 | 1 12/01/64 PAGE 7 | | 0 -| $IBFTC 7KOVA REGT0990 | 1 12/01/64 PAGE 8 | 7KOVA - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE KOVA (NV,NP,X,S,C,M,IFT,IT,IB,IZ,Y) REGT1000 | C SUBROUTINE TO COMPUTE A COVARIANCE MATRIX. REGT1010 | C THE FIRST PART OF THE DATA IS READ INTO THE MATRIX C TO COMPUTE REGT1020 | C TEMPORARY MEANS,Y(I). USING THESE TEMPORARY MEANS IN COMPUTING REGT1030 | C THE COVARIANCE MATRIX REDUCES ROUNDING ERRORS. REGT1040 C MEANING OF ARGUMENTS REGT1050 | C NV = NUMBER OF VARIABLES REGT1060 | C NP = NUMBER OF PROBANDS IFT = NUMBER OF FORMAT-CARDS REGT1070 | C X = ERASABLE ARRAY, SIZE NV IT = TAPE NUMBER REGT1080 | C S = MEAN VALUE IB BINARY DATA IF UNEQUAL 0 REGT1090 C C = COVARIANCE MATRIX IZ IF IZ.GT.0 WRITE DATA ON TAPE REGT1100 | C M = DIMENSION OF C Y = ERASABLE ARRAY, SIZE NV REGT1110 | C REGT1120 | DIMENSION X(1),S(1),C(M,M),FT(108),Y(1) REGT1130 | IF(IFT.EQ.0 .AND. IB.EQ.0) JFT=12 REGT1140 IF(IFT.GT.0 .AND. IB.EQ.0) JFT=12*IFT REGT1150 | JT=5 REGT1160 | IF(IT.GT.0) JT=IT REGT1170 | IF(JFT.GT.0) READ (5,1) (FT(I),I=1,JFT) REGT1180 9 | DO 20 I=1,NV REGT1190 20 S(I)=0. REGT1200 | L=MIN0(M,NP) REGT1210 | C REGT1220 | C ** PROVISIONAL MEAN VALUES REGT1230 | C REGT1240 DO 30 K=1,L REGT1250 | IF(IB.EQ.0) READ (JT,FT) (C(K ,I),I=1,NV) REGT1260 28 | IF(IB.NE.0) READ (JT) (C(K,I),I=1,NV) REGT1270 34 | IF(IZ.GT.0) WRITE (IZ) (C(K,I),I=1,NV) REGT1280 40 | DO 30 I=1,NV REGT1290 30 S(I)=S(I)+ C(K,I) REGT1300 | AL=L REGT1310 | DO 40 I=1,NV REGT1320 | Y(I)=S(I)/AL REGT1330 | 40 S(I)=0 REGT1340 C REGT1350 | C ** CROSS PRODUCTS FROM THE DATA STORED IN C REGT1360 | C REGT1370 | DO 70 I=1,NV REGT1380 | DO 50 J=1,NV REGT1390 50 X(J)=0. REGT1400 | DO 60 K=1,L REGT1410 | S(I)=S(I) + (C(K,I)-Y(I)) REGT1420 | DO 60 J=I,NV REGT1430 | 60 X(J)=X(J) + (C(K,I)-Y(I))*(C(K,J)-Y(J)) REGT1440 DO 70 J=1,NV REGT1450 | 70 C(J,I)=X(J) REGT1460 | IF(L.EQ.NP) GO TO 85 REGT1470 | MP=M+1 REGT1480 | C REGT1490 C ** CROSS PRODUCT FROM REMAINING DATA REGT1500 | C REGT1510 | DO 80 K=MP,NP REGT1520 | IF(IB.EQ.0) READ (JT,FT) (X(I),I=1,NV) REGT1530 104 | IF(IB.NE.0) READ (JT) (X(I),I=1,NV) REGT1540 110 1 12/01/64 PAGE 9 | 7KOVA - EFN SOURCE STATEMENT - IFN(S) - | 0 IF(IZ.GT.0) WRITE (IZ) (X(I),I=1,NV) REGT1550 116 | DO 80 I=1,NV REGT1560 | S(I)=S(I) + (X(I)-Y(I)) REGT1570 DO 80 J=I,NV REGT1580 | 80 C(J,I)=C(J,I) + (X(I)-Y(I))*(X(J)-Y(J)) REGT1590 | C REGT1600 | C ** COMPUTATION OF COVARIANCES REGT1610 | C REGT1620 85 P=NP REGT1630 | DO 90 I=1,NV REGT1640 | S(I)=S(I)/P REGT1650 | DO 90 J=1,I REGT1660 | C(I,J)=(C(I,J)- P*S(I)*S(J))/(P-1.) REGT1670 90 C(J,I)=C(I,J) REGT1680 | DO 100 I=1,NV REGT1690 | 100 S(I)=S(I)+Y(I) REGT1700 | RETURN REGT1710 | 1 FORMAT (12A6) REGT1720 END REGT1730 | 1 12/01/64 PAGE 10 | | 0| - $IBFTC 7MATIN REGT1740 | 1 12/01/64 PAGE 11 | 7MATIN - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE MATIN (A,N,IN,V,U,M) REGT1750 C SUBROUTINE MATIN SEPTEMBER 1, 1967 REGT1760 | C PROGRAM WRITTEN BY HEALTH SCIENCES REGT1770 | C COMPUTING FACILITY, UCLA (BMD 04 V) REGT1780 | C ALTERED (VARIABLE DIMENSION) BY DR. F. GEBHARDT REGT1790 | C MATIN INVERTS AN N*N SYMETRIC MATRIX. REGT1800 C DETERMINANT WILL BE WRITTEN INTO V(1). REGT1810 | C MEANING OF ARGUMENTS REGT1820 | C A = MATRIX REGT1830 | C N = NUMBER OF LINES OF THE MATRIX REGT1840 | C IN,V,U = STORAGES ,SIZE N REGT1850 C M = DIMENSION OF A REGT1860 | C REGT1870 | DIMENSION A(M,M),IN(1),V(1),U(1) REGT1880 | DO 1 I=1,N REGT1890 | V(I)=A(I,I) REGT1900 1 IN(I)=0 REGT1910 | D=1.0 REGT1920 | K=1 REGT1930 | DO 7 L=1,N REGT1940 | DO 2 I=1,K REGT1950 U(I)=A(K,I) REGT1960 | 2 A(K,I)=0.0 REGT1970 | P=U(K) REGT1980 | DO 3 I=K,N REGT1990 | U(I)=A(I,K) REGT2000 3 A(I,K)=0.0 REGT2010 | T=H REGT2020 | H=-1.E20 REGT2030 | IN(K)=1 REGT2040 | U(K)=-1.0 REGT2050 D=D*P REGT2060 | DO 7 I=1,N REGT2070 | Y=U(I)/P REGT2080 | DO 4 J=1,I REGT2090 | 4 A(I,J)=A(I,J)-Y*U(J) REGT2100 IF(IN(I))5,5,7 REGT2110 | 5 IF(H-A(I,I)/V(I))6,7,7 REGT2120 | 6 H=A(I,I)/V(I) REGT2130 | K=I REGT2140 | 7 CONTINUE REGT2150 DO 8 I=1,N REGT2160 | DO 8 J=1,I REGT2170 | A(I,J)=-A(I,J) REGT2180 | 8 A(J,I)=A(I,J) REGT2190 | V(1)=D REGT2200 RETURN REGT2210 | END REGT2220 | 1 12/01/64 PAGE 12 | | 0 -| $IBFTC NORMT REGT2230 | 1 12/01/64 PAGE 13 | NORMT - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE NRMT(NV,NA,A,NRA,NB,B,NRB,NBEOB,C,NABH,KTAPE,S2) REGT2240 | C TEST ON NORMALITY REGT2250 | C WRITTEN BY DR. F. GEBHARDT, DEUTSCHES RECHENZENTRUM REGT2260 | C F. GEBHARDT, VERTEILUNG UND SIGNIFIKANZ-SCHRANKEN DES REGT2270 | C 3. UND 4. STICHPROBENMOMENTES BEI NORMALVERTEILTEN VARIABLEN, REGT2280 C BIOMETRISCHE ZEITSCHRIFT (1966) 8,219-241 REGT2290 | C REGT2300 | DIMENSION A(1),B(1),C(1),NRA(1),NRB(1) REGT2310 | DIMENSION ZW(6),ZR(6), II(6),Q(4) REGT2320 | DATA Q /1H ,6H NOT,6HWEAKLY,6HHIGHLY / REGT2330 C REGT2340 | C * SUMS OF POWERS REGT2350 | C REGT2360 | S1 = 0.0 REGT2370 | S2 = 0.0 REGT2380 S3 = 0.0 REGT2390 | S4 = .0 REGT2400 | NVAR=NV-1 REGT2410 | C(NV)=1. REGT2420 | DO 40 IJ =1,NBEOB REGT2430 READ (KTAPE) (C(I),I=1,NVAR) REGT2440 6 | Y = C(NABH) REGT2450 | IF(NA) 12,20,12 REGT2460 | 12 DO 15 I=1,NA REGT2470 | IA = NRA(I) REGT2480 15 Y = Y - C(IA)*A(I) REGT2490 | 20 IF(NB) 25,30,25 REGT2500 | 25 DO 28 I =1,NB REGT2510 | IA = NRB(I) REGT2520 | 28 Y = Y - C(IA)*B(I) REGT2530 30 S1 = S1 + Y REGT2540 | S2 = S2 + Y**2 REGT2550 | S3 = S3 + Y**3 REGT2560 | 40 S4 = S4 + Y**4 REGT2570 | REWIND KTAPE REGT2580 C REGT2590 | C * TEST STATISTICS REGT2600 | C REGT2610 34 | FN = NBEOB REGT2620 | SM1 = S1 / FN REGT2630 SM2 = S2 / FN - SM1**2 REGT2640 | SM3 = S3/FN - 3.0*S2*SM1/FN + 2.0*SM1**3 REGT2650 | SM4 = S4/FN - 4.0*S3*SM1/FN + 6.0*S2*SM1**2/FN - 3.0*SM1**4 REGT2660 | WRITE (6,2) S2 REGT2670 35 | ST3 = SM3/SM2**1.5 * SQRT( ((FN + 4.0)*FN + 3.0)/(6.0*FN - 12.0) )REGT2680 36 37 ST4 = (SM4/SM2**2 -3.0+6.0/(FN + 1.0) ) * SQRT( ((FN + 8.0)*FN + 1REGT2690 | 15.0) / (((24.0*FN - 120.0)*FN + 144.0)*FN) ) * (FN + 1.0) REGT2700 38 | U = ABS(ST3) REGT2710 | Y=SQRT(1.0/FN) REGT2720 39 | TX=Q(2) REGT2730 IF( AMIN1( ((-3.3646*Y+1.7309)*Y-.0797)*Y+1.9628 - U, ((-6.1463*Y-REGT2740 | 16.9178)*Y+4.4435)*Y+2.0056 - ST4, ((-12.6973*Y+17.6666)*Y-6.3611)*REGT2750 | 2Y+1.9508 + ST4) ) 53,60,60 REGT2760 | 53 TX=Q(3) REGT2770 | IF (AMIN1( ((-11.5468*Y+3.6195)*Y+ .3909)*Y+2.3158 - U, ((-2.7158*REGT2780 1 12/01/64 PAGE 14 | NORMT - EFN SOURCE STATEMENT - IFN(S) - | 0 1Y-21.4460)*Y+9.9154)*Y+2.3735 - ST4, ((-24.7735*Y+27.5382)*Y-9.029REGT2790 | 24)*Y+2.3133 + ST4) ) 54,60,60 REGT2800 | 54 TX=Q(1) REGT2810 IF (AMIN1( ((-26.8301*Y+6.6524)*Y+1.3764)*Y+2.7708 - U, ((8.4573*YREGT2820 | 1-53.6171)*Y+21.3461)*Y+2.7977 - ST4, ((-39.7062*Y+39.9197)*Y-12.52REGT2830 | 208)*Y+2.7815 + ST4) ) 55,60,60 REGT2840 | 55 TX=Q(4) REGT2850 | 60 WRITE (6,1) ST3,ST4,TX REGT2860 46 RETURN REGT2870 | 1 FORMAT (1H0/40H0TEST OF RESIDUES ON NORMAL DISTRIBUTION / REGT2880 | 2 49H TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = F8.3/ REGT2890 | 1 49H TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = F8.3/ REGT2900 | 2 10H0DEVIATION REGT2910 3 40H OF RESIDUES FROM NORMAL DISTRIBUTION IS ,A6,12H SIGNIFICANT )REGT2920 | 2 FORMAT (1H0//44H0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM, REGT2930 | 1 1PE15.5) REGT2940 | END REGT2950 | 1 12/01/64 PAGE 15 | 0| -| $IBFTC REDU REGT2960 | 1 12/01/64 PAGE 16 REDU - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE REDUX(D,B,M,N) REGT2970 | C SUBROUTINE TO COMPUTE THE INVERSE OF A MATIX IF THE INVERSE OF REGT2980 | C THE MATRIX BORDERED BY AN ADDITIONAL ROW AND COLUMN IS GIVEN. REGT2990 C REGT3000 | DIMENSION D(M,M),B(1) REGT3010 | R=1./B(N) REGT3020 | L=N-1 REGT3030 | DO 10 I=1,L REGT3040 DO 10 J=1,L REGT3050 | 10 D(I,J)= D(I,J) - B(I)*R*B(J) REGT3060 | RETURN REGT3070 | END REGT3080 | 1 12/01/64 PAGE 17 | 0| -| $IBFTC REGRE REGT3090 | 1 12/01/64 PAGE 18 REGRE - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE REGR(ITYPE,K1,IEL) REGT3100 | C REGRESSION SUBROUTINE FOR REGT. REGT3110 | C WRITTEN BY DR. F. GEBHARDT, DEUTSCHES RECHENZENTRUM. REGT3120 COMMON A(100,100),C(100),CON(100),NUN(31),TX1(4), REGT3130 | 1 BEOB,FQS,IUNE,IV,LTAPE,MAX,MIN,NBEOB,NHYP,NORT,NV,SSERR,TMIN REGT3140 | DIMENSION BETA(31),D(31,31),KONV(100),M(100),TX(4,2),Z(31) REGT3150 | 2 FORMAT(12A6) REGT3160 | 3 FORMAT(I2,3I1,34I2) REGT3170 4 FORMAT( 7(I2,F8.0) ) REGT3180 | 5 FORMAT(1H1,50X,19HREGRESSION ANALYSIS// REGT3190 | 1 24H DEUTSCHES RECHENZENTRUM,84X, 4A6/ REGT3200 | 2 10H DARMSTADT, 98X, 4A6 // REGT3210 | 3 24H UNRESTRICTED HYPOTHESIS / REGT3220 4 10X,18HDEPENDENT VARIABLE , I11/ REGT3230 | 5 10X,21HINDEPENDENT VARIABLES,4X,24I4/35X,6I4) REGT3240 | 6 FORMAT (30H0REGRESSION THROUGH THE ORIGIN ) REGT3250 | 7 FORMAT (41H0ABSOLUTE TERM IS SPECIFIED (VAR. NO. 0) ) REGT3260 | 8 FORMAT (37H0COEFFICIENTS SPECIFIED IN HYPOTHESIS / REGT3270 1 6(12H NO VALUE, 8X) / 6(I3, 1PE13.5, 4X) ) REGT3280 | 9 FORMAT (1H0/1H0/40X,31HINVERSE OF CROSS PRODUCT MATRIX //) REGT3290 | 10 FORMAT(I4, 1P 7E16.6 / (4X 7E16.6) ) REGT3300 | 11 FORMAT (1H0/ 23H0NUMBER OF OBSERVATIONS ,13X,I6/ REGT3310 | 1 21H TOTAL SUM OF SQUARES ,17X,1PE14.5,I8,5H D.F. ) REGT3320 12 FORMAT (21H ERROR SUM OF SQUARES, 17X,1PE14.5,I8,5H D.F. / REGT3330 | 1 18H MEAN SQUARE ERROR,15X E19.5 //) REGT3340 | 13 FORMAT (1H0/63H0VAR.NO. MEAN STAND.DEV. REGRESS. REGT3350 | 1STAND.DEV.OF6X,1HT,5X,4HD.F.,4X,25HPART. ERROR PROB. MULT. REGT3360 | 2 5X,5HCORR./25X,7HOF MEAN,8X,6HCOEFF.6X10HREG.COEFF.,21X, REGT3370 3 5HCORR., 3X, 9H(2-SIDED),3X,5HCORR.,5X,6HWITH Y //) REGT3380 | 14 FORMAT (I6,2F13.5,2F14.5, F12.3,I6,4F10.4 ) REGT3390 | 15 FORMAT ( 68H0SUBSTANTIAL DEVIATION OF THE ERROR SUM OF SQUARES FROREGT3400 | 1M THE CHECKSUM /70H PRESUMABLY INDICATES THAT THE COVARIANCE MATRIREGT3410 | 2X IS (NEARLY) SINGULAR. ) REGT3420 16 FORMAT(1H1,107X,4A6/ REGT3430 | 1 I5,24H. RESTRICTED HYPOTHESIS ,79X,4A6/ REGT3440 | 2 10X,18HDEPENDENT VARIABLE,I11,69X,4A6/ REGT3450 | 3 10X,21HINDEPENDENT VARIABLES,4X,24I4/35X,6I4) REGT3460 | 17 FORMAT (1H0/1H0/47H ANALYSIS OF VARIANCE FOR TESTING THE RESTRICTEREGT3470 113HD HYPOTHESIS // REGT3480 | 2 7H SOURCE,20X,10HSUM OF SQ.,4X,4HD.F.,4X,8HMEAN SQ.,6X,1HF// REGT3490 | 3 6H TOTAL,1PE32.5,I7,E14.5/ REGT3500 | 4 20H UNRESTR. HYPOTHESIS,1PE18.5, I7,E14.5/ REGT3510 | 5 20H RESTR. HYPOTHESIS 1PE18.5, I7,E14.5/ REGT3520 6 11H DIFFERENCE ,1PE 27.5,I7,E14.5,0PF10.3 / REGT3530 | 7 6H ERROR, 1PE32.5, I7, E14.5) REGT3540 | 18 FORMAT (14H NO REGRESSION ) REGT3550 | 19 FORMAT (1X,I5,26X,2F14.5,F12.3, I6,F20.4 ) REGT3560 | 20 FORMAT (38H MEAN VALUE OF THE DEPENDENT VARIABLE , 1PE14.5/ REGT3570 1 32H STANDARD-DEV. OF THE MEAN VALUE ,E20.5 ) REGT3580 | 21 FORMAT(1H0/1H0) REGT3590 | 22 FORMAT(21H MULTIPLE CORRELATION, 17X, F10.5/ REGT3600 | 1 36H SQUARE OF THE MULTIPLE CORRELATION F12.5,7X,2HF(,I2,1H,, REGT3610 | 2 I4,3H) =, F8.3) REGT3620 23 FORMAT (/72H THESE RESULTS REFER TO THE LINEAR COMBINATION OF THE REGT3630 | 1DEPENDENT VARIABLE/73H WITH THOSE VARIABLES WHOSE COEFFICIENTS AREREGT3640 | 1 12/01/64 PAGE 19 | REGRE - EFN SOURCE STATEMENT - IFN(S) - | 0 2 SPECIFIED BY THE HYPOTHESIS. ) REGT3650 25 FORMAT (34H0STEPWISE REGRESSION. VARIABLE NO., I3, REGT3660 | 1 21H HAS BEEN ELIMINATED. ) REGT3670 | 26 FORMAT (54H0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE/REGT3680 | 1 57H RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES.) REGT3690 | DATA QHL /1H / REGT3700 ITYP=ITYPE REGT3710 | DO 30 I=1,4 REGT3720 | 30 TX(I,2)=QHL REGT3730 | ITP=0 REGT3740 | IF (ITYP-3) 60,50,40 REGT3750 40 ITYP=1 REGT3760 | IUNE=IUNE-1 REGT3770 | ITP=1 REGT3780 | ITY=2 REGT3790 | GO TO 110 REGT3800 50 IUNE=IUNE-1 REGT3810 | ITP=1 REGT3820 | CALL ELIM (IV,MIN,D,BETA,NUN) REGT3830 20 | ITY=2 REGT3840 | IF( IUNE ) 123,122,123 REGT3850 C REGT3860 | C * SELECTION CARD REGT3870 | C REGT3880 | 60 READ (5,3) NHYP,ITX,INV,NULLP,KONS,KABH,IUNE,(NUN(I),I=1,IUNE) REGT3890 24 | IF (ITYP-1) 90,70,90 REGT3900 70 NABH=KABH REGT3910 | DO 80 I=1,4 REGT3920 | 80 TX(I,1)=QHL REGT3930 | 90 IF(ITX.NE.0) READ (5,2) (TX(I,ITYP), I=1,4) REGT3940 49 | ITY=1 REGT3950 IF (ITYP.EQ.2) GO TO 123 REGT3960 | 110 WRITE (6,5) TX1,(TX(I,1),I=1,4),NABH,(NUN(I),I=1,IUNE) REGT3970 61 | GO TO 125 REGT3980 | 122 WRITE (6,16) TX1,K1,(TX(I,1),I=1,4),NABH,(TX(I,2),I=1,4) REGT3990 76 | GO TO 124 REGT4000 123 WRITE (6,16) TX1,K1,(TX(I,1),I=1,4),NABH,(TX(I,2),I=1,4), REGT4010 | 1 (NUN(I),I=1,IUNE) REGT4020 91 | 124 IF (ITYP.EQ.3) WRITE (6,25) IEL REGT4030 | C REGT4040 | C * SPECIFIED COEFFICIENTS REGT4050 C REGT4060 112 | 125 IF (KONS.EQ.0) GO TO 135 REGT4070 | 132 IF (ITP.EQ.0) READ (5,4) (KONV(I),CON(I), I=1,KONS) REGT4080 118 | WRITE (6,8)(KONV(I),CON(I), I=1,KONS) REGT4090 127 | IF(KONV(1)) 135,137,135 REGT4100 135 IF (NULLP) 136,145,136 REGT4110 | 136 WRITE (6,6) REGT4120 138 | GO TO 140 REGT4130 | 137 WRITE (6,7) REGT4140 140 | KONV(1) = NV REGT4150 140 IV = IUNE REGT4160 | GO TO 150 REGT4170 | 145 IV = IUNE + 1 REGT4180 | NUN(IV) = NV REGT4190 | C REGT4200 1 12/01/64 PAGE 20 | REGRE - EFN SOURCE STATEMENT - IFN(S) - | 0C * CROSS PRODUCT MATRIX REGT4210 | C REGT4220 | 150 IF(ITYP.EQ.3 .AND. NORT.NE.0 .AND. RUND.LT..001) GO TO 157 REGT4230 151 DO 155 I=1,IV REGT4240 | IA = NUN(I) REGT4250 | DO 155 J=1,IV REGT4260 | JA = NUN(J) REGT4270 | 155 D(I,J) = A(IA,JA) REGT4280 CALL MATIN (D,IV,M,C,C(51),31) REGT4290 164 | GO TO 160 REGT4300 | 157 CALL REDUX(D,BETA,31,IV+1) REGT4310 | C REGT4320 | C * WRITE INVERSE MATRIX REGT4330 C REGT4340 167 | 160 IF (INV) 1605,1615,1605 REGT4350 | 1605 WRITE (6,9) REGT4360 170 | DO 161 I=1,IV REGT4370 | NA = MOD(NUN(I),NV) REGT4380 161 WRITE (6,10)NA,(D(I,J), J=1,IV) REGT4390 | C REGT4400 | C * REGRESSION COEFFICIENTS REGT4410 | C REGT4420 175 | 1615 DO 162 I=1,IV REGT4430 IA = NUN(I) REGT4440 | 162 Z(I) = A(IA,NABH) REGT4450 | IF(KONS) 163,167,163 REGT4460 | 163 DO 165 I=1,IV REGT4470 | IA = NUN(I) REGT4480 DO 165 J=1,KONS REGT4490 | JA = KONV(J) REGT4500 | 165 Z(I) = Z(I) - A(IA,JA)*CON(J) REGT4510 | 167 DO 169 I=1,IV REGT4520 | 169 BETA(I) = 0.0 REGT4530 DO 170 J=1,IV REGT4540 | DO 170 I=1,IV REGT4550 | 170 BETA(I) = BETA(I) + D(I,J)*Z(J) REGT4560 | C REGT4570 | C * SUM OF SQUARES REGT4580 C REGT4590 | SS = A(NABH,NABH) REGT4600 | IF(KONS) 175,182,175 REGT4610 | 175 DO 180 I=1,KONS REGT4620 | IA = KONV(I) REGT4630 SS = SS - 2.0*A(NABH,IA)*CON(I) REGT4640 | DO 180 J=1,KONS REGT4650 | JA = KONV(J) REGT4660 | 180 SS = SS + CON(J)*A(JA,IA)*CON(I) REGT4670 | 182 SSERR = SS REGT4680 IF(IV) 184,192,184 REGT4690 | 184 DO 185 I=1,IV REGT4700 | 185 SSERR = SSERR - BETA(I)*Z(I) REGT4710 | NFG = NBEOB - IV REGT4720 | FG=NFG REGT4730 SSERRM = SSERR / FLOAT(NFG) REGT4740 | C REGT4750 | C * WRITE REGRESSION COEFFICIENTS REGT4760 | 1 12/01/64 PAGE 21 | REGRE - EFN SOURCE STATEMENT - IFN(S) - 0C REGT4770 | SSW = SQRT(SSERRM) REGT4780 263 | WRITE (6,13) REGT4790 264 | TMIN=1.E20 REGT4800 | SMAX=0. REGT4810 CORY=A(NABH,NABH)*BEOB - A(NABH,NV)**2 REGT4820 | DO 190 I=1,IV REGT4830 | NA= MOD(NUN(I),NV) REGT4840 | ERRB = SQRT( D(I,I) )*SSW REGT4850 274 | T=BETA(I) / ERRB REGT4860 CALL TVRT(T,NFG,P) REGT4870 276 | P=1.0- 2.0*ABS(.5-P) REGT4880 | IF(NA) 186,189,186 REGT4890 | 186 XQ=A(NA,NV)/BEOB REGT4900 | STXQ=SQRT((A(NA,NA)/BEOB-XQ**2)/(BEOB-1.0) ) REGT4910 282 RM=0. REGT4920 | IF (IUNE.LE.1 .OR. IV.EQ.IUNE) GO TO 1865 REGT4930 | RM= SQRT( 1. - 1./(D(I,I)*STXQ**2*(BEOB-1.)*BEOB) ) REGT4940 287 | 1865 RPART = T/SQRT(FG+T**2) REGT4950 289 | CORR = (A(NA,NABH)*BEOB - A(NA,NV)*A(NABH,NV)) / REGT4960 1 SQRT( (A(NA,NA)*BEOB - A(NA,NV)**2) * CORY ) REGT4970 296 | WRITE (6,14) NA,XQ,STXQ, BETA(I),ERRB,T,NFG, RPART,P,RM,CORR REGT4980 297 | IF (IV.EQ.IUNE) RM=A(NA,NA)*D(I,I) REGT4990 | IF (ABS(T).GE.TMIN) GO TO 187 REGT5000 | TMIN=ABS(T) REGT5010 MIN=I REGT5020 | 187 IF (RM .LT. SMAX) GO TO 190 REGT5030 | SMAX=RM REGT5040 | MAX=I REGT5050 | GO TO 190 REGT5060 189 WRITE (6,19)NA,BETA(I),ERRB,T,NFG,P REGT5070 310 | 190 CONTINUE REGT5080 | WRITE (6,26) REGT5090 | C REGT5100 | C * WRITE SUMS OF SQUARES REGT5110 C REGT5120 314 | 192 SSM = SS/FLOAT(NBEOB) REGT5130 | SSH=SS-SSERR REGT5140 | SSHM=SSH / FLOAT(IV) REGT5150 | IF (IV.EQ.0) SSHM=0. REGT5160 YQ=A(NABH,NV) REGT5170 | IF(KONS) 1900,1902,1900 REGT5180 | 1900 DO 1901 I=1,KONS REGT5190 | IA=KONV(I) REGT5200 | 1901 YQ=YQ-CON(I)*A(IA,NV) REGT5210 1902 YQ=YQ/BEOB REGT5220 | SSY=(SS/BEOB-YQ**2)/(BEOB-1.0) REGT5230 | SSYSQ=SQRT(SSY) REGT5240 331 | IF(IV-IUNE) 1903,1904,1903 REGT5250 | 1903 CMULTQ=1.0-SSERR/(SS-YQ**2*BEOB) REGT5260 CMULT=SIGN(SQRT(ABS(CMULTQ)),CMULTQ) REGT5270 334 | F=(SS-YQ**2*BEOB-SSERR)*FLOAT(NFG) / (SSERR*FLOAT(IUNE)) REGT5280 | 1904 IF (ITYP-2) 200,230,230 REGT5290 | 200 WRITE (6,11) NBEOB,SS,NBEOB REGT5300 337 | WRITE (6,12) SSERR,NFG,SSERRM REGT5310 338 WRITE (6,20) YQ,SSYSQ REGT5320 339 | 1 12/01/64 PAGE 22 | REGRE - EFN SOURCE STATEMENT - IFN(S) - | 0 IF(KONS.NE.0) WRITE (6,23) REGT5330 340 | IF (IV.NE.IUNE) WRITE (6,22) CMULT,CMULTQ,IUNE,NFG,F REGT5340 342 SSERRU=SSERR REGT5350 | SSERUM=SSERRM REGT5360 | NFGU=NFG REGT5370 | IVU=IV REGT5380 | GO TO 250 REGT5390 230 SSU=SS-SSERRU REGT5400 | SSUM = SSU / FLOAT(IVU) REGT5410 | SSD = SSU - SSH REGT5420 | NFGD= IVU - IV REGT5430 | SSDM= SSD / FLOAT(NFGD) REGT5440 FD=SSDM / SSERUM REGT5450 | WRITE (6,21) REGT5460 348 | WRITE (6,20)YQ,SSYSQ REGT5470 349 | IF (IV.NE.IUNE) WRITE (6,22) CMULT,CMULTQ,IUNE,NFG,F REGT5480 350 | WRITE (6,17) SS,NBEOB,SSM, SSU,IVU,SSUM, SSH,IV,SSHM, REGT5490 1 SSD,NFGD,SSDM,FD, SSERRU,NFGU,SSERUM REGT5500 351 | IF (KONS.NE.0) WRITE (6,23) REGT5510 352 | 250 IF (NORT.EQ.0) RETURN REGT5520 | CALL NRMT(NV,IV,BETA,NUN,KONS,CON,KONV,NBEOB,C, REGT5530 | 1 NABH,LTAPE,FQS) REGT5540 357 RUND=ABS(1.-SSERR/FQS) REGT5550 | IF(RUND.GE..01) WRITE (6,15) REGT5560 359 | RETURN REGT5570 | END REGT5580 | 1 12/01/64 PAGE 23 | 0| -| $IBFTC 7TVRT REGT5590 | 1 12/01/64 PAGE 24 7TVRT - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE TVRT(T,N,P) REGT5600 | C TVRT COMPUTES THE PROBABILITY P THAT A RANDOM VARIABLE DISTRIBUTEDREGT5610 | C LIKE STUDENT'S T WITH N DEGREES OF FREEDOM IS LESS THAN T. REGT5620 C WRITTEN BY DR. F. GEBHARDT, DEUTSCHES RECHENZENTRUM. REGT5630 | C REGT5640 | APROX(Q )=1.0+ 1.0/(12.0*Q ) + 1.0/(288.0*Q **2) - .00REGT5650 | 12681/Q **3 REGT5660 | IF (N) 4,4,5 REGT5670 4 P=.5 REGT5680 | GO TO 80 REGT5690 | 5 AN=N REGT5700 | RN=1.0/AN REGT5710 | QN=SQRT(AN) REGT5720 8 IF(N-15) 10,30,30 REGT5730 | 10 K=100 REGT5740 | X=ABS(T)+10.0 REGT5750 | IF(N-6) 15,15,20 REGT5760 | 15 P= (QN/X)**N /QN REGT5770 14 GO TO 40 REGT5780 | 30 X=ABS(T)+6.0 REGT5790 | K=60 REGT5800 | 20 P=0.0 REGT5810 | 40 EN= -(AN+1.0)*.5 REGT5820 S = (1.0+X**2*RN)**EN*.5 REGT5830 19 | DO 50 I=1,K REGT5840 | S=S+ 2.0* (1.0+(X-.05)**2*RN)**EN REGT5850 23 | X=X-.1 REGT5860 | SN= (1.0+X**2*RN)**EN REGT5870 24 50 S=S+SN REGT5880 | P= (S- SN*.5)/30.0+P REGT5890 | IF(N-4) 55,65,65 REGT5900 | 55 IF(N-2) 60,61,62 REGT5910 | 60 PR=.31830989 REGT5920 GOTO 68 REGT5930 | 61 PR=.35355339 REGT5940 | GO TO 68 REGT5950 | 62 PR=.36755260 REGT5960 | GO TO 68 REGT5970 65 PR=.564190*((AN-1.0)/(AN-2.0))**((AN-1.0)*.5) * SQRT(RN*(AN-1.0) /REGT5980 | 15.4365637) * APROX((AN-1.0)*.5) / APROX(AN*.5-1.0) REGT5990 38 39 | 68 P=P*PR REGT6000 | IF(T) 80,70,75 REGT6010 | 70 P=.5 REGT6020 GO TO 80 REGT6030 | 75 P=1.0-P REGT6040 | 80 RETURN REGT6050 | END REGT6060 | 1 12/01/64 PAGE 25 REGT STORAGE MAP | 0 MAIN PROGRAM | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 24010 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE A 00000 R C 23420 R CON 23564 R | NUN 23730 I TX1 23767 R BEOB 23773 R | FQS 23774 R IUNE 23775 I IV 23776 I | LTAPE 23777 I MAX 24000 I MIN 24001 I | NBEOB 24002 I NHYP 24003 I NORT 24004 I NV 24005 I SSERR 24006 R TMIN 24007 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | KTAPE 24011 I NVAR 24012 I NFMT 24013 I | MTAPE 24014 I NAUS 24015 I IBIN 24016 I I 24017 I J 24020 I K 24021 I | KHYP 24022 I IEL 24023 I MAXP 24024 I | - ENTRY POINTS | 0 ...... SECTION 3 | - SUBROUTINES CALLED 0 DATV SECTION 4 .FRDD. SECTION 5 .FSLI. SECTION 6 | .FVIO. SECTION 7 .FRWT. SECTION 8 KOVA SECTION 9 | REGR SECTION 10 .FWRD. SECTION 11 .UN05. SECTION 12 | .FRTN. SECTION 13 .FCNV. SECTION 14 .UN06. SECTION 15 | .FFIL. SECTION 16 E.1 SECTION 17 E.2 SECTION 18 E.3 SECTION 19 E.4 SECTION 20 CC.1 SECTION 21 | CC.2 SECTION 22 CC.3 SECTION 23 CC.4 SECTION 24 | SYSLOC SECTION 25 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 1 FORMAT 24050 2 FORMAT 24054 100 4A 24071 | 102 18A 24154 103 25A 24177 110 46A 24324 | 200 101A 24557 170 70A 24423 120 61A 24400 | 130 65A 24412 180 89A 24510 175 85A 24465 | 173 81A 24461 190 97A 24544 THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 24602.| 1 12/01/64 PAGE 26 | DATV7 STORAGE MAP | 0 SUBROUTINE DATV | - ENTRY POINTS 0 DATV SECTION 2 | - SUBROUTINES CALLED | 0 SYSLOC SECTION 3 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00022.| 1 12/01/64 PAGE 27 | ELIMI STORAGE MAP | 0 SUBROUTINE ELIM - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | I 00001 I B 00002 R MINM 00003 I | MINP 00004 I | - ENTRY POINTS 0 ELIM SECTION 2 | - SUBROUTINES CALLED | 0 SYSLOC SECTION 3 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 10 5A 00023 30 34A 00137 20 26A 00106 | | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00235.| 1 12/01/64 PAGE 28 | 7KOVA STORAGE MAP 0 SUBROUTINE KOVA | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | FT 00001 R | - UNDIMENSIONED PROGRAM VARIABLES 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | JFT 00155 I JT 00156 I I 00157 I | L 00160 I K 00161 I AL 00162 R | MP 00163 I P 00164 R | - ENTRY POINTS 0 KOVA SECTION 2 | - SUBROUTINES CALLED | 0 .FRDD. SECTION 3 .FVIO. SECTION 4 .FRDB. SECTION 5 | .FWRB. SECTION 6 .UN05. SECTION 7 .FRTN. SECTION 8 | .FCNV. SECTION 9 .FRLR. SECTION 10 .FBLT. SECTION 11 .FBDT. SECTION 12 .FWLR. SECTION 13 E.1 SECTION 14 | E.2 SECTION 15 E.3 SECTION 16 E.4 SECTION 17 | CC.1 SECTION 18 CC.2 SECTION 19 CC.3 SECTION 20 | CC.4 SECTION 21 SYSLOC SECTION 22 | - EFN IFN CORRESPONDENCE 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 1 FORMAT 00203 20 20A 00264 30 48A 00450 | 40 60A 00500 70 91A 00633 50 69A 00532 | 60 80A 00604 85 141A 01036 80 130A 01016 | 90 153A 01123 100 162A 01137 THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 01366.| 1 12/01/64 PAGE 29 | 7MATIN STORAGE MAP | 0 SUBROUTINE MATIN | - UNDIMENSIONED PROGRAM VARIABLES 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | I 00001 I D 00002 R K 00003 I | L 00004 I P 00005 R T 00006 R | H 00007 R Y 00010 R | - ENTRY POINTS 0 MATIN SECTION 2 | - SUBROUTINES CALLED | 0 SYSLOC SECTION 3 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 1 7A 00034 7 51A 00227 2 19A 00103 | 3 27A 00127 4 37A 00176 5 44A 00213 | 6 48A 00222 8 62A 00271 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00430.| 1 12/01/64 PAGE 30 NORMT STORAGE MAP | 0 SUBROUTINE NRMT | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | ZW 00001 R ZR 00007 R II 00015 I Q 00023 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | S1 00027 R S3 00030 R S4 00031 R | NVAR 00032 I IJ 00033 I Y 00034 R IA 00035 I FN 00036 R SM1 00037 R | SM2 00040 R SM3 00041 R SM4 00042 R | ST3 00043 R ST4 00044 R U 00045 R | TX 00046 R | - ENTRY POINTS 0 NRMT SECTION 2 | - SUBROUTINES CALLED | 0 .FVIO. SECTION 3 .FRDB. SECTION 4 .FRWT. SECTION 5 | .FWRD. SECTION 6 SQRT SECTION 7 .XP3. SECTION 8 | .FRLR. SECTION 9 .FBLT. SECTION 10 .FBDT. SECTION 11 .UN06. SECTION 12 .FFIL. SECTION 13 .FCNV. SECTION 14 | E.1 SECTION 15 E.2 SECTION 16 E.3 SECTION 17 | E.4 SECTION 18 CC.1 SECTION 19 CC.2 SECTION 20 | CC.3 SECTION 21 CC.4 SECTION 22 SYSLOC SECTION 23 | - EFN IFN CORRESPONDENCE 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 40 32A 00400 12 13A 00327 20 21A 00344 | 15 17A 00332 25 23A 00346 30 31A 00363 | 28 27A 00351 2 FORMAT 00232 53 41A 00733 | 60 46A 01100 54 43A 01015 55 45A 01076 1 FORMAT 00144 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 01216.| 1 12/01/64 PAGE 31 | REDU STORAGE MAP | 0 SUBROUTINE REDUX - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | R 00001 R L 00002 I I 00003 I | - ENTRY POINTS | 0 REDUX SECTION 2 - SUBROUTINES CALLED | 0 SYSLOC SECTION 3 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 10 9A 00050 THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00143.| 1 12/01/64 PAGE 32 | REGRE STORAGE MAP | 0 SUBROUTINE REGR | - COMMON VARIABLES 0 COMMON BLOCK // ORIGIN 00001 LENGTH 24010 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | A 00000 R C 23420 R CON 23564 R | NUN 23730 I TX1 23767 R BEOB 23773 R | FQS 23774 R IUNE 23775 I IV 23776 I LTAPE 23777 I MAX 24000 I MIN 24001 I | NBEOB 24002 I NHYP 24003 I NORT 24004 I | NV 24005 I SSERR 24006 R TMIN 24007 R | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE BETA 24011 R D 24050 R KONV 25751 I | M 26115 I TX 26261 R Z 26271 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | ITYP 26330 I I 26331 I QHL 26332 R ITP 26333 I ITY 26334 I ITX 26335 I | INV 26336 I NULLP 26337 I KONS 26340 I | KABH 26341 I NABH 26342 I RUND 26343 R | IA 26344 I J 26345 I JA 26346 I | NA 26347 I SS 26350 R NFG 26351 I FG 26352 R SSERRM 26353 R SSW 26354 R | SMAX 26355 R CORY 26356 R ERRB 26357 R | T 26360 R P 26361 R XQ 26362 R | STXQ 26363 R RM 26364 R RPART 26365 R | CORR 26366 R SSM 26367 R SSH 26370 R SSHM 26371 R YQ 26372 R SSY 26373 R | SSYSQ 26374 R CMULTQ 26375 R CMULT 26376 R | F 26377 R SSERRU 26400 R SSERUM 26401 R | NFGU 26402 I IVU 26403 I SSU 26404 R | SSUM 26405 R SSD 26406 R NFGD 26407 I SSDM 26410 R FD 26411 R | - ENTRY POINTS | 0 REGR SECTION 3 | - SUBROUTINES CALLED | 0 ELIM SECTION 4 .FRDD. SECTION 5 .FWRD. SECTION 6 .FSLO. SECTION 7 MATIN SECTION 8 REDUX SECTION 9 | SQRT SECTION 10 TVRT SECTION 11 NRMT SECTION 12 | .UN05. SECTION 13 .FRTN. SECTION 14 .FCNV. SECTION 15 | 1 12/01/64 PAGE 33 | REGRE STORAGE MAP 0 .UN06. SECTION 16 .FFIL. SECTION 17 E.1 SECTION 18 | E.2 SECTION 19 E.3 SECTION 20 E.4 SECTION 21 | CC.1 SECTION 22 CC.2 SECTION 23 CC.3 SECTION 24 | CC.4 SECTION 25 SYSLOC SECTION 26 | - EFN IFN CORRESPONDENCE 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 2 FORMAT 26443 3 FORMAT 26444 4 FORMAT 26447 | 5 FORMAT 26452 6 FORMAT 26550 7 FORMAT 26556 | 8 FORMAT 26566 9 FORMAT 26610 10 FORMAT 26621 | 11 FORMAT 26626 12 FORMAT 26651 13 FORMAT 26671 14 FORMAT 26742 15 FORMAT 26750 16 FORMAT 27002 | 17 FORMAT 27052 18 FORMAT 27173 19 FORMAT 27177 | 20 FORMAT 27205 21 FORMAT 27227 22 FORMAT 27231 | 23 FORMAT 27261 25 FORMAT 27313 26 FORMAT 27332 | 30 6A 27362 60 24A 27427 50 17A 27405 40 12A 27373 110 61A 27526 123 91A 27613 | 122 76A 27561 90 48A 27475 70 39A 27466 | 80 44A 27471 125 114A 27667 124 111A 27654 | 135 136A 27740 132 117A 27673 137 140A 27751 | 136 138A 27742 145 143A 27764 140 141A 27761 150 146A 27772 157 166A 30061 151 149A 30005 | 155 158A 30031 160 168A 30073 1605 170A 30075 | 1615 182A 30151 161 175A 30124 162 187A 30157 | 163 193A 30173 167 210A 30236 165 202A 30217 | 169 214A 30242 170 223A 30271 175 233A 30310 182 251A 30370 180 244A 30350 184 253A 30374 | 192 315A 31066 185 257A 30400 190 312A 31055 | 186 279A 30545 189 310A 31033 1865 288A 30646 | 187 306A 31021 1900 320A 31131 1902 330A 31154 | 1901 325A 31137 1903 333A 31205 1904 335A 31274 200 337A 31300 230 346A 31407 250 354A 31573 | | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 31675.| 1 12/01/64 PAGE 34 | 7TVRT STORAGE MAP 0 SUBROUTINE TVRT | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | AN 00001 R RN 00002 R QN 00003 R | K 00004 I X 00005 R EN 00006 R S 00007 R SN 00010 R PR 00011 R | - ENTRY POINTS | 0 TVRT SECTION 2 | - SUBROUTINES CALLED | 0 SQRT SECTION 3 .XP2. SECTION 4 .XP3. SECTION 5 E.1 SECTION 6 E.2 SECTION 7 E.3 SECTION 8 | E.4 SECTION 9 CC.1 SECTION 10 CC.2 SECTION 11 | CC.3 SECTION 12 CC.4 SECTION 13 SYSLOC SECTION 14 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 4 5A 00115 5 7A 00120 80 45A 00421 | 10 10A 00144 30 16A 00171 15 13A 00156 | 20 17A 00177 40 18A 00200 50 25A 00271 | 55 29A 00312 65 37A 00327 60 31A 00316 | 61 33A 00321 62 35A 00324 68 40A 00404 70 42A 00413 75 44A 00416 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00467. 1 12/01/64 PAGE 35 | | 0$DATA 1 IBLDR 12/01/64 PAGE 36 | * MEMORY MAP * | 0 | 0SYSTEM 00000 THRU 02717| FILE BLOCK ORIGIN 02720 FILES 1. UNIT01 | 2. UNIT02 | 3. UNIT03 | 4. UNIT04 | 5. UNIT05 6. UNIT06 | 7. UNIT07 | 8. UNIT08 | FILE LIST ORIGIN 03060| PRE-EXECUTION INITIALIZATION 03100 CALL ON OBJECT PROGRAM 03137| OBJECT PROGRAM 03144 THRU 32607| 0 DECK ORIGIN CONTROL SECTIONS (/NAME/=NON 0 LENGTH, (LOC)=DELETED, 'LOC'=MOVED, *=NOT REFERENCED)| | 1. REGT 03144 /// /(53770) ...... 03720 * 2. DATV7 03736 DATV 03746 | 3. ELIMI 03760 ELIM 04135 | 4. 7KOVA 04215 KOVA 05425 | 5. 7MATIN 05603 MATIN 06155 | 6. NORMT 06233 NRMT 07352 7. REDU 07451 REDUX 07570 | 8. REGRE 07614 /// /(53770) REGR 15454 | 9. 7TVRT 15501 TVRT 16125 | 10. .LXCON 16170 .LXSTR 16170 * .LXSTP 16173 .LXOUT 16241 * .LXERR 16250 .LXCAL 16253 * | .LXRTN 16253 IBEXIT 16253 * .DBCLS 16435 * .LXARG 16604 * .L0 16627 * .CLSE 16635 .LFBL 16636 * .LUNB 16637 .DFOUT 16640 | 11. .IODEF 16644 .DEFIN 16644 .ATTAC 16650 * .CLOSE 16652 .OPEN 16654 .READ 16656 | .WRITE 16660 .BSR 16670 * .READR 16700 .RELES 16702 * .LAREA 16713 | .LFBLK 16731 .LTSX 16734 * .RLHLD 16744 * .AREA1 16747 .LUNBL 16755 | .ENTRY 16761 .GOA 17016 .GO 17022 .DERR 17036 .NOPXI 17037 .COMXI 17041 .EX34 17063 | 12. .LXSL 17070 .LXSLB 17070 * .LXSEL 17071 .LXSL1 17072 .LXTST 17075 * .LXOVL 17140 * | .LXMOD 17202 * .LXIND 17226 * .LXDIS 17231 * .LXFLG 17232 * .LTCH 17233 | 13. .FPTRP 17242 .FFPT. 17242 * .FPOUT 17371 * .FPARG 17377 * /.COUNT/ 17401 * OVFLOW 17445 * | 14. .ERAS. 17452 E.1 17452 E.2 17453 E.3 17454 E.4 17455 15. .XCC. 17456 CC.1 17456 CC.2 17457 CC.3 17460 CC.4 17461 | 16. FXEM 17462 .FXEM. 17462 TRACE 17470 * /.OPTW./ 20025 * | 17. FCOM 20041 .FCOM. 20041 /FORCOM/ 20144 * | 18. FRDB 20156 .FRDB. 20156 .FWRB. 20160 | 19. FRDD 20210 .FRDD. 20210 .FWRD. 20213 20. FRWT 20325 .FRWT. 20325 | 21. FIOS 20413 .FIOS. 20413 .FSEL. 20602 .FRTB. 20617 .FRTD. 20634 /.FBF1./ 21066 * | 22. FIOB 21072 .FIOB. 21072 .FBLT. 21134 .FBDT. 21162 .FRLR. 21233 .FWLR. 21233 | 23. FIOH 21527 .FIOH. 21527 .FCNV. 22306 .FFIL. 22403 .FRTN. 22403 | 24. FCNV 22522 .FCON. 22522 25. FVIO 25610 .FVIO. 25610 ..FVIO 25610 * | 26. FIOT 25613 .FIOT. 25613 .FIOE. 25704 | 27. UN01 25732 .UN01. 25732 | 28. UN02 25733 .UN02. 25733 | 29. UN03 25734 .UN03. 25734 30. UN04 25735 .UN04. 25735 | 1 IBLDR 12/01/64 PAGE 37 | * MEMORY MAP * | 0 31. UN05 25736 .UN05. 25736 | 32. UN06 25737 .UN06. 25737 ..UN06 25737 * 33. UN07 25740 .UN07. 25740 | 34. UN08 25741 .UN08. 25741 | 35. FLOG 25742 ALOG10 25742 * ALOG 25743 | 36. FXPF 26146 EXP 26146 | 37. FSQR 26267 SQRT 26267 38. FXP2 26342 .XP2. 26342 | 39. FXP3 26460 .XP3. 26460 | 40. FSLDI 26605 .FSLI. 26623 .FSDI. 26631 * | 41. FSLI 26642 .SLI. 26642 .SLI1. 26647 .SDI. 26655 .SDI1. 26663 | 42. FSLDO 26701 .FSLO. 26717 .FSDO. 26725 * 43. FSLO 26736 .SLO. 26736 .SLO2. 26744 .SDO. 26751 .SDO2. 26763 | 44. .IOCS 26775 .L(0) 26775 .MONSW 27015 .TEOR 27064 .DEFI. 27144 .JOINX 27210 * | .CLOS. 27227 .ATTC. 27242 .SH1 27454 * .SH9 27516 * .OPEN. 27537 | .OP4 27565 * .OP7 27616 * .OP9.2 27632 * .RLSE. 27704 .RER2. 27704 | .READ. 27705 .RER1. 27730 .WRIT. 27732 .MNT1A 30122 * .EOFEX 30203 * .FEEIT 30253 .GTIOX 30274 .RW7 30412 * .RE7 31031 * .ENDTR 31472 | .SEL59 31474 * .BSR. 32113 .EOTOF 32240 .ETOF3 32246 * .SWITC 32275 | .TCHEX 32602 .BASIO 32605 * | 45. .IOCSM 32610 | 46. // 53770 0I/O BUFFERS 32610 THRU 53743| 0UNUSED CORE 53744 THRU 53767| 0 1 REGRESSION ANALYSIS DEUTSCHES RECHENZENTRUM TEST DATA FOR ZF REGT DARMSTADT STEPWISE REGRESSION UNRESTRICTED HYPOTHESIS DEPENDENT VARIABLE 2 INDEPENDENT VARIABLES 4 7 9 1 5 8 0 0VAR.NO. MEAN STAND.DEV. REGRESS. STAND.DEV.OF T D.F. PART. ERROR PROB. MULT. CORR. OF MEAN COEFF. REG.COEFF. CORR. (2-SIDED) CORR. WITH Y 4 1.52000 0.32069 0.37281 0.14251 2.616 3 0.8338 0.0793 0.9136 0.2418 7 3.23000 0.37596 0.12724 0.10028 1.269 3 0.5910 0.2940 0.8700 0.8349 9 8.42000 2.01168 0.11717 0.03555 3.296 3 0.8852 0.0459 0.9656 0.6146 1 0.55000 0.09574 -1.79129 0.72199 -2.481 3 -0.8200 0.0892 0.9632 0.0490 5 5.42000 0.44692 -0.19694 0.11297 -1.743 3 -0.7094 0.1797 0.9298 -0.7978 8 4.74000 2.17982 -0.05499 0.01785 -3.081 3 -0.8717 0.0541 0.8785 0.3578 0 1.72900 0.68267 2.533 3 0.0852 0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES. 0 0NUMBER OF OBSERVATIONS 10 TOTAL SUM OF SQUARES 2.41000E 01 10 D.F. ERROR SUM OF SQUARES 9.32624E-02 3 D.F. MEAN SQUARE ERROR 3.10875E-02 MEAN VALUE OF THE DEPENDENT VARIABLE 1.38000E 00 STANDARD-DEV. OF THE MEAN VALUE 2.37018E-01 MULTIPLE CORRELATION 0.99073 SQUARE OF THE MULTIPLE CORRELATION 0.98155 F( 6, 3) = 26.606 0 0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM 9.31791E-02 0 0TEST OF RESIDUES ON NORMAL DISTRIBUTION TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = 0.383 TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = -0.564 0DEVIATION OF RESIDUES FROM NORMAL DISTRIBUTION IS NOT SIGNIFICANT 1 TEST DATA FOR ZF REGT 1. RESTRICTED HYPOTHESIS STEPWISE REGRESSION DEPENDENT VARIABLE 2 INDEPENDENT VARIABLES 4 9 1 5 8 0STEPWISE REGRESSION. VARIABLE NO. 7 HAS BEEN ELIMINATED. 0 0VAR.NO. MEAN STAND.DEV. REGRESS. STAND.DEV.OF T D.F. PART. ERROR PROB. MULT. CORR. OF MEAN COEFF. REG.COEFF. CORR. (2-SIDED) CORR. WITH Y 4 1.52000 0.32069 0.50009 0.10865 4.603 4 0.9172 0.0100 0.8199 0.2418 9 8.42000 2.01168 0.12835 0.03697 3.471 4 0.8665 0.0256 0.9634 0.6146 1 0.55000 0.09574 -1.91204 0.76820 -2.489 4 -0.7795 0.0676 0.9625 0.0490 5 5.42000 0.44692 -0.26433 0.10702 -2.470 4 -0.7772 0.0690 0.9089 -0.7978 8 4.74000 2.17982 -0.06723 0.01612 -4.172 4 -0.9017 0.0140 0.8232 0.3578 0 2.34215 0.51758 4.525 4 0.0106 0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES. 0 0 MEAN VALUE OF THE DEPENDENT VARIABLE 1.38000E 00 STANDARD-DEV. OF THE MEAN VALUE 2.37018E-01 MULTIPLE CORRELATION 0.98573 SQUARE OF THE MULTIPLE CORRELATION 0.97166 F( 5, 4) = 27.433 0 0 ANALYSIS OF VARIANCE FOR TESTING THE RESTRICTED HYPOTHESIS SOURCE SUM OF SQ. D.F. MEAN SQ. F TOTAL 2.41000E 01 10 2.41000E 00 UNRESTR. HYPOTHESIS 2.40067E 01 7 3.42953E 00 RESTR. HYPOTHESIS 2.39567E 01 6 3.99279E 00 DIFFERENCE 5.00009E-02 1 5.00009E-02 1.608 ERROR 9.32624E-02 3 3.10875E-02 0 0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM 1.43228E-01 0 0TEST OF RESIDUES ON NORMAL DISTRIBUTION TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = 0.363 TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = -0.716 0DEVIATION OF RESIDUES FROM NORMAL DISTRIBUTION IS NOT SIGNIFICANT 1 TEST DATA FOR ZF REGT 2. RESTRICTED HYPOTHESIS STEPWISE REGRESSION DEPENDENT VARIABLE 2 INDEPENDENT VARIABLES 4 9 1 8 0STEPWISE REGRESSION. VARIABLE NO. 5 HAS BEEN ELIMINATED. 0 0VAR.NO. MEAN STAND.DEV. REGRESS. STAND.DEV.OF T D.F. PART. ERROR PROB. MULT. CORR. OF MEAN COEFF. REG.COEFF. CORR. (2-SIDED) CORR. WITH Y 4 1.52000 0.32069 0.55253 0.15145 3.648 5 0.8526 0.0148 0.8119 0.2418 9 8.42000 2.01168 0.20462 0.02890 7.081 5 0.9536 0.0009 0.8730 0.6146 1 0.55000 0.09574 -3.50632 0.59200 -5.923 5 -0.9355 0.0020 0.8659 0.0490 8 4.74000 2.17982 -0.07607 0.02234 -3.406 5 -0.8359 0.0191 0.8129 0.3578 0 1.10629 0.18815 5.880 5 0.0020 0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES. 0 0 MEAN VALUE OF THE DEPENDENT VARIABLE 1.38000E 00 STANDARD-DEV. OF THE MEAN VALUE 2.37018E-01 MULTIPLE CORRELATION 0.96356 SQUARE OF THE MULTIPLE CORRELATION 0.92845 F( 4, 5) = 16.220 0 0 ANALYSIS OF VARIANCE FOR TESTING THE RESTRICTED HYPOTHESIS SOURCE SUM OF SQ. D.F. MEAN SQ. F TOTAL 2.41000E 01 10 2.41000E 00 UNRESTR. HYPOTHESIS 2.40067E 01 7 3.42953E 00 RESTR. HYPOTHESIS 2.37382E 01 5 4.74765E 00 DIFFERENCE 2.68500E-01 2 1.34250E-01 4.318 ERROR 9.32624E-02 3 3.10875E-02 0 0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM 3.61721E-01 0 0TEST OF RESIDUES ON NORMAL DISTRIBUTION TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = 0.525 TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = -1.049 0DEVIATION OF RESIDUES FROM NORMAL DISTRIBUTION IS NOT SIGNIFICANT 1 REGRESSION ANALYSIS DEUTSCHES RECHENZENTRUM TEST DATA FOR ZF REGT DARMSTADT SECOND GROUP OF TESTS UNRESTRICTED HYPOTHESIS DEPENDENT VARIABLE 1 INDEPENDENT VARIABLES 5 4 7 8 9 0 0VAR.NO. MEAN STAND.DEV. REGRESS. STAND.DEV.OF T D.F. PART. ERROR PROB. MULT. CORR. OF MEAN COEFF. REG.COEFF. CORR. (2-SIDED) CORR. WITH Y 5 5.42000 0.44692 0.10533 0.05785 1.821 4 0.6732 0.1428 0.8673 0.3532 4 1.52000 0.32069 0.09351 0.08691 1.076 4 0.4737 0.3427 0.8870 0.7889 7 3.23000 0.37596 -0.01831 0.06884 -0.266 4 -0.1318 0.8037 0.8676 0.2701 8 4.74000 2.17982 -0.00915 0.01148 -0.797 4 -0.3701 0.4703 0.8577 0.4770 9 8.42000 2.01168 0.04277 0.01220 3.506 4 0.8686 0.0248 0.8514 0.6963 0 -0.42065 0.42341 -0.993 4 0.3768 0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES. 0 0NUMBER OF OBSERVATIONS 10 TOTAL SUM OF SQUARES 3.85000E 00 10 D.F. ERROR SUM OF SQUARES 5.96402E-02 4 D.F. MEAN SQUARE ERROR 1.49101E-02 MEAN VALUE OF THE DEPENDENT VARIABLE 5.50000E-01 STANDARD-DEV. OF THE MEAN VALUE 9.57427E-02 MULTIPLE CORRELATION 0.96318 SQUARE OF THE MULTIPLE CORRELATION 0.92771 F( 5, 4) = 10.266 0 0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM 5.96372E-02 0 0TEST OF RESIDUES ON NORMAL DISTRIBUTION TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = -1.157 TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = -0.262 0DEVIATION OF RESIDUES FROM NORMAL DISTRIBUTION IS NOT SIGNIFICANT 1 TEST DATA FOR ZF REGT 1. RESTRICTED HYPOTHESIS SECOND GROUP OF TESTS DEPENDENT VARIABLE 1 FIRST RESTRICTED HYPOT H. INDEPENDENT VARIABLES 7 8 0COEFFICIENTS SPECIFIED IN HYPOTHESIS NO VALUE NO VALUE NO VALUE NO VALUE NO VALUE NO VALUE 4 1.50000E-01 0REGRESSION THROUGH THE ORIGIN 0 0 INVERSE OF CROSS PRODUCT MATRIX 7 1.440318E-02 -3.891588E-03 8 -3.891588E-03 2.584456E-03 0 0VAR.NO. MEAN STAND.DEV. REGRESS. STAND.DEV.OF T D.F. PART. ERROR PROB. MULT. CORR. OF MEAN COEFF. REG.COEFF. CORR. (2-SIDED) CORR. WITH Y 7 3.23000 0.37596 0.07480 0.02822 2.650 8 0.6838 0.0292 0. 0.2701 8 4.74000 2.17982 0.00974 0.01196 0.815 8 0.2768 0.4388 0. 0.4770 0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES. 0 0 MEAN VALUE OF THE DEPENDENT VARIABLE 3.22000E-01 STANDARD-DEV. OF THE MEAN VALUE 6.49153E-02 0 0 ANALYSIS OF VARIANCE FOR TESTING THE RESTRICTED HYPOTHESIS SOURCE SUM OF SQ. D.F. MEAN SQ. F TOTAL 1.41610E 00 10 1.41610E-01 UNRESTR. HYPOTHESIS 1.35646E 00 6 2.26077E-01 RESTR. HYPOTHESIS 9.73688E-01 2 4.86844E-01 DIFFERENCE 3.82772E-01 4 9.56929E-02 6.418 ERROR 5.96402E-02 4 1.49101E-02 THESE RESULTS REFER TO THE LINEAR COMBINATION OF THE DEPENDENT VARIABLE WITH THOSE VARIABLES WHOSE COEFFICIENTS ARE SPECIFIED BY THE HYPOTHESIS. 0 0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM 4.42412E-01 0 0TEST OF RESIDUES ON NORMAL DISTRIBUTION TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = -0.491 TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = -0.428 0DEVIATION OF RESIDUES FROM NORMAL DISTRIBUTION IS NOT SIGNIFICANT 1 TEST DATA FOR ZF REGT 2. RESTRICTED HYPOTHESIS SECOND GROUP OF TESTS DEPENDENT VARIABLE 1 SECOND RESTR. HYPOTHES IS INDEPENDENT VARIABLES 5 9 0 0VAR.NO. MEAN STAND.DEV. REGRESS. STAND.DEV.OF T D.F. PART. ERROR PROB. MULT. CORR. OF MEAN COEFF. REG.COEFF. CORR. (2-SIDED) CORR. WITH Y 5 5.42000 0.44692 0.14569 0.02857 5.099 7 0.8876 0.0014 0.3499 0.3532 9 8.42000 2.01168 0.04447 0.00635 7.005 7 0.9355 0.0002 0.3499 0.6963 0 -0.61407 0.18418 -3.334 7 0.0125 0'MULT. CORR.' IS THE MULTIPLE CORRELATION BETWEEN THE RESPECTIVE VARIABLE AND ALL OTHER INDEPENDENT VARIABLES. 0 0 MEAN VALUE OF THE DEPENDENT VARIABLE 5.50000E-01 STANDARD-DEV. OF THE MEAN VALUE 9.57427E-02 MULTIPLE CORRELATION 0.94378 SQUARE OF THE MULTIPLE CORRELATION 0.89072 F( 2, 7) = 28.529 0 0 ANALYSIS OF VARIANCE FOR TESTING THE RESTRICTED HYPOTHESIS SOURCE SUM OF SQ. D.F. MEAN SQ. F TOTAL 3.85000E 00 10 3.85000E-01 UNRESTR. HYPOTHESIS 3.79036E 00 6 6.31727E-01 RESTR. HYPOTHESIS 3.75985E 00 3 1.25328E 00 DIFFERENCE 3.05117E-02 3 1.01706E-02 0.682 ERROR 5.96402E-02 4 1.49101E-02 0 0SUM OF SQUARES ERROR + DIFFERENCE, CHECKSUM 9.01515E-02 0 0TEST OF RESIDUES ON NORMAL DISTRIBUTION TEST STATISTIC DERIVED FROM 4TH SAMPLE MOMENT = -0.049 TEST STATISTIC DERIVED FROM 3RD SAMPLE MOMENT = -0.083 0DEVIATION OF RESIDUES FROM NORMAL DISTRIBUTION IS NOT SIGNIFICANT 1 1278 LINES OUTPUT. $IBSYS $STOP | 1 PERIPHERAL UNIT POSITIONS AT END OF JOBS | 1 SYSPP1 IS B2 REC. 00296, FILE 00000 SYSOU1 IS B1 REC. 00466, FILE 00000 SYSIN1 IS A2 REC. 00002, FILE 00001| 1 END OF JOBS | 1