1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 PAGE: 1 0REQUESTED OPTIONS (EXECUTE): OPT(2) XREF SXM 0OPTIONS IN EFFECT: NOLIST NOMAP XREF GOSTMT NODECK SOURCE TERM OBJECT FIXED TRMFLG SRCFLG NOSYM NORENT SDUMP(ISN) SXM NOVECTOR IL(DIM) NOTEST SC(*) NODC NOEC NOEMODE NOICA NODIRECTIVE NODBCS NOSAA NOPARALLEL NOSAVE NOTABS OPT(2) LANGLVL(77) NOFIPS FLAG(I) AUTODBL(NONE) LINECOUNT(60) CHARLEN(500) 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C C CLUSPAC: Computer Programs for Mixture-Model Clustering C C C C COPYRIGHT (C) 1991, 1992 STANLEY L. SCLOVE. C C C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C$STATEMENTS=100000 C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C C CMS DSN = MIX1DT CLUSPAC C C C C "MIX1DT CLUSPAC" IS A PROGRAM FOR CLUSTERING UNIVARIATE DATA C C (DATA ON THE LINE) BY ITERATIVE MAXIMIZATION OF THE MIXTURE- C C MODEL LIKELIHOOD C C C C N K C C --- -- C C L = | | > P(C)*F(X(I)|C) C C | | -- C C I=1 C=1 C C C C C C REFERENCE: C C C C Wolfe, J. H. (1970). Pattern clustering by multivariate C C mixture analysis. Multivariate Behavioral Research 5, 329-350. C C C C MANUAL MODE: NUMBER OF CLUSTERS AND INITIAL MEANS ARE C C INPUT. (USE PROGRAM MIX1DTA FOR AUTOMATIC SETTING OF C C NUMBERS OF CLUSTERS AND INITIAL MEANS.) C C C C C C PROGRAMMED BY C C DR. STANLEY L. SCLOVE 312/996-2681 C C DEPARTMENT OF C C INFORMATION & DECISION SCIENCES M/C 294 C C COLLEGE OF BUSINESS ADMINISTRATION C C UNIVERSITY OF ILLINOIS AT CHICAGO C C BOX 4348 C C CHICAGO, IL 60680-4348 C C C C C C VERSION 1.6 15-MAY-92 C C C 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 2 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C C C C C C C RESTRICTIONS (CAN BE MODIFIED): C C N, SAMPLE SIZE, AT MOST 999; C C K, NUMBER OF CLUSTERS, AT MOST 29; C C ITER, MAXIMUM NUMBER OF ITERATIONS, 99. C C C C C C C C CONTROL CARDS: C C C C (1) DATASET TITLE C C (2) N, IN FORMAT (2X,I4) C C (3) FMT, IN FORMAT (18A4), E.G., (1X,F4.1). C C ALLOW AT LEAST ONE BLANK IN FMT: IT WILL ALSO BE USED C C FOR OUTPUT, WHERE CC1 IS FOR CARRIAGE CONTROL. C C ALLOW A CC FOR THE DECIMAL POINT ON OUTPUT, C C WHETHER OR NOT THERE IS ONE ON INPUT. C C (4) DATA, IN FORMAT SPECIFIED BY FMT C C (5) K, NUMBER OF CLUSTERS, IN FORMAT (2X,I1) C C (6) K INITIAL VALUES OF PRIOR PROBS, MEANS AND VARIANCES C C IN FORMAT (5X,F3.2,2X,F8.2,2X,F8.2). C C C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C 1 DIMENSION X(999),XMNDSQ(999),ICLUS(999),IOTA(999) 2 DIMENSION DSQ(29),C(29),SUM(29) 3 DIMENSION TITLE(18) 4 DIMENSION B(29),NC(29),XMEAN(29),XMEANOL(29) 5 DIMENSION FMT(18) 6 DIMENSION SS(29),SSD(29) 7 DIMENSION SD(29) 8 DIMENSION VAR(29) 9 DIMENSION ICLSOL(999) 10 DIMENSION F(999,29) 11 DIMENSION P(29),XNC(29) 12 DIMENSION PP(29,999) 13 DIMENSION XMXPR(999) 14 DIMENSION DENOM(999) C 15 DOUBLE PRECISION SUM,SS,F,P,PP,SUMSQS,TOTAL,SSD,XNC,VAR 16 DOUBLE PRECISION TEMP 17 DOUBLE PRECISION XMEAN,XMEANOL,XMN2LL,XMN2LLOL,WGSS,SSDEVS C 18 READ (5,64000) TITLE C C WRITE PROGRAM INFORMATION. 19 WRITE (6,56000) 20 WRITE (6,32000) 21 WRITE (6,34000) C 22 WRITE (6,70000) TITLE C 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 3 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C READ SAMPLE SIZE, N. 23 READ (5,12000) N 24 XN = N 25 WRITE (6,66000) N C C READ DATA FORMAT. 26 READ (5,64000) FMT C C READ DATA AND C COMPUTE STATISTICS OF WHOLE SAMPLE: 27 TOTAL=0.0 28 SUMSQS=0.0 29 SSDEVS=0.0 C 30 DO 300 I = 1,N 1 31 READ (5,FMT) X(I) 1 32 TOTAL = TOTAL + X(I) 1 33 SUMSQS = SUMSQS + X(I)*X(I) 1 34 IF (I .EQ. 1) GO TO 100 1 35 GO TO 200 1 36 100 XMAX = X(1) 1 37 XMIN = X(1) 1 38 200 CONTINUE 1 39 IF (X(I) .LT. XMIN) XMIN=X(I) 1 41 IF (X(I) .GT. XMAX) XMAX=X(I) 1 43 300 CONTINUE C WRITE DATA: 44 WRITE (6,24000) 45 IF ( N .GE. 31) GO TO 350 46 WRITE (6,FMT) (X(I),I=1,N) 47 350 CONTINUE C WRITE SUMMARY STATISTICS FOR WHOLE SAMPLE: 48 WRITE (6,44000) 49 WRITE (6,FMT) XMIN 50 WRITE (6,46000) 51 WRITE (6,FMT) XMAX 52 XBAR = TOTAL/XN 53 SSDEVS = SUMSQS- TOTAL*TOTAL/XN 54 VARHAT = SSDEVS/XN 55 WRITE (6,50000) XBAR 56 WRITE (6,52000) VARHAT C 57 PI = 3.141593 58 TEMP = 2.0*PI*SSDEVS/N 59 XMN2LL = N*(1.0 + DLOG(TEMP)) 60 STDDEV = SQRT(VARHAT) 61 WRITE (6,54000) SSDEVS, XMN2LL, STDDEV 62 NOPARM = 1 + 1 63 AIC = XMN2LL + 2.0*NOPARM 64 WRITE (6,82000) AIC 65 SCH = XMN2LL + ALOG(XN)*NOPARM 66 XKASH = SCH - ALOG(2*VARHAT**3) 67 WRITE (6,48000) XKASH C C C READ K, NUMBER OF CLUSTERS. 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 4 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 68 READ (5,10000) K 69 WRITE (6,36000) K C C READ INITIAL PRIOR PROBS, MEANS AND VARIANCES: 70 DO 400 IC=1,K 1 71 READ (5,14000) P(IC),XMEAN(IC),VAR(IC) 1 72 400 CONTINUE C C C C WRITE INITIAL PRIOR PROBS, MEANS AND VARIANCES: 73 WRITE (6,26000) 74 DO 500 IC = 1,K 1 75 WRITE (6,28000) IC, P(IC), XMEAN(IC), VAR(IC) 1 76 500 CONTINUE C C SET CONSTANTS. 77 DO 600 INTGER=1,N 1 78 IOTA(INTGER) = INTGER 1 79 600 CONTINUE C 80 ITER = 1 C IF THE INITIAL PRIOR PROBABILITIES ARE EQUAL, THEN THE C FIRST ITERATION IS EQUIVALENT TO MINIMUM DISTANCE CLUSTERING C TO INITIAL MEANS (I.E., "ISODATA"). C IN GENERAL, THE CLUSTERING IS BY MAXIMUM POSTERIOR C PROBABILITY CLUSTERING. C 81 700 CONTINUE C 82 IF (ITER .EQ. 1) GO TO 900 C C STORE OLD CLUSTERING RESULTS FOR TESTING CONVERGENCE: 83 XMN2LLOL = XMN2LL 84 DO 800 I = 1,N 1 85 ICLSOL(I) = ICLUS(I) 1 86 800 CONTINUE 87 DO 850 IC = 1,K 1 88 XMEANOL(IC) = XMEAN(IC) 1 89 850 CONTINUE C C COMMENCE DISTANCE COMPUTATIONS. 90 900 CONTINUE 91 DO 1100 I = 1,N 1 92 DO 1100 IC = 1,K 2 93 DSQ(IC) = ( XMEAN(IC) - X(I) )**2 2 94 ZSQ = DSQ(IC)/VAR(IC) 2 95 IF ( ZSQ .LE. 174.673 ) GO TO 1000 2 96 F(I,IC) = 0.0 2 97 GO TO 1100 2 98 1000 CONTINUE C C NOTE THAT A PROB. DENSITY FUNCTION OTHER THAN THE GAUSSIAN C COULD BE USED HERE: 2 99 F(I,IC) =(EXP(-ZSQ/2.0))/(SQRT(2*PI*VAR(IC))) 2 100 1100 CONTINUE 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 5 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C XMNDSQ(I) = MIN SQ. DISTANCE FROM X(I) TO ANY MEAN C C C 101 1200 CONTINUE C C COMPUTE POSTERIOR PROBABILITIES OF GROUP MEMBERSHIP: 102 DO 1300 I = 1,N 1 103 DENOM(I) = 0.0 1 104 DO 1300 IH=1,K 2 105 DENOM(I) = DENOM(I) + P(IH)*F(I,IH) 2 106 1300 CONTINUE 107 DO 1400 I = 1,N 1 108 DO 1400 IC=1,K C IF ( DENOM(I) .EQ. 0.0 ) DENOM(I)=0.0001 2 109 PP(IC,I)= P(IC)*F(I,IC)/DENOM(I) 2 110 1400 CONTINUE C C COMPUTE NEW LABELS BY MAX POSTERIOR PROBABILITY: 111 DO 1600 I = 1,N 1 112 XMXPR(I) = PP(1,I) 1 113 ICLUS(I) = 1 1 114 DO 1600 IC = 2,K 2 115 IF ( PP(IC,I) .GT. XMXPR(I) ) GO TO 1500 2 116 GO TO 1600 2 117 1500 XMXPR(I) = PP(IC,I) 2 118 ICLUS(I) = IC 2 119 1600 CONTINUE C 120 IF (N .GE. 31) GO TO 1650 C C WRITE NEW LABELS: 121 WRITE (6,16000) 122 WRITE (6,18000) (IOTA(I), I=1,N) 123 WRITE (6,20000) (ICLUS(I), I=1,N) 124 1650 CONTINUE C C UPDATE CLUSTER PRIOR PROBABILITIES P(IC), MEANS XMEAN(IC) AND C VARIANCES VAR(IC): 125 WGSS = 0.0 126 DO 1800 IC = 1,K 1 127 XNC(IC) = 0.0 C XNC(IC) WILL BE THE SUM OVER ALL N OBSERVATIONS OF THEIR C POSTERIOR PROBABILITIES OF MEMBERSHIP IN CLUSTER IC. 1 128 SUM(IC) = 0.0 1 129 SS(IC) = 0.0 1 130 DO 1700 I = 1,N 2 131 XNC(IC) = XNC(IC) + PP(IC,I) 2 132 SUM(IC) = SUM(IC) + PP(IC,I)*X(I) 2 133 SS(IC) = SS(IC) + PP(IC,I)*X(I)*X(I) 2 134 1700 CONTINUE 1 135 IF ( XNC(IC) .EQ. 0.0 ) GO TO 2700 1 136 XMEAN(IC) = SUM(IC)/XNC(IC) 1 137 SSD(IC) = SS(IC) - SUM(IC)*SUM(IC)/XNC(IC) 1 138 VAR(IC) = SSD(IC)/(XNC(IC)) C IF ( VAR(IC) .LE. 0.0 ) VAR(IC) = 0.0001 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 6 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 1 139 SD(IC) = SQRT(VAR(IC)) 1 140 P(IC) = XNC(IC)/XN 1 141 WGSS = WGSS + SSD(IC) 1 142 1800 CONTINUE C C COUNT NUMBERS IN CLUSTERS: 143 DO 1900 IC = 1,K 1 144 NC(IC) = 0 1 145 1900 CONTINUE 146 DO 2000 I = 1,N 1 147 IGROUP = ICLUS(I) 1 148 NC(IGROUP) = NC(IGROUP) + 1 1 149 2000 CONTINUE C C C 150 VARHAT = WGSS/XN C 151 SUMLNF = 0.0 152 DO 2200 I=1,N 1 153 SUMPXF = 0.0 1 154 DO 2100 IC=1,K 2 155 SUMPXF = SUMPXF + P(IC)*F(I,IC) 2 156 2100 CONTINUE 1 157 SUMLNF = SUMLNF + ALOG(SUMPXF) 1 158 2200 CONTINUE C 159 XMN2LL = -2.0*SUMLNF C C 160 WGMS = WGSS/(N-K) 161 WRITE (6,40000) WGSS, XMN2LL, WGMS 162 STDERR = SQRT(WGMS) 163 WRITE (6,74000) STDERR C 164 KM1 = K-1 165 DO 2300 IC=1,KM1 1 166 ICP1 = IC+1 C B(IC) IS BOUNDARY BETWEEN G-TH AND G+1-ST CLASSES. 1 167 B(IC) = ( XMEAN(IC) + XMEAN(ICP1) )/2.0 1 168 BDYADJ = DLOG(P(ICP1)) - DLOG(P(IC)) 1 169 BDYADJ = BDYADJ/(XMEAN(ICP1) - XMEAN(IC)) 1 170 BDYADJ = VARHAT*BDYADJ 1 171 B(IC) = B(IC) - BDYADJ 1 172 2300 CONTINUE 173 WRITE (6,58000) ITER 174 WRITE (6,42000) (B(IC),IC=1,KM1) 175 WRITE (6,38000) (XMEAN(IC),IC=1,K) 176 IF (ITER .EQ. 1) GO TO 2500 177 DO 2400 I = 1,N C If any case changes cluster, another iteration will be performed. 1 178 IF (ICLUS(I) .EQ. ICLSOL(I)) GO TO 2400 1 179 GO TO 2500 1 180 2400 CONTINUE 181 WRITE (6,79000) C Even if no case changes cluster, another iteration will be 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 7 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C performed if the parameter estimates have changed much. 182 DO 2450 IC = 1,K 1 183 TEST=(XMEAN(IC)-XMEANOL(IC))/XMEANOL(IC) 1 184 TOLERANS = 0.01 1 185 IF (ABS(TEST) .LE. TOLERANS) GO TO 2450 1 186 GO TO 2500 1 187 2450 CONTINUE 188 GO TO 2900 189 2500 CONTINUE 190 ITER = ITER + 1 191 IF (ITER.GE.100) GO TO 2600 192 GO TO 700 193 2600 WRITE (6,80000) 194 2700 CONTINUE 195 WRITE (6,88000) C C Next the results will be written: C 196 2900 CONTINUE C 197 WRITE (6,60000) ( P(IC), IC = 1,K ) 198 WRITE (6,62000) (NC(IC),IC=1,K) 199 WRITE (6,72000) (VAR(IC),IC=1,K) 200 WRITE (6,68000) (SD(IC), IC=1,K) C VARHAT IS MLE OF VARIANCE. 201 VARHAT = WGSS/N 202 WRITE (6,30000) VARHAT C C 203 IF (N .GE. 31) GO TO 3050 204 WRITE (6,78000) 205 DO 3000 I=1,N 1 206 WRITE (6,76000) X(I),ICLUS(I),(PP(IC,I),IC=1,K) 1 207 3000 CONTINUE C 208 3050 CONTINUE C C COMPUTE MODEL-SELECTION CRITERIA: C NO. PARAMETERS = K MEANS + K VARIANCES + (K-1) PROBS. 209 NOPARM = K + K + (K-1) C 210 WRITE (6,86000) NOPARM 211 AIC = XMN2LL + 2.0*NOPARM 212 SCH = XMN2LL + ALOG(XN)*NOPARM C C SCHWARZ' CRITERION IS FIRST-DEGREE EXPANSION OF C LOG POSTERIOR PROBABILITY OF THE MODEL. C KASHYAP'S CRITERION IS SECOND-DEGREE EXPANSION OF SAME. 213 XKASH = SCH - ALOG(2*VARHAT**3) C 214 WRITE (6,82000) AIC 215 WRITE (6,84000) SCH 216 WRITE (6,48000) XKASH C C C 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 8 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 217 STOP 218 10000 FORMAT(2X,I1) 219 12000 FORMAT(2X,I4) 220 14000 FORMAT(5X,F3.2,2X,F8.2,2X,F8.2) 221 16000 FORMAT(1X,'CLUSTERING') 222 18000 FORMAT(1X,'POINT: '/, (1X,40I3)) 223 20000 FORMAT(1X,'CLUSTER: '/, (1X,40I3)) 224 24000 FORMAT(1X,'DATA:'/) 225 26000 FORMAT(/1X,'INITIAL PRIOR PROBS, MEANS AND VARIANCES:') 226 28000 FORMAT(1X,I2,2X,3F11.2) 227 30000 FORMAT(1X, 'M.L. ESTIMATE OF COMMON VARIANCE = ',F14.5/) 228 32000 FORMAT(1X,'CMS DSN = MIX1DT CLUSPAC') 229 34000 FORMAT(1X,'COPYRIGHT (C) 1991, 1992 STANLEY L. SCLOVE'/) 230 36000 FORMAT('1',1X,'K = ',I1,' CLUSTERS') 231 38000 FORMAT(1X,'MEANS: ',9F13.3) 232 40000 FORMAT(/1X,'WGSS = ',F14.4,' MINUS 2 LOG LIKELIHOOD = ', XF14.4, ' WGMS = ',F14.4/) 233 42000 FORMAT(1X,'BOUNDARIES:', 8X, 9F13.3) 234 44000 FORMAT(1X,'MINIMUM OF SAMPLE: ') 235 46000 FORMAT(1X,'MAXIMUM OF SAMPLE: ') 236 48000 FORMAT( 1X,'KASHYAP CRITERION = ', F14.4) 237 50000 FORMAT(/,1X,'MEAN = ', F14.4) 238 52000 FORMAT(1X, 'M.L. ESTIMATE OF VARIANCE = ',F14.5/) 239 54000 FORMAT(/1X,'SSDEVS = ',F14.4,' MINUS 2 LOG LIKELIHOOD = ', XF14.4, ' STDDEV = ',F14.4/) 240 56000 FORMAT('1','.................................................', X/,1X,'PROGRAM MIX1DT CLUSPAC'/1X,'MIXTURE MODEL CLUSTERING'/ X,1X,'FOR UNIVARIATE DATA (DATA ON THE LINE)'/ X,1X,'WITH UNEQUAL CLASS VARIANCES '/ X1X,'DEVELOPED AND PROGRAMMED BY DR. STANLEY L. SCLOVE' X/,1X,'VERSION 1.6 15-MAY-92 '/) 241 58000 FORMAT(1X,'ITERATION ', I2) 242 60000 FORMAT(1X,'PROBS:',5X,9F13.3) 243 62000 FORMAT(1X,'NUMBERS:',3X,9(I10,3X)) 244 64000 FORMAT(18A4) 245 66000 FORMAT(1X,'N = ',I3/) 246 68000 FORMAT(1X,'STD.DEVS.: ',9F13.3) 247 70000 FORMAT(1X,18A4/) 248 72000 FORMAT(1X,'VARIANCES: ',9F13.3) 249 74000 FORMAT(1X,'STD.ERROR=SQRT(WGMS) = ',F13.4/) 250 76000 FORMAT(1X,F8.2,2X,I2,3X,(8F6.3/)) 251 78000 FORMAT(1X/,1X,'DATA, LABELS AND PROBS OF CLUSTER MEMBERSHIP'/) 252 79000 FORMAT(/,1X,'NO CASE CHANGED CLUSTERS IN THIS ITERATION.') 253 80000 FORMAT(1X,'PROGRAM HAS NOT CONVERGED IN 99 ITERATIONS. STOP') 254 82000 FORMAT(1X,'AIC = ', F14.4 ) 255 84000 FORMAT(1X,'SCHWARZ CRITERION = ', F14.4 ) 256 86000 FORMAT(/1X,'NUMBER OF PARAMETERS = ',I4 ) 257 88000 FORMAT(1X,'NO OBSERVATIONS IN GROUP ',I3,'. STOP') 258 END 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 9 0SYMBOL CROSS REFERENCE DICTIONARY 0PROGRAM NAME: MAIN# 0TAGS: A-ARRAY I-INTRINSIC FUNCTION S-ASSIGNED C-COMMON K-NAMED CONSTANT T-EXPLICITLY TYPED D-DUMMY ARGUMENT N-ENTRY V-INITIAL VALUE E-EQUIVALENCED P-PROMOTED X-EXTERNAL SUBPROGRAM F-STATEMENT FUNCTION Q-PADDED Y-DYNAMIC COMMON G-GENERIC NAME R-SUBPROGRAM NAME Z-EXTENDED COMMON 0NAME TYPE TAG DECLARED REFS (F:REFD S:SET B:REFD/MAY BE SET) +________ ____ ______ ________ ______________________________________ 0ABS GI 185 AIC R*4 63S 64F 211S 214F ALOG R*4 I 65 66 157 212 213 B R*4 A 4 167S 171F 171S 174F BDYADJ R*4 168S 169F 169S 170F 170S 171F C R*4 A 2 UNREFERENCED DENOM R*4 A 14 103S 105F 105S 109F DLOG R*8 I 59 168 168 DSQ R*4 A 2 93S 94F EXP GI 99 F R*8 AT 10 15 96S 99S 105F 109F 155F FMT R*4 A 5 26S 31 46F 49F 51F I I*4 30S 31 32F 33F 33F 34F 39F 40F 41F 42B 46S 46S 84S 85F 85B 91S 93F 96F 99B 102S 103F 105F 105F 105B 107S 109F 109F 109B 111S 112F 112F 113F 115F 115F 117F 117F 118B 122S 122S 123S 123S 130S 131F 132F 132F 133F 133F 133B 146S 147B 152S 155B 177S 178F 178B 205S 206F 206F 206B IC I*4 70S 71F 71S 71F 74S 75F 75F 75F 75B 87S 88F 88B 92S 93F 93F 94F 94F 96F 99F 99B 108S 109F 109F 109B 114S 115F 117F 118B 126S 127F 128F 129F 131F 131F 131F 132F 132F 132F 133F 133F 133F 135F 136F 136F 136F 137F 137F 137F 137F 137F 138F 138F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 10 0TAGS: A-ARRAY I-INTRINSIC FUNCTION S-ASSIGNED C-COMMON K-NAMED CONSTANT T-EXPLICITLY TYPED D-DUMMY ARGUMENT N-ENTRY V-INITIAL VALUE E-EQUIVALENCED P-PROMOTED X-EXTERNAL SUBPROGRAM F-STATEMENT FUNCTION Q-PADDED Y-DYNAMIC COMMON G-GENERIC NAME R-SUBPROGRAM NAME Z-EXTENDED COMMON 0NAME TYPE TAG DECLARED REFS (F:REFD S:SET B:REFD/MAY BE SET) +________ ____ ______ ________ ______________________________________ 0 138F 139F 139F 140F 140F 141B 143S 144B 154S 155F 155B 165S 166F 167F 167F 168F 169F 171F 171B 174S 174S 175S 175S 182S 183F 183F 183B 197S 197S 198S 198S 199S 199S 200S 200S 206S 206S ICLSOL I*4 A 9 85S 178F ICLUS I*4 A 1 85F 113S 118S 123F 147F 178F 206F ICP1 I*4 166S 167F 168F 169F IGROUP I*4 147S 148F 148F IH I*4 104S 105F 105B INTGER I*4 77S 78F 78B IOTA I*4 A 1 78S 122F ITER I*4 80S 82F 173F 176F 190F 190S 191F K I*4 68S 69F 70F 74F 87F 92F 104F 108F 114F 126F 143F 154F 160F 164F 175F 182F 197F 198F 199F 200F 206F 209F 209F 209F KM1 I*4 164S 165F 174F N I*4 23S 24F 25F 30F 45F 46F 58F 59F 77F 84F 91F 102F 107F 111F 120F 122F 123F 130F 146F 152F 160F 177F 201F 203F 205F NC I*4 A 4 144S 148F 148S 198F NOPARM I*4 62S 63F 65F 209S 210F 211F 212F P R*8 AT 11 15 71S 75F 105F 109F 140S 155F 168F 168F 197F PI R*4 57S 58F 99F PP R*8 AT 12 15 109S 112F 115F 117F 131F 132F 133F 206F SCH R*4 65S 66F 212S 213F 215F SD R*4 A 7 139S 200F SQRT GI 60 99 139 162 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 11 0TAGS: A-ARRAY I-INTRINSIC FUNCTION S-ASSIGNED C-COMMON K-NAMED CONSTANT T-EXPLICITLY TYPED D-DUMMY ARGUMENT N-ENTRY V-INITIAL VALUE E-EQUIVALENCED P-PROMOTED X-EXTERNAL SUBPROGRAM F-STATEMENT FUNCTION Q-PADDED Y-DYNAMIC COMMON G-GENERIC NAME R-SUBPROGRAM NAME Z-EXTENDED COMMON 0NAME TYPE TAG DECLARED REFS (F:REFD S:SET B:REFD/MAY BE SET) +________ ____ ______ ________ ______________________________________ 0SS R*8 AT 6 15 129S 133F 133S 137F SSD R*8 AT 6 15 137S 138F 141F SSDEVS R*8 T 17 29S 53S 54F 58F 61F STDDEV R*4 60S 61F STDERR R*4 162S 163F SUM R*8 AT 2 15 128S 132F 132S 136F 137F 137F SUMLNF R*4 151S 157F 157S 159F SUMPXF R*4 153S 155F 155S 157F SUMSQS R*8 T 15 28S 33F 33S 53F TEMP R*8 T 16 58S 59F TEST R*4 183S 185F TITLE R*4 A 3 18S 22F TOLERANS R*4 184S 185F TOTAL R*8 T 15 27S 32F 32S 52F 53F 53F VAR R*8 AT 8 15 71S 75F 94F 99F 138S 139F 199F VARHAT R*4 54S 56F 60F 66F 150S 170F 201S 202F 213F WGMS R*4 160S 161F 162F WGSS R*8 T 17 125S 141F 141S 150F 160F 161F 201F X R*4 A 1 31S 32F 33F 33F 36F 37F 39F 40F 41F 42F 46F 93F 132F 133F 133F 206F XBAR R*4 52S 55F XKASH R*4 66S 67F 213S 216F XMAX R*4 36S 41F 42S 51F XMEAN R*8 AT 4 17 71S 75F 88F 93F 136S 167F 167F 169F 169F 175F 183F XMEANOL R*8 AT 4 17 88S 183F 183F XMIN R*4 37S 39F 40S 49F XMNDSQ R*4 A 1 UNREFERENCED XMN2LL R*8 T 17 59S 61F 63F 65F 83F 159S 161F 211F 212F XMN2LLOL R*8 T 17 83S XMXPR R*4 A 13 112S 115F 117S XN R*4 24S 52F 53F 54F 65F 140F 150F 212F XNC R*8 AT 11 15 127S 131F 131S 135F 136F 137F 138F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 12 0TAGS: A-ARRAY I-INTRINSIC FUNCTION S-ASSIGNED C-COMMON K-NAMED CONSTANT T-EXPLICITLY TYPED D-DUMMY ARGUMENT N-ENTRY V-INITIAL VALUE E-EQUIVALENCED P-PROMOTED X-EXTERNAL SUBPROGRAM F-STATEMENT FUNCTION Q-PADDED Y-DYNAMIC COMMON G-GENERIC NAME R-SUBPROGRAM NAME Z-EXTENDED COMMON 0NAME TYPE TAG DECLARED REFS (F:REFD S:SET B:REFD/MAY BE SET) +________ ____ ______ ________ ______________________________________ 0 140F ZSQ R*4 94S 95F 99F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 13 0LABEL CROSS REFERENCE DICTIONARY 0TAGS: A-USED AS ARGUMENT F-FORMAT S-USED IN ASSIGN STATEMENT B-OBJECT OF BRANCH N-NON-EXECUTABLE 0 LABEL TAG DEFINED REFERENCED +____________ ___ _______ ___________________________________________ 0 100 B 36 34 200 B 38 35 300 43 30 350 B 47 45 400 72 70 500 76 74 600 79 77 700 B 81 192 800 86 84 850 89 87 900 B 90 82 1000 B 98 95 1100 B 100 91 92 97 1200 101 UNREFERENCED 1300 106 102 104 1400 110 107 108 1500 B 117 115 1600 B 119 111 114 116 1650 B 124 120 1700 134 130 1800 142 126 1900 145 143 2000 149 146 2100 156 154 2200 158 152 2300 172 165 2400 B 180 177 178 2450 B 187 182 185 2500 B 189 176 179 186 2600 B 193 191 2700 B 194 135 2900 B 196 188 3000 207 205 3050 B 208 203 10000 NF 218 68 12000 NF 219 23 14000 NF 220 71 16000 NF 221 121 18000 NF 222 122 20000 NF 223 123 24000 NF 224 44 26000 NF 225 73 28000 NF 226 75 30000 NF 227 202 32000 NF 228 20 34000 NF 229 21 36000 NF 230 69 38000 NF 231 175 40000 NF 232 161 42000 NF 233 174 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN MAY 15, 1992 23:40:56 NAME:MAIN# PAGE: 14 0TAGS: A-USED AS ARGUMENT F-FORMAT S-USED IN ASSIGN STATEMENT B-OBJECT OF BRANCH N-NON-EXECUTABLE 0 LABEL TAG DEFINED REFERENCED +____________ ___ _______ ___________________________________________ 0 44000 NF 234 48 46000 NF 235 50 48000 NF 236 67 216 50000 NF 237 55 52000 NF 238 56 54000 NF 239 61 56000 NF 240 19 58000 NF 241 173 60000 NF 242 197 62000 NF 243 198 64000 NF 244 18 26 66000 NF 245 25 68000 NF 246 200 70000 NF 247 22 72000 NF 248 199 74000 NF 249 163 76000 NF 250 206 78000 NF 251 204 79000 NF 252 181 80000 NF 253 193 82000 NF 254 64 214 84000 NF 255 215 86000 NF 256 210 88000 NF 257 195 0 NUMBER MODULE LEVEL ISN VS FORTRAN ERROR MESSAGES + ______ ______ _____ ___ __ _______ _____ ________ 0ILX3470I VPCQ 0(I) "XMAX,XMIN" MAY BE UNINITIALIZED WHEN USED 0*STATISTICS* SOURCE STATEMENTS: 256, PROGRAM SIZE: 499232 BYTES, PROGRAM NAME: MAIN#, PAGE: 1 *STATISTICS* 1 DIAGNOSTIC GENERATED. SEVERITY CODE IS 0. **MAIN#** END OF COMPILATION 1 ****** TIME STAMP: 92.13623.40.56