C*********************************************************************** C THIS SUBROUTINE IS DESIGNED TO SERVE TWO PURPOSES: 1) SUMMARIZE C THE ESTIMATION RESULTS FOR PERPENDICULAR DISTANCE ESTIMATORS AND C SIGHTING DISTANCE-ANGLE ESTIMATORS FOR ONE ANALYSIS AND 2) OB- C TAIN AN AVERAGE DENSITY ESTIMATE FOR ALL OF THE REPLICATE LINES C FOR EACH ESTIMATOR AND SUMMARIZE THEM. C C FUNCTIONS CALLED: TVALUE C*********************************************************************** SUBROUTINE SUMARY (ITIME) C*********************************************************************** C DECLARATIONS C*********************************************************************** INCLUDE 'PARMTR.INC' REAL LOGL, LCL LOGICAL GRP, POOL, PEST, SEST, DESC, DEF, CUTP, TRUNC INTEGER STATUS, PDEST, SDEST, CNT CHARACTER*1 DASH CHARACTER*24 ESTNAM(8) LOGICAL STRT, SET, WARN, HELP DIMENSION DNSTY(13,MAXLIN) C*********************************************************************** C COMMON STATEMENTS C*********************************************************************** COMMON /OPTION/ GRP, POOL, PEST, SEST, DESC, DEF, CUTP, TRUNC, 1 HELP COMMON /INTER/ KCUT, CUT(MAXCEL), FREQ(MAXCEL), NCUT, NKC(5), 1 NK(5), RFREQ(5,MAXCEL), CCUT(5,MAXCEL) COMMON /LINE/ SXL COMMON /STORE/ D(13), DL(13), DU(13), STD(13), LOGL(10), PB(5,10) COMMON /NUM/ XL(MAXLIN), WIDTH, N, CNT, CONV(3), VARN, IDF, WARN COMMON /PDOPT/ STRT, SET, PSTRT(10,MAXPAR), NPSET(10) COMMON /ESTM/ PDEST(10), SDEST(3), NSD, NPD COMMON /IND/ IC, II, IREP, STATUS C*********************************************************************** C DATA STATEMENTS C*********************************************************************** DATA DASH /'-'/ DATA ESTNAM/'Fourier Series','Exponential Power Series', 1'Exponential Polynomial','Negative Exponential','Half-normal', 2'Generalized Hayne','Modified Hayne','Hayne'/ C*********************************************************************** C FIRST STORE DENSITY ESTIMATES FOR THIS REPLICATION C*********************************************************************** IF (ITIME.EQ.2) GO TO 60 DO 10 JJ=1,13 10 DNSTY(JJ,II)=D(JJ) C*********************************************************************** C NEXT SUMMARIZE ESTIMATION RESULTS FOR THIS REPLICATION C*********************************************************************** IF (NPD+NSD.LE.1) RETURN CALL HEADER (1) WRITE (6,120) CALL HEADER (2) WRITE (6,130) WRITE (6,140) DO 40 I=1,NPD+NSD IF (I.LE.NPD) THEN IN=PDEST(I) ELSE IN=I-NPD IN=SDEST(IN)+5 ENDIF CV=0.0 IF (D(I).NE.0.0) CV=100.*STD(I)/D(I) WRITE (6,230) ESTNAM(IN),NPSET(I),D(I),STD(I),CV,DL(I),DU(I) 40 CONTINUE WRITE (6,140) C*********************************************************************** C FOR THE PERPENDICULAR DISTANCE ESTIMATORS SUMMARIZE THE LOG- C LIKELIHOOD VALUES AND CHI-SQUARE VALUES. C*********************************************************************** IF (.NOT.PEST) RETURN WRITE (6,240) NDASH=MAX0(21,NCUT*7) WRITE (6,250) (DASH,I=1,NDASH) DO 50 I=1,NPD IF (LOGL(I).EQ.0.) THEN WRITE (6,320) ESTNAM(PDEST(I)),NPSET(I),(PB(J,I),J=1,NCUT) ELSE WRITE (6,270) ESTNAM(PDEST(I)),NPSET(I),LOGL(I), 1 (PB(J,I),J=1,NCUT) ENDIF 50 CONTINUE WRITE (6,250) (DASH,I=1,NDASH) RETURN C*********************************************************************** C IF THE NPOL OPTION WAS USED THEN A CALL IS MADE TO HERE IF ALL C OF THE REPS WERE PROCESSED. AFTER ESTIMATES HAVE BEEN CALCULATED C FOR EACH REP THEN AN AVERAGE ESTIMATE OF DENSITY IS CALCULATED. C THOSE REPS IN WHICH NO OBSERVATIONS WERE MADE ARE INCLUDED BY C TREATING THEM AS ESTIMATING A ZERO DENSITY. C*********************************************************************** 60 CALL HEADER (1) WRITE (6,300) WRITE (6,280) WRITE (6,130) WRITE (6,140) DO 110 I=1,NPD+NSD VARD=0.0 DBAR=0.0 SX=0.0 DO 70 J=1,IREP SX=SX+XL(J) 70 DBAR=DBAR+DNSTY(I,J)*XL(J) DBAR=DBAR/SXL DO 80 J=1,IREP 80 VARD=VARD+XL(J)*((DNSTY(I,J)-DBAR)**2) VARD=VARD+(SXL-SX)*(DBAR**2) VARD=VARD/(SXL*FLOAT(IDF)) VARD=SQRT(VARD) CV=100.*VARD/DBAR XMULT=TVALUE(IDF) LCL=DBAR-XMULT*VARD UCL=DBAR+XMULT*VARD IF (I.LE.NPD) THEN IN=PDEST(I) ELSE IN=I-NPD IN=SDEST(IN)+5 ENDIF WRITE (6,230) ESTNAM(IN),NPSET(I),DBAR,VARD,CV,LCL,UCL 110 CONTINUE WRITE (6,140) WRITE (6,290) IDF,XMULT RETURN C*********************************************************************** C FORMAT STATEMENTS C*********************************************************************** C 120 FORMAT (//10X,60('*'),/10X,'*',58X,'*'/10X,'*',15X, 1 'Summary of Estimation Results',14X,'*'/10X,'*',58X,'*'/ 2 10X,60('*')) 130 FORMAT (///10X,'Estimator No. Density ', 1'Standard % Coef. 95 Percent' 2 /10X,' Name Par. Estimate ', 3' Error Var. Confidence Interval') 140 FORMAT (1X,24('-'),1X,4('-'),1X,10('-'),1X,10('-'),1X,7('-'), 1 1X,19('-')) 230 FORMAT (1X,A24,2X,I2,2X,G10.4,1X,G10.4,1X,F6.2,1X,G10.4,1X,G9.4) 240 FORMAT (///10X,'Table of chi-square probability levels (after', 1' pooling'/ 2 10X,'if any) and the log-likelihood value (if computable).'/ 2 /8X,'Estimator No. Log-lik. Chi-square prob. for' 4/8X,' Name Par. Value Each Cut Point Set') 250 FORMAT (1X,24('-'),1X,4('-'),1X,10('-'),1X,35A1) 270 FORMAT (1X,A24,2X,I2,2X,G10.4,1X,5(F5.3,2X)) 320 FORMAT (1X,A24,2X,I2,3X,'*****',5X,5(F5.3,2X)) 280 FORMAT (///' This analysis is done to provide an alternative', 1' method of calculating'/ 2' a variance for density. An average estimate of density is', 3' calculated by'/ 4' weighting each estimate by its line length. The average', 5' estimate and the'/ 6' individual estimates can then be used to calculate a variance', 7' for density.'/ 8' The average estimate of density is given by:'/ 8' D = (L(1)*D(1) + L(2)*D(2) + ... + L(R)*D(R))/(L(1) + L(2)', 9' + ... + L(r)).'/ A' The variance formula is given in Burnham et al. (1980).', B' The user should make'/ C' sure solutions were found for each replication before using', D' this analysis.') 290 FORMAT (//' The confidence interval for density was constructed', 1' with a t distribution'/ 2' with the degrees of freedom equal to one less the number of', 3' replicate density'/ 4' estimates. The t-value with ',I2,' degrees of freedom is ', 5 F5.3,'.') 300 FORMAT (//10X,60('*'),/10X,'*',58X,'*',/10X,'*',7X, 1 'Average Density Estimation from Replications',7X,'*'/10X, 2 '*',58X,'*',/10X,60('*')) END