Final Report on the Status of Minorities and Women in the Department of Labor

427.

428.

429.

430.

129 FORMAT('1PLOT OF STANDARDIZED RESIDUALS'///15x,'-5S',7X,'-45',7x,
1-3 S',7X,'-2S',
25S/16X,5('I.

DO 1580 I=1,M

7X,'-1S',8x,'0',9×‚'1S',8x,'2S',8X,'35',8X,'45',8X, I'))

458.

459.

460.

461.

462.

463.

464.

131 FORMAT (///, (' PERCENTAGE OF DEVIATIONS WITHIN',12,' STANDARD ERROR 1S OF ZERO IS',F10.4))

GO TO 590

1590 WRITE(K0,1591)

1591 FORMAT('1PLOT OF STANDARDIZED RESIDUALS' /5X, ALL RESIDUALS ARE ZER 10--NO PLOTTING WILL BE DONE')

590 CONTINUE

515.

516.

517.

518.

519.

520.

521.

522.

523.

524.

525.

2121 KGPH (J)=0

DO 620KKK=1,NPLTS

READ PLOT CARD

READ(K5,111)

GET REQUISITE DATA FROM K8

FIND(K8 KCOL)

READ(K8 KCOL)(X(I,LCOL),I=1,M)

DO 625 KK=1,4

TAMP=0.

DO 626 I=1,M

626 TAMP=TAMP+X (I,LCOL)

TAMP TAMP/M

DO 625 I=1, M

625 XII,LCOL)=X(I,LCOL)-TAMP

STDE=VPROD(LCOL,LCOL,1.000)
STDE=DSQRT(STDE/(M-1.)}

WRITE(KO,631) NAME (KCOL)

631 FORMAT(1PLOT OF STANDARDIZED RESIDUALS (Y-AXIS)' ,1X, VS STANDARDI
IZED,A8, (X-AXIS) '//12X,'-5',8×,'-4', 8X,'-3', 8X,
•
'-2',8x,'-1',9×,
2'0',9x,'1',9×‚'2',9×‚'3',9x,'4',9×, '5'/13x,101('

KAMP=6

TEMP=26.*TUMF

IROW=1

DO 2120 K=1,51

ISTRT=IROW

TEMP=TEMP-TUMP

DO 2121 J=1,101

IF(IROW.GT.M) GO TO 2125

DO 2130 I=I STRT,M

IF(X(ICOL(I), NCOL).LT.TEMP) GO TO 2125

IPLT=10.*X(ICOL(I),LCOL)/STDE+51.5

IF(IPLT.GT.101) IPLT=101

IF(IPLT.LT.1) IPLT=1

KGPH (IPLT)=KGPH(IPLT)+1

IROW=IROW+1

2130 CONTINUE

2125 DO 2126 J=1,101

IF(KGPH(J).GT.5) KGPH(J)=5

NCODE=KGPH(J) +1

GO TO (2140,2141,2142,2143,2144,2145),NCODE

2140 KGPH(J) =KBL NK

GO TO 2126

2141 KGPH(J)=NUMBR (1)

GO TO 2126

545.

2143 KGPH(J)=NUMBR (3)

702 FORMAT('1THESE ARE THE DATA USED TO MAKE PROJECTIONS //

703 FORMAT(/'OPROJECTED VALUE OF, A8, FOR PROJECTION PERIOD ',14,' I 1S,G17.6//3X, 'THE STANDARD ERROR OF FORECAST FOR THIS PROJECTION I 2S,G17.6)

IF(NPROJ.EQ.C) GO TO 7CO

READ AND PRINT PROJECTION DATA

WRITE(KO,702)

DO 710 J=1,N

JJ=JCOL (J)

READ(K5,701) (X(I, JJ),I=1,NPROJ)

710 WRITE(KO, 1051)JJ, NAME (JJ), (X(I,JJ),I=1, NPROJ)

1051 FORMAT('OVARIABLE', 14, '--', A 8/ (1X, 7616.7))

COMPUTE X-TRANSPOSE-X INVERSE

DO 720 I=1, N

IF(I.EQ.N) GO TO 720

II = I + 1

DO 721 J=II, N

TEMP=0.

DO 725 K=J,N

725 TEMP=TEMP+X (M1+JCOL (I),JCOL(K))*X(M1+JCOL (J),JCOL(K))/D(JCOL(K)} X (M1+JCOL(J), JCOL (I))=TEMP

721 CONTINUE

720 CONTINUE

MAKE PROJECTIONS

DO 730 K=1, NPROJ

PROJ=0.

IF(CSPQ.NE. YES) PROJ=-X(M1,NCOL)

DO 731 J=1,N

731 PROJ=PROJ-X (K,JCOL (J))*X(M1+JCOL (J), NCOL)
COMPUTE STANDARD ERROR OF FORECAST
SEF=0.

IF(CSPQ.EQ.YES) GO TO 735

TAKE OUT MEANS IF CONSTANT IN EQUATION
DO 732 J=1,N

732 X (K, JCOL{J}}=X(K, JCOL (J))-XBAR (JCOL (J))
735 DO 736 I=1,N

605.

II=JCOL (I)

628.

629.

630.

631.

632.

633.

634.

635.

636.

637.

638.

639.

640.

641.

642.

643.

644.

645.

646.

647.

648.

649.

650.

651.

652.

653.

654.

655.

656.

REAL*8 NAME

657.

658.

659.

660.

661.

662.

663.

664.

665.

666.

667.

668.

669.

670.

READ(K5,103) ENDQ

IF(ENDQ.EQ.XMOR) GO TO 5000

STOP

END

SUBROUTINE ELIM(I, J)

IMPLICIT INTEGER*2 (I-N), REAL*8(A-H,0-2)
INTEGER 4 KO

REAL 8 NAME

DIMENSION X(820,10),XINV(11),SS(10),D(10),JCOL(10), NAME (10)
COMMON NAME,X,XINV,SS,D,B,STDE,FS, PCC,TV,BETA,ELAS,CSPQ, TOL
COMMON KO,JCOL,M,M1, MN, NN,NSTOP,NT,NIN, NEND,N,NREQ, NF
TEMP=VPROD(I, J,D(I))

00 10 K=1,MN

10 X (K,J) = X(K,J)-TEMP*X(K,I)

D(J)=VPROD(J,J,1.0)

RETURN

END

REAL FUNCTION VPROD*8( I, J, DUMMY)

IMPLICIT INTEGER*2(I-N), REAL 8(A-H,O-Z)

INTEGER 4 KO

REAL*8 NAME

DIMENSION X (820,10),XINV(11),SS(10),D(10), JCOL (10),NAME(10)
COMMON NAME,X,XINV,SS,D,B,ST DE,FS, PCC,TV, BETA, ELAS, CSPQ, TOL
COMMON KO,JCOL,M,M1,MN,NN,NSTOP,NT,NIN,NEND,N,NREQ, NF
TEMP=0.000

TUMP DUMMY

DO 10 K=1,M
TAMP=X(K,I)

TOMP=X(K, J)

10 TEMP=TEMP+T AMP/TUMP*TOMP

VPROD=TEMP

RETURN

END

SUBROUTINE WRT(NCODE, I)

IMPLICIT INTEGER *2(I-N), REAL*8(A-H,O-Z)

INTEGER 4 KO,NCODE

DIMENSION X (820,10),XINV(11),SS(10),D(10), JCOL(10), NAME (10)
COMMON NAME,X,XINV,SS,D,B,ST DE,FS,PCC,TV,BETA, ELAS,CSPQ, TOL
COMMON KO,JCOL,M,M1, MN, NN,NSTOP,NT,NIN, NEND,N,NREQ, NF
DATA CONS/ CONSTANT'/

GO TO(10,20,30,40,50,60,70,80,90,100,110,120), NCODE

10 WRITE(KO, 1) NAME(JCOL(I))

GO TO 1100

20 WRITE(KO,2) I,NAME(I), (X(J, I), J=1,M)

GO TO 1100

30 WRITE(K0,3) I,(X( J, I), J= 1,M)

GO TO 1100

40 WRITE(KO,4) I,(X( M1 + I, J), J=1,NF)

GO TO 1100

50 WRITE(K0,5)

671.

672. 673.

674. 675.

676.

677.

678.

679.

680.

681.

682.

683.

684.

685.

686.

687.

688.

689.

690. 691. 692.

693.

694.

695.

696.

697.

GO TO 1100

60 WRITE(K0,6)

GO TO 1100

70 WRITE(KO,7) CONS, B, STDE,FS

GO TO 1100

80 WRITE(KO,7) NAME(JCOL(I)),B,STDE,FS,TV,PCC
GO TO 1100

90 WRITE(KO,7) NAME(JCOL(I)), BETA, ELAS
GO TO 1100

100 WRITE(KO,7) NAME (JCOL (I)),B,STDE,FS
GO TO 1100

110 WRITE(K0,8)

GO TO 1100

120 WRITE(KO,9) NAME(JCOL(I)),TV,FS, PCC 1100 RETURN

1 FORMAT(15X,A8)

2 FORMAT('OVARIABLE',13,'--', A 8/ (1X, 7616.7))

3 FORMAT('OVARIABLE, 13/(1X,7616.7))

4 FORMAT('OROW',14/(5X, 7616.7))

5 FORMAT('ONOTE--CONSTANT SUPPRESS OPTION HAS BEEN CHOSEN/10X, ALL 1STATISTICS GIVEN ARE NOT CENTERED ABOUT THE MEANS'///Y

8 FORMAT(2X, VARIABLE',5X,'REGRESSION COEF',6X,'STANDARD ERROR, 5X, 1F-STATISTIC '/)

9 FORMAT (2X,A8,5X,F9.6,5X,G11.4,5X,F9.6)

END

//GO.FT08F001 DD DSNAME=&JUNK, DISP=(NEW, DELETE),UNIT=SYSDA,
// SPACE=(TRK, (30, 30))

//GO.FT05F001 DD *

TITLE SUPRB TEST--EEO TASK FORCE PROJECT

FORMAT(4F1.0,2F6.0,F2.0,1X,2F6.0,1X,9F1.0,1X,8F1.0)

OPTIONS DPNT, TPNT, MNSD, CORR, PLOT, BETA

Books