SUBROUTINE MODLA2(TITLE,CHIA,IDFA) IMPLICIT DOUBLE PRECISION (A-Z) INTEGER IDIM,IDIM2 PARAMETER (IDIM=50,IDIM2=IDIM*IDIM) INTEGER I,J,K,IDFA,LM1,LM2,LAST,IP1,IM1,JP1,T,MAXITR,YR(IDIM), , NN,IPRNT CHARACTER*128 TITLE,UNDEF*6 DIMENSION K(IDIM),PVARM1(IDIM),PVARM2(IDIM),PVARM3(IDIM) DIMENSION PHIY(IDIM),PHI(IDIM),XM(IDIM),P(IDIM),SEPHI(IDIM), , PM1(IDIM),SEPHIY(IDIM),ESTN(IDIM),SDP(IDIM) DIMENSION SEPHY2(IDIM),SEPHI2(IDIM),SDN2(IDIM) DIMENSION SDN(IDIM),ACOV(IDIM,IDIM), , ESTB0(IDIM), ESTB (IDIM,IDIM) DIMENSION COVAA(IDIM),COVYA(IDIM),COVYA1(IDIM) COMMON /BLK1/M(IDIM),N(IDIM),NY(IDIM),B(IDIM),R(IDIM),RY(IDIM), , S(IDIM),SY(IDIM),YR COMMON /BLK3/LAST,AVFI,FIMD,AVFY,FYMD,EPSLON,MAXITR,IPRNT COMMON /BLK2/CA(IDIM,2,2),CM(IDIM,2,2),CN(IDIM,2,2),Z(IDIM,2,2) DATA AVN,AVP,AVB,AVPM1,AVPM2,AVPM3,AVPM1V,AVPM2V,AVPM3V, , AVSEN,AVSEFI,AVSEFY,ASEFI2,ASEFY2,AVSEN2, , AVSEP,AVSEB/17*.0D0/ DATA XM,PVARM3/IDIM*.0D0,IDIM*.0D0/ DATA PM1,ESTN,PHI,PHIY,P,SEPHI,SEPHIY,SEPHY2,SEPHI2,SDN2, ,SDN,ESTB0,COVAA,COVYA,COVYA1/IDIM*.12345D55, ,IDIM*.12345D55,IDIM*.12345D55,IDIM*.12345D55,IDIM*.12345D55, ,IDIM*.12345D55,IDIM*.12345D55,IDIM*.12345D55,IDIM*.12345D55, ,IDIM*.12345D55,IDIM*.12345D55,IDIM*.12345D55,IDIM*.12345D55, ,IDIM*.12345D55,IDIM*.12345D55/ DATA ACOV,ESTB/IDIM2*.12345D55,IDIM2*.12345D55/ DATA UNDEF/'Div/0.'/ C U(I)=ESTN(I)-XM(I) Q(I)=.1D1-P(I) AVFI=.0D0 AVFY=.0D0 LM1=LAST-1 LM2=LAST-2 XLM2=LM2 C COMPUTE ESTIMATORS OF C # MARKED ADULTS - XM(I), C CAPTURE Probability - P(I),PY(I) C # IN POPULATION - ESTN(I),ESTNY(I) C AVG. ESTIMATES - AVM,AVP,AVN C C -------------- COMPUTE Variances --------- C Variances OF ADULT AND YOUNG SURVIVAL RATES ARE COMPUTED USING C POLLOCKS FORMULA FOR "EXCLUDING NON-SAMPLING ERROR Variance", C ( ADULTS - SEPHI2(I), YOUNG - SEPHY2(I) ) C AND FOR "INCLUDING NON-SAMPLING ERROR Variance". C ( ADULTS - SEPHI(I), YOUNG - SEPHIY(I) ) C DO 20 I = 1,LM2 J=I+1 XM(J)=M(J)+B(J)*(S(J)+.1D1)/(R(J)+.1D1) P(J)=M(J)/XM(J) ESTN(J)=(N(J)+.1D1)*XM(J)/(M(J)+.1D1) IF(ESTN(J).LT.N(J))ESTN(J)=N(J) T7=.12345D55 IF(R(J)*S(J).GT..0D0)T7=.1D1/R(J)-.1D1/S(J) C C (L) (L+1) (L) (L+1) C PT1=T PT2=T PM1=M PM2(J)=M C I+1 I+1 I+1 I+1 C PT1=CM(J,2,2)+Z(J,2,2) IF(PT1.LE..0D0)PT1=RY(I) PT2=M(J)+B(J)-PT1 IF(PT1+PT2.NE..0D0)PM1(J)=PT1/(PT1+PT2)*XM(J) PM2=XM(J)-PM1(J) PC=PM1(J)*PM2*(XM(J)-PT1-PT2)/(XM(J)*(PT1+PT2)) PVARM3(J)=(XM(J)-M(J))*(XM(J)-M(J)+S(J))*T7 PVARM1(J)=PM1(J)/XM(J)*PM1(J)/XM(J)*PVARM3(J)+PC IF(SY(I).LE..0D0)GO TO 10 PHIY(I)=PM1(J)/SY(I) IF(PM1(J).LE..0D0)GO TO 10 PVARFY=PHIY(I)*PHIY(I)*PVARM1(J)/PM1(J)/PM1(J) SEPHY2(I)=XSQRT(PVARFY) SEPHIY(I)=XSQRT(PVARFY+PHIY(I)*(.1D1-PHIY(I))/SY(I)) 10 IF(S(I)*R(I).LE..0D0)GO TO 20 PHI(I)=PM2/(S(I)*B(I)/R(I)+S(I)) PVM3I=(XM(I)-M(I))*(XM(I)-M(I)+S(I))*(.1D1/R(I)-.1D1/S(I)) PVARM2(J)=PM2/XM(J)*PM2/XM(J)*PVARM3(J)+PC PVARFI=PHI(I)**2*(PVARM2(J)/PM2**2+ + PVM3I/(XM(I)-M(I)+S(I))**2) SEPHI2(I)=XSQRT(PVARFI) SEPHI(I)=XSQRT(PVARFI+PHI(I)*(.1D1-PHI(I))/ / (S(I)*B(I)/R(I)+S(I))) C IF(SY(I).GT..0D0) ) COVYA(I)=(PM1(J)*PM2/(PM1(J)+PM2)**2*PVARM3(J)-PC)/ / (SY(I)*(XM(I)-M(I)+S(I))) IF(I.LE.1)GO TO 20 COVAA(I-1)=-PHI(I)*PHI(I-1)*(XM(I)-M(I))/XM(I)* * (.1D1/R(I)-.1D1/S(I)) IF(PHI(I-1).NE..0D0)COVYA1(I-1)=COVAA(I-1)/PHI(I-1)*PHIY(I-1) 20 CONTINUE C 100 WRITE(7,110)CHAR(12),TITLE 110 FORMAT(A1,A128/ //' Model A2 - Time-specific survival rates for adults and young, ' 2,'time-specific adult capture probabilities.'/ 3/' Parameters: PHI(i) = Adult survival rate in period i'/ 5 12X,6H PHI','(i) = Young survival rate in period i'/ 6 ' P(i) = Adult capture probability in period i'/) WRITE(7,120) 120 FORMAT(//' i',6X,4HPHI',5X,8HSE(PHI'),5X,9HSE'(PHI'), ,7X,'95% confidence interval'/1X,73('-')) C C PRINT ESTIMATES FOR YOUNG C NN=0 DO 140 I = 1,LM2 IF(PHIY(I).GE..1D50)THEN WRITE(7,129)YR(I),(UNDEF,J=1,5) 129 FORMAT(1X,I4,A9,A13,A14,A18,' - ',A7) ELSE AVFY=AVFY+PHIY(I) AVSEFY=AVSEFY+SEPHIY(I)*SEPHIY(I) ASEFY2=ASEFY2+SEPHY2(I)*SEPHY2(I) NN=NN+1 CI1=PHIY(I)-.196D1*SEPHIY(I) CI2=PHIY(I)+.196D1*SEPHIY(I) WRITE(7,130)YR(I),PHIY(I),SEPHIY(I),SEPHY2(I),CI1,CI2 130 FORMAT(1X,I4,F9.4,F13.4,F14.4,F18.4,' - ',F7.4) ENDIF IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 140 CONTINUE C C COMPUTE MEAN ESTIMATES C IF(NN.GT.0)THEN AVFY=AVFY/FLOAT(NN) AVSEFY=XSQRT(AVSEFY)/FLOAT(NN) ASEFY2=XSQRT(ASEFY2)/FLOAT(NN) CI1=AVFY-.196D1*AVSEFY CI2=AVFY+.196D1*AVSEFY WRITE(7,150)AVFY,AVSEFY,ASEFY2,CI1,CI2 150 FORMAT(1X,73('-')/' MEAN',F9.4,F13.4,F14.4,F18.4,' - ',F7.4/ / 1X,73('-')) ENDIF C C --------- WRITE STATS FOR ADULTS ------------------------ C WRITE(7,170) 170 FORMAT(//' I',6X,4HPHI ,5X,8HSE(PHI ),5X,9HSE'(PHI ), ,7X,'95% confidence interval'/1X,73('-')) NN=0 COVAA(LM2)=.0D0 DO 180 I=1,LM2 IF(PHI(I).GE..1D50)THEN WRITE(7,179)YR(I),(UNDEF,J=1,5) 179 FORMAT(1X,I4,A9,A13,A14,A18,' - ',A7) ELSE NN=NN+1 AVFI=AVFI+PHI(I) AVSEFI=AVSEFI+SEPHI(I)*SEPHI(I)+.2D1*COVAA(I) ASEFI2=ASEFI2+SEPHI2(I)*SEPHI2(I)+.2D1*COVAA(I) CI1=PHI(I)-.196D1*SEPHI(I) CI2=PHI(I)+.196D1*SEPHI(I) WRITE(7,130)YR(I),PHI(I),SEPHI(I),SEPHI2(I),CI1,CI2 ENDIF IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 180 CONTINUE IF(NN.GT.0)THEN AVFI=AVFI/FLOAT(NN) AVSEFI=XSQRT(AVSEFI)/FLOAT(NN) ASEFI2=XSQRT(ASEFI2)/FLOAT(NN) CI1=AVFI-.196D1*AVSEFI CI2=AVFI+.196D1*AVSEFI WRITE(7,150)AVFI,AVSEFI,ASEFI2,CI1,CI2 ENDIF WRITE(7,190) 190 FORMAT(//' I M1 SE(M1) 95% conf. interval', ,' M2 SE(M2) 95% conf. interval', ,' M SE(M) 95% conf. interval'/1X,132('-')) NN=0 DO 210 I=2,LM1 AVPM3=AVPM3+XM(I) SEM3=XSQRT(PVARM3(I)) AVPM3V=AVPM3V+PVARM3(I) CI1M3=XM(I)-.196D1*SEM3 CI2M3=XM(I)+.196D1*SEM3 IF(PM1(I).GE..1234D55)THEN WRITE(7,191)YR(I),(UNDEF,J=1,8),XM(I),SEM3,CI1M3,CI2M3 191 FORMAT(1X,I4,A9,A11,A10,'-',A9,A12,A11,A10,'-',A9, , F12.2,F11.2,F10.2,'-',F9.2) ELSE NN=NN+1 SEM1=XSQRT(PVARM1(I)) SEM2=XSQRT(PVARM2(I)) AVPM1=AVPM1+PM1(I) PM2=XM(I)-PM1(I) AVPM2=AVPM2+PM2 AVPM1V=AVPM1V+PVARM1(I) AVPM2V=AVPM2V+PVARM2(I) CI1M1=PM1(I)-.196D1*SEM1 CI2M1=PM1(I)+.196D1*SEM1 CI1M2=PM2-.196D1*SEM2 CI2M2=PM2+.196D1*SEM2 WRITE(7,200)YR(I),PM1(I),SEM1,CI1M1,CI2M1,PM2,SEM2,CI1M2, , CI2M2,XM(I),SEM3,CI1M3,CI2M3 200 FORMAT(1X,I4,F9.2,F11.2,F10.2,'-',F9.2, , 2(F12.2,F11.2,F10.2,'-',F9.2)) ENDIF IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 210 CONTINUE AVPM3=AVPM3/XLM2 SEM3=XSQRT(AVPM3V)/XLM2 CI1M3=AVPM3-.196D1*SEM3 CI2M3=AVPM3+.196D1*SEM3 IF(NN.GT.0)THEN SEM1=XSQRT(AVPM1V)/FLOAT(NN) SEM2=XSQRT(AVPM2V)/FLOAT(NN) AVPM1=AVPM1/FLOAT(NN) AVPM2=AVPM2/FLOAT(NN) CI1M1=AVPM1-.196D1*SEM1 CI2M1=AVPM1+.196D1*SEM1 CI1M2=AVPM2-.196D1*SEM2 CI2M2=AVPM2+.196D1*SEM2 WRITE(7,220)AVPM1,SEM1,CI1M1,CI2M1,AVPM2,SEM2,CI1M2,CI2M2, , AVPM3,SEM3,CI1M3,CI2M3 220 FORMAT(1X,132('-')/' MEAN',F9.2,F11.2,F10.2,'-',F9.2, , 2(F12.2,F11.2,F10.2,'-',F9.2)/1X,132('-')) ELSE WRITE(7,221)AVPM1,SEM1,CI1M1,CI2M1,(UNDEF,J=1,8) 221 FORMAT(1X,132('-')/' MEAN',F9.2,F11.2,F10.2,'-',F9.2, , 2(A12,A11,A10,'-',A9)/1X,132('-')) ENDIF C COMPUTE # ADULTS ENTERING POP. IN TIME I, C & ALIVE (AS ADULTS) IN TIME J. = ESTB(I,J) C = ESTB0(J) FOR I=0 ESTN(1)=ESTN(2) ESTB0(1)=ESTN(1) DO 40 I = 2,LAST IF(PHI(I-1).GE..1234D55)GO TO 40 ESTB0(I)=PHI(I-1)*ESTB0(I-1) ESTB(I,I-1)=ESTN(I)-(ESTN(I-1)+N(I)-S(I))*PHI(I-1) IF (I.GE.LAST) GO TO 40 DO 30 J = I+1,LAST 30 IF(PHI(J-1).LT..1234D55)ESTB(J,I-1)=PHI(J-1)*ESTB(J-1,I-1) 40 CONTINUE C COMPUTE Variances & Covariances ON ESTN(I) - ADULTS C DO 90 I = 2,LM1 IF(ESTN(1).EQ..0D0)GO TO 90 VSUM1=ESTB0(I)*ESTB0(I)/ESTN(1) IM1=I-1 DO 50 J = 1,IM1 IF(ESTB(J+1,J).GE..1234D55)THEN VSUM1=.1234D55 ELSE VSUM1=VSUM1+ESTB(I,J)*ESTB(I,J)/ESTB(J+1,J) ENDIF 50 CONTINUE T1=XM(I) T2=S(I)-M(I) IF(R(I)*N(I)*S(I)*XM(I).EQ..0D0)THEN VARN=.12345D55 SDN(I)=.12345D55 SDN2(I)=.12345D55 ELSE T3=.1D1/R(I)-.1D1/S(I) T4=ESTN(I)-N(I) ALFA=M(I)/N(I) VARN=ESTN(I)*T4*((T1+T2)/T1*T3+(.1D1-ALFA)/M(I)) IF(VARN.GE..0D0)THEN SDN2(I)=XSQRT(VARN) AVSEN2=AVSEN2+VARN ENDIF IF(VSUM1.GE..1234D55)THEN SDN(I)=.12345D55 AVSEN=.12345D55 ELSE VARN=VARN+ESTN(I)-VSUM1 IF(VARN.GT..0D0)THEN AVSEN=AVSEN+VARN SDN(I)=XSQRT(VARN) ELSE SDN(I)=.0D0 ENDIF ENDIF ENDIF 60 IF(ESTB0(1).EQ..0D0.OR.ESTB0(1).GE..1234D55)GO TO 90 DO 80 T = I+1,LAST IF(ESTB0(T).GE..1234D55.OR.ESTB0(I).GE..1234D55)GO TO 80 ACOV(I,T)=ESTB0(T)*(.1D1-ESTB0(I)/ESTB0(1)) DO 70 J = 1,IM1 IF(ESTB(J+1,J).EQ..0D0.OR.ESTB(J+1,J).GE..1234D55.OR. . ESTB(T,J).GE..1234D55.OR.ESTB(I,J).GE..1234D55)THEN ACOV(I,T)=.12345D55 AVSEN=.12345D55 ELSE ACOV(I,T)=ACOV(I,T)+ESTB(T,J)*(.1D1-ESTB(I,J)/ESTB(J+1,J)) ENDIF 70 CONTINUE IF(AVSEN.LT..1234D55.AND.ACOV(I,T).LT..1234D55) . AVSEN=AVSEN+ACOV(I,T)*.2D1 80 CONTINUE 90 CONTINUE WRITE(7,230) 230 FORMAT(//' I',6X,4H N ,5X,8HSE( N ),5X,9HSE'( N ), ,7X,'95% confidence interval'/1X,73('-')) NN=0 DO 250 I=2,LM1 IF(ESTN(I).GE..1234D55)THEN WRITE(7,237)YR(I),(UNDEF,J=1,5) 237 FORMAT(1X,I4,A9,A13,2A14,' - ',A14) ELSE NN=NN+1 AVN=AVN+ESTN(I) IF(SDN(I).GE..1234D55.AND.SDN2(I).LT..1234D55)THEN WRITE(7,238)YR(I),ESTN(I),UNDEF,SDN2(I),UNDEF,UNDEF 238 FORMAT(1X,I4,F9.2,A13,F14.2,A14,' - ',A14) ELSEIF(SDN2(I).GE..1234D55)THEN WRITE(7,239)YR(I),ESTN(I),(UNDEF,J=1,4) 239 FORMAT(1X,I4,F9.2,A13,2A14,' - ',A14) ELSE CI1=ESTN(I)-.196D1*SDN(I) CI2=ESTN(I)+.196D1*SDN(I) WRITE(7,240)YR(I),ESTN(I),SDN(I),SDN2(I),CI1,CI2 240 FORMAT(1X,I4,F9.2,F13.2,2F14.2,' - ',F14.2) ENDIF ENDIF IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 250 CONTINUE IF(NN.GT.0)THEN AVN=AVN/FLOAT(NN) IF(AVSEN.GE..1234D55.AND.AVSEN2.GE..1234D55)THEN WRITE(7,258)AVN,(UNDEF,J=1,4) 258 FORMAT(1X,73('-')/' MEAN',F9.2/1X,73('-')) ELSEIF(AVSEN2.LT..1234D55.AND.AVSEN.GE..1234D55)THEN WRITE(7,259)AVN,UNDEF,AVSEN2,UNDEF,UNDEF 259 FORMAT(1X,73('-')/' MEAN',F9.2,A13,F14.2,A14,' - ',A14/ / 1X,73('-')) ELSE AVSEN=XSQRT(AVSEN)/FLOAT(NN) AVSEN2=XSQRT(AVSEN2)/FLOAT(NN) CI1=AVN-.196D1*AVSEN CI2=AVN+.196D1*AVSEN WRITE(7,260)AVN,AVSEN,AVSEN2,CI1,CI2 260 FORMAT(1X,73('-')/' MEAN',F9.2,F13.2,2F14.2,' - ',F14.2/ / 1X,73('-')) ENDIF ENDIF WRITE(7,270) 270 FORMAT(//' I',6X,4H p ,5X,8HSE( P ), ,7X,'95% confidence interval'/1X,59('-')) DO 290 I=2,LM1 AVP=AVP+P(I) IF(R(I)*S(I)*B(I)*M(I).LE..0D0)THEN SDP(I)=.12345D55 AVSEP=.12345D55 WRITE(7,279)YR(I),P(I),(UNDEF,J=1,3) 279 FORMAT(1X,I4,F9.4,A13,A18,' - ',A7) ELSE SDP(I)=P(I)*(.1D1-P(I))*XSQRT(.1D1/R(I)-.1D1/S(I)+.1D1/B(I)+ + .1D1/M(I)) CI1=P(I)-.196D1*SDP(I) CI2=P(I)+.196D1*SDP(I) AVSEP=AVSEP+SDP(I)*SDP(I) WRITE(7,280)YR(I),P(I),SDP(I),CI1,CI2 280 FORMAT(1X,I4,F9.4,F13.4,F18.4,' - ',F7.4) ENDIF IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 290 CONTINUE AVP=AVP/XLM2 IF(AVSEP.GE..1234D55)THEN WRITE(7,299)AVP 299 FORMAT(1X,59('-')/' MEAN',F9.4/1X,59('-')) ELSE AVSEP=XSQRT(AVSEP)/XLM2 CI1=AVP-.196D1*AVSEP CI2=AVP+.196D1*AVSEP WRITE(7,300)AVP,AVSEP,CI1,CI2 300 FORMAT(1X,59('-')/' MEAN',F9.4,F13.4,F18.4,' - ',F7.4/ / 1X,59('-')) ENDIF WRITE(7,310) 310 FORMAT(//' I',6X,4H B ,5X,8HSE( B ), ,7X,'95% confidence interval'/1X,59('-')) NN=0 DO 340 I=2,LM2 COVB=.0D0 IF(R(I)*S(I+1).EQ..0D0.OR.PHI(I).GE..1234D55)THEN COVFP=.12345D55 ELSE COVFP=P(I)*Q(I)*PHI(I)*(.1D1/R(I)-.1D1/S(I+1)) ENDIF IF(I.EQ.2)GO TO 320 IF(R(I)*S(I).EQ..0D0.OR.PHI(I-1).GE..1234D55)THEN COVFF1=.12345D55 ELSE COVFF1=-PHI(I-1)*Q(I)*PHI(I)*(.1D1/R(I)-.1D1/S(I)) ENDIF IF(R(I)*S(I).EQ..0D0.OR.PHIY(I-1).GE..1234D55)THEN COVFF1=.12345D55 ELSE CVFFY1=-PHIY(I-1)*Q(I)*PHI(I)*(.1D1/R(I)-.1D1/S(I)) ENDIF IF(PHI(I).GE..1234D55.OR.P(I).EQ..0D0.OR.COVFF1.GE..12345D55 . .OR.COVFP.GE..1234D55.OR.SDP(I).GE..1234D55)THEN COVB=.12345D55 ELSE COVB=-PHI(I)*Q(I)*U(I)/P(I)+Q(I-1)*U(I-1)*(Q(I)*U(I)* * COVFF1-PHI(I)*U(I)/P(I)*COVFP1)+U(I)/P(I)*(Q(I)*U(I)*COVFP- - PHI(I)*U(I)/P(I)*SDP(I)**2)+SY(I-1)*(-Q(I)*U(I)*CVFFY1+ + PHI(I)*U(I)/P(I)*CVFYP1) ENDIF 320 IF(PHI(I).GE..1234D55.OR.R(I+1)*S(I+1)*B(I+1).EQ..0D0)THEN COVFP1=.12345D55 COVFFY=.12345D55 ELSE COVFP1=-PHI(I)*P(I+1)*Q(I+1)**2*(.1D1/R(I+1)-.1D1/S(I+1)+ + (S(I+1)-R(I+1))/S(I+1)/B(I+1)) IF(PHIY(I).GE..1234D55)THEN COVFFY=.12345D55 ELSE COVFFY=PHI(I)*PHIY(I)*Q(I+1)**2*(.1D1/R(I+1)-.1D1/S(I+1)+ + ((S(I+1)-R(I+1))/S(I+1))**2*M(I+1)/B(I+1)/(M(I+1)+B(I+1))) ENDIF ENDIF IF(PHIY(I).GE..1234D55.OR.R(I+1)*S(I+1)*B(I+1).EQ..0D0)THEN CVFYP1=.12345D55 ELSE CVFYP1=-PHIY(I)*P(I+1)*Q(I+1)**2*(.1D1/R(I+1)-.1D1/S(I+1)+ + (S(I+1)-R(I+1))/S(I+1)/B(I+1)) ENDIF IF(ESTB(I+1,I).GE..1234D55)THEN AVSEB=.12345D55 WRITE(7,328)YR(I),(UNDEF,J=1,4) 328 FORMAT(1X,I4,A9,A13,A18,' - ',A10) ELSE AVB=AVB+ESTB(I+1,I) NN=NN+1 IF(P(I+1)*P(I).EQ..0D0.OR.PHI(I)+SDP(I+1)+SEPHI(I)+COVFP1+ + COVFP+SEPHIY(I)+COVFFY+CVFYP1.GE..1234D55)THEN AVSEB=.12345D55 WRITE(7,329)YR(I),ESTB(I+1,I),(UNDEF,J=1,3) 329 FORMAT(1X,I4,F9.2,A13,A18,' - ',A10) ELSE VARB=U(I+1)*Q(I+1)/P(I+1)+PHI(I)*Q(I)*U(I)*(.1D1+PHI(I)*Q(I)/ / P(I))+(U(I+1)/P(I+1)*SDP(I+1))**2+(PHI(I)*U(I)/P(I))**2+ + (Q(I)*U(I)*SEPHI(I))**2+.2D1*Q(I)*U(I)/U(I+1)/P(I+1)*COVFP1- - .2D1*PHI(I)*Q(I)*U(I)**2/P(I)*COVFP+(SY(I)*SEPHIY(I))**2- - .2D1*SY(I)*(Q(I)*U(I)*COVFFY+U(I+1)/P(I+1)*CVFYP1) SDB=XSQRT(VARB) CI1=ESTB(I+1,I)-.196D1*SDB CI2=ESTB(I+1,I)+.196D1*SDB AVSEB=AVSEB+SDB*SDB+COVB*.2D1 WRITE(7,330)YR(I),ESTB(I+1,I),SDB,CI1,CI2 330 FORMAT(1X,I4,F9.2,F13.2,F18.2,' - ',F10.2) ENDIF ENDIF IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 340 CONTINUE IF(NN.GT.0)THEN AVB=AVB/FLOAT(NN) IF(AVSEB.LT..1234D55)THEN AVSEB=XSQRT(AVSEB)/(XLM2-.1D1) CI1=AVB-.196D1*AVSEB CI2=AVB+.196D1*AVSEB WRITE(7,350)AVB,AVSEB,CI1,CI2 350 FORMAT(1X,59('-')/' MEAN',F9.2,F13.2,F18.2,' - ',F10.2/ / 1X,59('-')) ELSE WRITE(7,351)AVB 351 FORMAT(1X,59('-')/' MEAN',F9.2/1X,59('-')) ENDIF ENDIF C C WRITE Covariance MATRICES C IF(IPRNT.LE.1)THEN WRITE(7,160)CHAR(12),TITLE 160 FORMAT(A1,A128/) WRITE(7,360) 360 FORMAT(//40X,'Model A2 Covariances of survival rates'// / ' i Adlt(i)-Adlt(i+1) ', , ' Yng(i)-Adlt(i+1) Yng(i)-Adlt(i)') DO 370 I = 1,LM2 WRITE(7,380)YR(I),COVAA(I),COVYA1(I),COVYA(I) IF(MOD(YR(I),5).EQ.0)WRITE(7,130) 370 CONTINUE 380 FORMAT(1X,I4,4(6X,F16.9)) 390 WRITE(7,400) 400 FORMAT(//' Covariances of N:') 410 FORMAT(14(1X,F8.2)) DO 420 I = 1,LAST 420 WRITE(7,410) (ACOV(I,J),J=1,LM1) ENDIF CALL LRTEST(TITLE,S,SY,R,RY,LAST,YR) CALL GOF2(TITLE,LAST,YR,CHIA,IDFA) RETURN END