SUBROUTINE MDL2(N,IYEAR,ID,K,KS,THCHI,HEADER) IMPLICIT DOUBLE PRECISION (A-G,O-Z) INTEGER N(20), IYEAR(21), ID(21,21) DOUBLE PRECISION THCHI(100) CHARACTER HEADER*80, FORMT*30 DOUBLE PRECISION LOWF,HIF,LOWAVF,HIAVF,LOWS,HIS DOUBLE PRECISION MLS,LOWMLS,HIMLS INTEGER SUMDF COMMON/SNEAK/TEST,TEST2,TEST3,IDF,IDF2,IDF3 COMMON/PATCH/TEMP1,TEMP2 ,TEMP3 COMMON /MDL2CM/ THETA(20),RHO(20), 1TPRIM(20),TDPRIM(20),RPRIM(20),RDPRIM(20),XI(20),DELTA(21), 2DIFF(21),SPROD(20,20),T(20),F(21), 3S(20),DR(20),RB(20),RSAVE(20),SP(20), 4V(21,21),CORR(21,21),LV(21),MV(21), 5SDEVF(21),LOWF(21),HIF(21),AHAT(20),EE(20,20),CHISQS(20,20), 6ISAVE(20) INTEGER H DOUBLE PRECISION R(20), C(20) INTEGER IR(20), IC(20) DATA IR,IC/40*0/,R,C,FHAT,SHAT/42*0.0D0/ TEST=0.0 TEST2=0.0 TEST3=0.0 FHAT=0.0 SHAT=0.0 CHISUM=0.0 IDF=0 IDF2=0 IDF3=0 DO 5 I=1,20 R(I)=0.0 C(I)=0.0 IR(I)=0 C(I)=0 THETA(I)=0.0 RHO(I)=0.0 TPRIM(I)=0.0 TDPRIM(I)=0.0 RPRIM(I)=0.0 RDPRIM(I)=0.0 XI(I)=0.0 DELTA(I)=0.0 DIFF(I)=0.0 T(I)=0.0 F(I)=0.0 S(I)=0.0 DR(I)=0.0 RB(I)=0.0 RSAVE(I)=0.0 DO 5 J=1,20 SPROD(I,J)=0.0 V(I,J)=0.0 CORR(I,J)=0.0 5 CONTINUE C CALL MOD2W1 WRITE (6,55) HEADER 55 FORMAT (' ',20X,A///) WRITE (6,60) 60 FORMAT (' ','YEAR NUMBER',35X,'RECOVERY MATRIX') WRITE (6,65) 65 FORMAT (' ',5X,'BANDED') WRITE (6,70) 70 FORMAT ('+','____ ______',35X,'_______________'/) DO 80 I=1,K ISS=(I-1)*6+1 WRITE(FORMT,75) ISS WRITE (6,FORMT) IYEAR(I),N(I),(ID(I,JJ),JJ=I,KS) 75 FORMAT ('(1X,I4,1X,I5,',I2,'X,20I6)') 80 CONTINUE C L=K-1 LL=K-2 KK=K+1 KJ=K+1 NR=0 C DO 85 I=1,K DO 85 J=1,KS IR(I)=IR(I)+ID(I,J) R(I)=R(I)+ID(I,J) 85 CONTINUE DO 90 I=1,KS DO 90 J=1,K C(I)=C(I)+ID(J,I) IC(I)=IC(I)+ID(J,I) 90 CONTINUE C C IS=KS-K T(1)=R(1) DONE=0.0D0 DTHREE=0.0D0 DO 95 I=1,K DONE=DONE+R(I) DTHREE=DTHREE+T(I) F(I)=R(I)*C(I)/(N(I)*T(I)) FHAT=FHAT+F(I) IF (I.EQ.K) GO TO 95 S(I)=R(I)*(1.0D0-C(I)/T(I))*N(I+1)/(N(I)*R(I+1)) SHAT=SHAT+S(I) T(I+1)=R(I+1)+T(I)-C(I) 95 CONTINUE FHAT=FHAT/K SHAT=SHAT/(K-1) DIFF(1)=1.0D0 DIFF(2)=1.0D0 IF (IS.LE.1) TPRIM(KS-1)=1.0D0 IF (IS.LE.2) TDPRIM(KS-2)=2.0D0 THETA(KS)=1.0D0 D1=0.0D0 SUM=DTHREE-DONE IF (KS.EQ.K) GO TO 105 DO 100 J=KJ,KS T(J)=T(J-1)-C(J-1) F(J)=FHAT 100 D1=D1+T(J) SUM=SUM+D1 105 CONTINUE XI(1)=N(1) C *** ITERATION FOR MODEL 3 STARTS HERE 110 CONTINUE DO 115 I=2,K 115 XI(I)=SHAT*XI(I-1)+N(I) RHO(KS)=FHAT DO 120 J=2,KS I=KS-J+1 THETA(I)=SHAT*THETA(I+1)+1.0D0 120 RHO(I)=FHAT*THETA(I) DO 125 I=1,2 DO 125 J=1,2 125 V(I,J)=0.0D0 DTWO=0.0D0 DFOUR=0.0D0 VTWO=0.0D0 DO 130 I=1,K DR(I)=1.0D0-RHO(I) DTWO=DTWO+THETA(I)*(N(I)-R(I))/DR(I) V(1,1)=V(1,1)+N(I)*THETA(I)/DR(I) IF (I.GE.KS-1) GO TO 130 TPRIM(I)=(THETA(I)-(KS-I+1)*SHAT**(KS-I))/(1.0D0-SHAT) IF (I.GE.KS-2) GO TO 130 TDPRIM(I)=(2*TPRIM(I)-(KS-I+1)*(KS-I)*SHAT**(KS-I-1))/(1.0D0-SHAT) 130 CONTINUE DO 135 I=1,K IF (I.GT.KS-1) GO TO 135 DFOUR=DFOUR+(N(I)-R(I))*TPRIM(I)/DR(I) V(1,2)=V(1,2)+N(I)*TPRIM(I)/DR(I) IF (I.GT.KS-2) GO TO 135 VTWO=VTWO+N(I)*(TDPRIM(I)+(FHAT*TPRIM(I)**2)/DR(I)) 135 CONTINUE DELTA(1)=DONE/FHAT-DTWO DELTA(2)=SUM/SHAT-FHAT*DFOUR ESUM=0.0D0 DO 140 I=1,L 140 ESUM=ESUM+XI(I)*RHO(I+1) IF (K.EQ.KS) GO TO 150 DO 145 J=1,IS 145 ESUM=ESUM+XI(K)*RHO(K+J)*SHAT**(J-1) 150 CONTINUE V(1,1)=V(1,1)/FHAT V(2,1)=V(1,2) V(2,2)=ESUM/SHAT+FHAT*VTWO IF (IS.LE.1) V(2,2)=V(2,2)+N(KS-1)*FHAT**2/(1.0D0-RHO(KS-1)) D11=V(1,1)*V(2,2)-V(1,2)*V(2,1) A=V(2,2)/D11 B=V(1,1)/D11 V(1,2)=-V(1,2)/D11 V(2,1)=V(1,2) V(1,1)=A V(2,2)=B IF (NR.EQ.20) GO TO 170 DO 155 I=1,2 IF (DABS(DIFF(I)).GE.0.1D-7) GO TO 160 155 CONTINUE GO TO 170 160 DIFF(1)=0.0D0 DIFF(2)=0.0D0 DO 165 I=1,2 DO 165 J=1,2 165 DIFF(I)=DIFF(I)+V(I,J)*DELTA(J) FHAT=FHAT+DIFF(1) SHAT=SHAT+DIFF(2) NR=NR+1 GO TO 110 170 CONTINUE TEMP1=(DSQRT(V(1,1)))*100.0 TEMP2=(DSQRT(V(2,2)))*100.0 TEMP3=V(1,2) C *** LIKELIHOOD RATIO TEST OF MODEL 3 VS. MODEL 1 TONE=0.D0 TTWO=0.D0 TTHREE=0.D0 T4=0.D0 DO 175 I=1,K RSAVE(I)=DR(I) TONE=TONE+C(I)*DLOG(F(I)) TTHREE=TTHREE+(N(I)-R(I))*DLOG(DR(I)/(1.D0-R(I)/N(I))) IF (I.EQ.K) GO TO 175 TTWO=TTWO+(T(I+1)-R(I+1))*DLOG(S(I)) 175 CONTINUE FSAVE=FHAT SSAVE=SHAT IF (IS.EQ.0) GO TO 185 DO 180 J=1,IS 180 T4=T4+C(K+J)*DLOG(R(K)*C(K+J)/(N(K)*T(K))) 185 TEST=DLOG(FHAT)*DONE+DLOG(SHAT)*SUM-TONE-TTWO+TTHREE IF (IS.GT.0) TEST=TEST-T4 TEST=-2*TEST IDF=2*K+IS-3 KP=KS+1 KL=KS-1 DO 190 I=1,KP 190 DIFF(I)=1.0D0 C *** ITERATION FOR MODEL 2 STARTS HERE NR=0 XI(1)=N(1) 195 CONTINUE DO 200 I=2,K 200 XI(I)=SHAT*XI(I-1)+N(I) SP(1)=SHAT DO 205 I=2,KS 205 SP(I)=SHAT*SP(I-1) RHO(KS)=F(KS) DO 210 I=1,KL J=KS-I 210 RHO(J)=SHAT*RHO(J+1)+F(J) K1=K IF (IS.EQ.0) K1=K-1 DO 220 I=1,K1 RPRIM(I)=0.D0 RDPRIM(I)=0.D0 M=0 SS=1.0D0 DO 220 J=I,KL RPRIM(I)=RPRIM(I)+F(J+1)*(M+1)*SS IF (J.EQ.KL) GO TO 215 IF (I.GT.KS-2) GO TO 215 RDPRIM(I)=RDPRIM(I)+F(J+2)*(M+2)*(M+1)*SS 215 M=M+1 SS=SS*SHAT 220 CONTINUE DJ=0.0D0 DTWO=0.0D0 ESUM=0.0D0 VJ=0.0D0 VKP=0.0D0 DO 240 J=1,K DR(J)=1.0D0-RHO(J) RB(J)=(N(J)-R(J))/DR(J) DJ=DJ*SHAT+RB(J) DELTA(J)=C(J)/F(J)-DJ VJ=VJ*SHAT*SHAT+N(J)/DR(J) V(J,J)=XI(J)/F(J)+VJ IF (J.LE.KS-2) VKP=VKP+N(J)*(RDPRIM(J)+RPRIM(J)*RPRIM(J)/DR(J)) IF (J.GT.KL) GO TO 240 DTWO=DTWO+RB(J)*RPRIM(J) ESUM=ESUM+XI(J)*RHO(J+1) V(J,KP)=N(J)*RPRIM(J)/DR(J) IF (J.GT.1) V(J,KP)=V(J,KP)+N(J-1)*(1.D0+SHAT*RPRIM(J-1)/DR(J-1)) IF (J.LE.2) GO TO 230 JM=J-2 DO 225 I=1,JM 225 V(J,KP)=V(J,KP)+N(I)*SP(J-I-1)*(J-I+SHAT*RPRIM(I)/DR(I)) 230 CONTINUE IF (J.EQ.K) GO TO 240 JP=J+1 DO 235 H=JP,KS 235 V(J,H)=VJ*SP(H-J) 240 CONTINUE IF (IS.EQ.0) GO TO 260 DO 255 J=1,IS DELTA(K+J)=C(K+J)/F(K+J)-DJ*SP(J) V(K+J,K+J)=XI(K)*SP(J)/F(K+J)+VJ*SP(J)*SP(J) DO 245 H=J,IS 245 V(K+J-1,K+H)=VJ*SP(H+J-1) IF (J.EQ.1) GO TO 255 ESUM=ESUM+XI(K)*SP(J-1)*RHO(K+J) VJJ=0.D0 DO 250 I=1,K 250 VJJ=VJJ+N(I)*SP(K+J-I-1)*(K+J-I+SHAT*RPRIM(I)/DR(I)) V(K+J,KP)=VJJ 255 CONTINUE 260 DELTA(KP)=SUM/SHAT-DTWO J=1 IF (IS.EQ.0) J=0 KM=K+J-1 VJP=N(KM)*(1.D0+SHAT*RPRIM(KM)/DR(KM)) H=KM-1 DO 265 I=1,H 265 VJP=VJP+N(I)*SP(KM-I)*(K+J-I+SHAT*RPRIM(I)/DR(I)) V(K+J,KP)=VJP V(KP,KP)=ESUM/SHAT+VKP IF (IS.LE.1) V(KP,KP)=V(KP,KP)+N(KL)*RPRIM(KL)*RPRIM(KL)/DR(KL) DO 270 I=1,KP DO 270 J=1,KP 270 V(J,I)=V(I,J) CALL INVERT(V,KP) IF (NR.EQ.20) GO TO 295 DO 275 I=1,KP IF (DABS(DIFF(I)).GE.0.1D-7) GO TO 280 275 CONTINUE GO TO 295 280 CONTINUE DO 285 I=1,KP DIFF(I)=0.0D0 DO 285 J=1,KP 285 DIFF(I)=DIFF(I)+V(I,J)*DELTA(J) DO 290 I=1,KS 290 F(I)=F(I)+DIFF(I) SHAT=SHAT+DIFF(KP) NR=NR+1 GO TO 195 295 CONTINUE C *** END OF ITERATION FOR MODEL 2 DO 305 I=1,KP V(I,I)=DABS(V(I,I)) DO 300 J=1,KP V(I,J)=DABS(V(I,J)) CORR(I,J)=V(I,J)/(DSQRT(V(I,I)*V(J,J))) 300 CONTINUE CORR(I,I)=1.0D0 305 CONTINUE C WRITE (6,310) 310 FORMAT (//' ',31X,'RECOVERY RATE F(I) (%)') WRITE (6,330) WRITE (6,315) 315 FORMAT (1X,6X,'I',8X,'ESTIMATE',14X,'STANDARD',14X,'95% CONFIDENC' 1,'E') WRITE (6,320) 320 FORMAT (' ',37X,' ERROR ',16X,'INTERVAL') WRITE (6,325) 325 FORMAT ('+',15X,'________',14X,'_________',13X,'______________'/) 330 FORMAT ('+',15X,60('_')) C********** JTRY=K IF (K.EQ.20) JTRY=K-1 DO 340 I=1,JTRY AAA=F(I)*100. BBB=DSQRT(V(I,I))*100. CCC=AAA-BBB*1.96 DDD=AAA+BBB*1.96 WRITE (6,335) I,AAA,BBB,CCC,DDD 335 FORMAT (' ',6X,I2,7X,F7.3,16X,F7.3,11X,F7.3,' - ',F7.3/) 340 CONTINUE IF (K.LT.20) GO TO 345 I=20 AAA=F(20)*100. BBB=DSQRT(V(K,K))*100. CCC=AAA-BBB*1.96 DDD=AAA+BBB*1.96 WRITE (6,335) I,AAA,BBB,CCC,DDD 345 CONTINUE C********** AVEF=0.0 DO 350 I=1,K AVEF=AVEF+(F(I)*100.0) 350 CONTINUE AVEF=AVEF/K WRITE (6,355) AVEF 355 FORMAT (/' ','ARITHMETIC MEAN RECOVERY RATE = ',F9.2) SUMONE=0.0 SUMTWO=0.0 DO 360 I=1,K DO 360 J=1,K SUMONE=SUMONE+V(I,J) 360 CONTINUE VARAVF=(SUMONE+2*SUMTWO)/(K )**2 SDAVF=DSQRT(VARAVF)*100. LOWAVF=AVEF-(1.96*SDAVF) HIAVF=AVEF+(1.96*SDAVF) WRITE (6,370) SDAVF 370 FORMAT (' ','STANDARD ERROR OF MEAN RECOVERY RATE = ',F9.2) WRITE (6,375) LOWAVF,HIAVF 375 FORMAT (' ','95% CONFIDENCE INTERVAL FOR MEAN RECOVERY RATE = ',F8 1.2,' - ',F8.2/) SHAT=SHAT*100. WRITE (6,380) SHAT 380 FORMAT (//' ','CONSTANT SURVIVAL RATE (%) = ',F7.2) SDSHAT=DSQRT(V(KS+1,KS+1))*100.0 LOWS=SHAT-1.96*SDSHAT HIS=SHAT+1.96*SDSHAT WRITE (6,385) SDSHAT 385 FORMAT (' ','STANDARD ERROR OF THE CONSTANT SURVIVAL RATE = ',F9.2 1) WRITE (6,390) LOWS,HIS 390 FORMAT (1X,'95% CONFIDENCE INTERVAL FOR THE CONSTANT SURVIVAL RAT' 1,'E',F8.2,' - ',F8.2//) SHAT=SHAT/100. MLS=-1.0/DLOG(SHAT) SDMLS=(MLS**2*(SDSHAT/100.))/(SHAT) LOWMLS=-1.0/DLOG(LOWS/100.) HIMLS=-1.0/DLOG(HIS/100.) IF (HIS.GT.100.0) HIMLS=99999999. WRITE (6,395) MLS 395 FORMAT (//' ','MEAN LIFE SPAN AS AN ADULT = ',F7.2) WRITE (6,400) SDMLS 400 FORMAT (' ','STANDARD ERROR OF THE MEAN LIFE SPAN = ',F7.2) IF(HIMLS.NE.99999999.) GO TO 404 WRITE (6,402) LOWMLS 402 FORMAT(' ','95% CONFIDENCE INTERVAL OF LIFE SPAN ', 1F7.2,' - *******'//) GO TO 407 404 WRITE (6,405) LOWMLS,HIMLS 405 FORMAT (' ','95% CONFIDENCE INTERVAL OF LIFE SPAN ',F7.2,' - ',F7. 12//) 407 CONTINUE TONE=0.D0 TSIX=0.D0 TFOUR=0.D0 TFIVE=0.D0 DO 410 I=1,K TONE=TONE+C(I)*DLOG(F(I)*N(I)*T(I)/(R(I)*C(I))) TSIX=TSIX+(N(I)-R(I))*DLOG(N(I)*DR(I)/(N(I)-R(I))) TFOUR=TFOUR+C(I)*DLOG(F(I)) 410 TFIVE=TFIVE+(N(I)-R(I))*DLOG(RSAVE(I)/DR(I)) TEST2=DLOG(FSAVE)*DONE-TFOUR+TFIVE+SUM*DLOG(SSAVE/SHAT) TEST3=TONE+SUM*DLOG(SHAT)-TTWO+TSIX IF (IS.EQ.0) GO TO 420 DO 415 J=1,IS TEST2=TEST2-C(K+J)*DLOG(F(K+J)) 415 TEST3=TEST3+C(K+J)*DLOG(F(K+J)*N(K)*T(K)/(R(K)*C(K+J))) 420 TEST2=-2*TEST2 TEST3=-2*TEST3 IDF2=KL IDF3=K-2 DO 425 JJ=1,K DO 425 II=JJ,KS EE(II,JJ)=0.0 CHISQS(II,JJ)=0.0 EE(JJ,II)=N(JJ)*SHAT**(II-JJ)*F(II) 425 CONTINUE DO 445 I=1,K DO 430 M=1,KS IF (I.EQ.K.AND.M.EQ.K.AND.K.EQ.KS) GO TO 445 IF (EE(I,KS-M+1).GT.2.0) GO TO 435 EE(I,KS)=EE(I,KS)+EE(I,KS-M) EE(I,KS-M)=0.0 IF (EE(I,KS).GT.2.0) GO TO 440 430 CONTINUE 435 ISAVE(I)=M-1 GO TO 445 440 ISAVE(I)=M C ISAVE(I)=NUMBER OF DEGREES OF FREEDOM LOST IN ROW I 445 CONTINUE WRITE (6,450) 450 FORMAT (///1X,'MATRIX OF EXPECTED VALUES -- ASSUMING A CONSTANT S' 1,'URVIVAL RATE AND TIME-SPECIFIC RECOVERY RATES (MODEL 2)'/) DO 460 I=1,K ISS=(I-1)*6+11 WRITE(FORMT,455) ISS WRITE (6,FORMT) (EE(I,JJ),JJ=I,KS) 455 FORMAT ('(',I2,'X,20F6.1)') 460 CONTINUE DO 470 I=1,K MM=ISAVE(I) DO 465 M=1,MM IF (ISAVE(I).NE.0) ID(I,KS)=ID(I,KS)+ID(I,KS-M) IF (ISAVE(I).NE.0) ID(I,KS-M)=0.0 465 CONTINUE 470 CONTINUE DO 475 I=1,K DO 475 JJ=1,KS IF (EE(I,JJ).GT.0.0) CHISQS(I,JJ)=((ID(I,JJ)-EE(I,JJ))**2)/EE(I,JJ 1) IF (EE(I,JJ).LT..01) CHISQS(I,JJ)=0.0 CHISUM=CHISUM+CHISQS(I,JJ) 475 CONTINUE WRITE (6,480) 480 FORMAT (//1X,'MATRIX OF CHI-SQUARE VALUES--ASSUMING A CONSTANT SU' 1,'RVIVAL RATE AND TIME-SPECIFIC RECOVERY RATE (MODEL 2)'/) DO 490 I=1,K ISS=(I-1)*6+11 WRITE(FORMT,485) ISS 485 FORMAT ('(',I2,'X,20F6.2)') WRITE (6,FORMT) (CHISQS(I,JJ),JJ=I,KS) 490 CONTINUE SUMDF=0.0 DO 495 I=1,K SUMDF=SUMDF+ISAVE(I) 495 CONTINUE IF (SUMDF.GT.0.0) WRITE (6,500) 500 FORMAT (/1X,'(FREQUENCIES WERE COMBINED WHERE EXPECTED VALUES WER' 1,'E LESS THAN 2.0)'/) IIDF=0.0 IIDF=(K*(K+1))/2+K*(KS-K)-(K+1)-(KS-K) IIDF=IIDF-SUMDF DO 580 I=1,K POOLR=0.0 POOLE=0.0 DO 585 M=1,KS POOLR=POOLR+ID(I,M) 585 POOLE=POOLE+EE(I,M) MORCHI=POOLR-POOLE MORCHI=MORCHI*MORCHI/(N(I)-POOLE) 580 CHISUM=CHISUM+MORCHI WRITE (6,505) 505 FORMAT (//1X,'TEST OF THE NULL HYPOTHESIS THAT THE DATA FIT THE M' 1,'ODEL ASSUMING A CONSTANT SURVIVAL AND TIME-SPECIFIC RECOVERY ', 2'RATES'/) WRITE (6,510) CHISUM 510 FORMAT (' ',' CHI-SQUARE VALUE (SAMPLE) = ',F8.2) WRITE (6,515) 515 FORMAT (' ',' THEORETICAL CHI-SQUARE VALUE AT THE 5% LEVEL = ') IF (IIDF-100) 520,520,535 520 WRITE (6,525) THCHI(IIDF) 525 FORMAT ('+',49X,F8.2) WRITE (6,530) IIDF 530 FORMAT (' ',' DEGREES OF FREEDOM = ',I6) GO TO 545 535 THECHI=.5*(1.64+SQRT((2.*IIDF)-1.))**2 WRITE (6,540) THECHI 540 FORMAT ('+',49X,F8.2) WRITE (6,530) IIDF 545 CONTINUE CALL CHI1(IIDF,CHISUM) IF (K.EQ.20) GO TO 567 WRITE (6,550) 550 FORMAT (///1X,' ESTIMATED COVARIANCE AND CORRELATION COEFFICIENTS' 1,':'//) WRITE (6,555) 555 FORMAT (1X,7X,'I',7X,6X,'COVAR(F(I),S)',11X,'CORR(F(I),S)',10X,'C' 1,'OVAR(F(I),F(I+1))',6X,'CORR(F(I),F(I+1))') L=K-1 COVFAS=0.0 DO 565 I=1,L WRITE (6,560) IYEAR(I),V(I,KP),CORR(I,KP),V(I,I+1),CORR(I,I+1) 560 FORMAT (1H ,5X,I4,4(12X,F12.9)) COVFAS=COVFAS+V(I,KP) 565 CONTINUE 567 QQ=100000.0 QQQ=1000000.0 IF(K.LT.KS) QQ=V(K,K+1) IF(K.LT.KS) QQQ=CORR(K,K+1) WRITE (6,560) IYEAR(K),V(K,KP),CORR(K,KP),QQ,QQQ IF(K.EQ.KS) WRITE (6,575) 575 FORMAT(/1X,'(QUANTITIES SHOWN AS ********** ARE NOT ESTIMABLE UND' 1,'ER MODEL 2)') COVFAS=COVFAS+V(K,KP) COVFAS=10000.*COVFAS/K CORFAS=COVFAS/(SDAVF*SDSHAT) COVFAS=COVFAS/10000. WRITE (6,561) CORFAS 561 FORMAT(//,1X,'CORR(SURVIVAL RATE, AVE RECOVERY RATE) = ',F10.7) WRITE (6,562) COVFAS 562 FORMAT(1X,'COVAR(SURVIVAL RATE, AVE RECOVERY RATE) = ',F10.7) 566 CONTINUE WRITE (6,570) NR 570 FORMAT (//' ',2X,I6,2X,'ITERATIONS') RETURN END SUBROUTINE MOD2W1 WRITE (6,10) 10 FORMAT ('0','MODEL 2') WRITE (6,15) 15 FORMAT (' ','********'/) WRITE (6,20) 20 FORMAT (1X,'ANALYSIS ASSUMING CONSTANT SURVIVAL BUT TIME-SPECIFIC' 1,' RECOVERY RATES'/) WRITE (6,25) 25 FORMAT (1X,'THIS MODEL WAS DEVELOPED BY DRS. BROWNIE AND ROBSON A' 1,'T THE CORNELL BIOMETRICS UNIT.') WRITE (6,30) 30 FORMAT (//' ','SPECIFICALLY, THE MODEL STRUCTURE IS: '/) WRITE (6,35) 35 FORMAT (1X,' F(1) SF(2) SSF(3) SSSF(4) SSS' 1,'SF(5)') WRITE (6,40) 40 FORMAT (1X,' F(2) SF(3) SSF(4) SS' 1,'SF(5)') WRITE (6,45) 45 FORMAT (1X,' F(3) SF(4) S' 1,'SF(5)') WRITE (6,50) 50 FORMAT (1X,' F(4) ' 1,'SF(5)'///) RETURN END