1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 PAGE: 1 0REQUESTED OPTIONS (EXECUTE): OPT(2) XREF SXM 0OPTIONS IN EFFECT: NOLIST NOMAP XREF NOGOSTMT 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 OPT(2) LANGLVL(77) NOFIPS FLAG(I) AUTODBL(NONE) LINECOUNT(60) CHARLEN(500) NAME(MAIN#) 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) 1998 Stanley Louis Sclove C C C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C$STATEMENTS=100000 C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C C Program MIX1BETA in CLUSPAC from ISOPAC C C C C ADAPTED FROM "MIX1DT CLUSPAC", C C 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 where the class-conditional probability density functions C C f(x|c) are betas. C C C C C The pdf is C r-1 s-1 C G(r + s) x (1-x) C f(x) = --------------------------------- r > 0, s > 0 C G(r) G(s) C C where G is the gamma function. C The beta function is B(r,s) = G(r)G(s)/G(r+s) . C so the constant in f(x) is 1/B(r,s) . C C C The mean is r/(r+s) = mu C C The variance is rs/(r+s)**2(r+s+1) = var . C C C C This gives r = c*mu, s = c*(1-mu), C where C c = [mu(1-mu)/var - 1] . C 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 2 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 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 601 S. MORGAN ST. C C CHICAGO, IL 60607-7124 C C C C C C VERSION 2.0 13-DEC-98 C C C C RESTRICTIONS (CAN BE MODIFIED): C C N, SAMPLE SIZE, AT MOST 999; C C K, NUMBER OF CLUSTERS, AT MOST 11; C C ITER, MAXIMUM NUMBER OF ITERATIONS, 4. C C (The program seems to wander off to an unwanted trivial, C global maximum as iterations go on.) C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C CONTROL STATEMENTS: 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(11),C(11),SUM(11),R(11),S(11),RPLUSS(11) 3 DIMENSION TITLE(18) 4 DIMENSION B(11),NC(11),XMEAN(11),XMEANOL(11),RECIPBET(11) 5 DIMENSION FMT(18) 6 DIMENSION SS(11),SSD(11) 7 DIMENSION SD(11) 8 DIMENSION VAR(11) 9 DIMENSION ICLSOL( 999) 10 DIMENSION F( 999,11) 11 DIMENSION P(11),XNC(11) 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 3 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 12 DIMENSION PP(11, 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,XLL,DENOM,SUMLNF,SUMPXF 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 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 C Initialize max and min: 30 XMAX = 0.0 31 XMIN = 1.0 32 DO 300 I = 1,N 1 33 READ (5,FMT) X(I) C Data cleaning: 1 34 IF ( X(I) .LE. 0.0001 ) X(I)=0.0001 1 36 IF ( X(I) .GE. 0.9999 ) X(I)=0.9999 C 1 38 TOTAL = TOTAL + X(I) 1 39 SUMSQS = SUMSQS + X(I)*X(I) 1 40 IF (X(I) .LT. XMIN) XMIN=X(I) 1 42 IF (X(I) .GT. XMAX) XMAX=X(I) 1 44 300 CONTINUE C WRITE DATA: 45 WRITE (6,24000) 46 IF ( N .GE. 31) GO TO 350 47 WRITE (6,FMT) (X(I),I=1,N) 48 350 CONTINUE C WRITE SUMMARY STATISTICS FOR WHOLE SAMPLE: 49 WRITE (6,44000) 50 WRITE (6,FMT) XMIN 51 WRITE (6,46000) 52 WRITE (6,FMT) XMAX 53 XBAR = TOTAL/XN 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 4 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 54 SSDEVS = SUMSQS- TOTAL*TOTAL/XN 55 VARHAT = SSDEVS/XN 56 WRITE (6,50000) XBAR 57 WRITE (6,52000) VARHAT C C Compute likelihood for a single beta. 58 XC = XBAR*(1-XBAR)/VARHAT - 1 59 XR = XBAR*XC 60 XS = (1-XBAR)*XC 61 XLL = 0.0 62 DO 105 I = 1,N 1 63 TEMP = (XR-1)*ALOG(X(I)) + (XS-1)*ALOG(1- X(I)) 1 64 XLL = XLL + TEMP 1 65 105 CONTINUE 66 XLL = XLL + ALOG( GAMMA(XR)*GAMMA(XS)/GAMMA(XR+XS) ) 67 XMN2LL = -2.0*XLL 68 STDDEV = SQRT(VARHAT) 69 WRITE (6,54000) SSDEVS, XMN2LL, STDDEV 70 NOPARM = 1 + 1 71 AIC = XMN2LL + 2.0*NOPARM 72 WRITE (6,82000) AIC 73 SCH = XMN2LL + ALOG(XN)*NOPARM 74 TEMP = 2*VARHAT**3 75 IF ( TEMP .LE. 0.0 ) TEMP = 0.0001 77 XKASH = SCH - DLOG(TEMP) 78 WRITE (6,48000) XKASH C C C READ K, NUMBER OF CLUSTERS. 79 READ (5,10000) K 80 WRITE (6,36000) K C C READ INITIAL PRIOR PROBS, MEANS AND VARIANCES: 81 DO 400 IC=1,K 1 82 READ (5,14000) P(IC),XMEAN(IC),VAR(IC) 1 83 RPLUSS(IC) = XMEAN(IC)*(1-XMEAN(IC))/VAR(IC) - 1 1 84 R(IC) = XMEAN(IC)*RPLUSS(IC) 1 85 S(IC) = (1-XMEAN(IC))*RPLUSS(IC) 1 86 400 CONTINUE C C C C WRITE INITIAL PRIOR PROBS, MEANS AND VARIANCES: 87 WRITE (6,26000) 88 DO 500 IC = 1,K 1 89 WRITE (6,28000) IC, P(IC), XMEAN(IC), VAR(IC) 1 90 WRITE (6,28100) IC, R(IC), S(IC) 1 91 500 CONTINUE C C SET CONSTANTS. 92 DO 600 INTGER=1,N 1 93 IOTA(INTGER) = INTGER 1 94 600 CONTINUE C 95 ITER = 1 C THE CLUSTERING IS BY MAXIMUM POSTERIOR 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 5 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C PROBABILITY CLUSTERING. C 96 700 CONTINUE C 97 IF (ITER .EQ. 1) GO TO 900 C because there are no old results to store. C C If it's not the first iteration, C STORE OLD CLUSTERING RESULTS FOR TESTING CONVERGENCE: 98 XMN2LLOL = XMN2LL 99 DO 800 I = 1,N 1 100 ICLSOL(I) = ICLUS(I) 1 101 800 CONTINUE 102 DO 850 IC = 1,K 1 103 XMEANOL(IC) = XMEAN(IC) 1 104 850 CONTINUE C 105 900 CONTINUE C 106 DO 1100 IC = 1,K C RECIPBET means 1/BETA. 1 107 RECIPBET(IC) = GAMMA(R(IC))*GAMMA(S(IC))/GAMMA(R(IC)+S(IC)) C 1 108 DO 1100 I = 1,N C C DSQ(IC) = ( XMEAN(IC) - X(I) )**2 C ZSQ = DSQ(IC)/VAR(IC) C IF ( ZSQ .LE. 174.673 ) GO TO 1000 C F(I,IC) = 0.0 C GO TO 1100 C1000 CONTINUE C C NOTE THAT in MIX1* the GAUSSIAN WAS USED HERE: C F(I,IC) =(EXP(-ZSQ/2.0))/(SQRT(2*PI*VAR(IC))) C 2 109 F(I,IC) = (X(I))**(R(IC)-1)*(1-X(I))**(S(IC)-1) C C CONST = RECIPBET(IC) 2 110 F(I,IC) = F(I,IC)*RECIPBET(IC) 2 111 1100 CONTINUE C 112 1200 CONTINUE C C COMPUTE POSTERIOR PROBABILITIES OF GROUP MEMBERSHIP: 113 DO 1300 I = 1,N 1 114 DENOM(I) = 0.0 1 115 DO 1300 IH=1,K 2 116 DENOM(I) = DENOM(I) + P(IH)*F(I,IH) 2 117 1300 CONTINUE 118 DO 1400 I = 1,N 1 119 DO 1400 IC=1,K C 2 120 IF ( DENOM(I) .EQ. 0.0 ) DENOM(I)=0.0001 C 2 122 PP(IC,I)= P(IC)*F(I,IC)/DENOM(I) 2 123 1400 CONTINUE 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 6 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C C COMPUTE NEW LABELS BY MAX POSTERIOR PROBABILITY: 124 DO 1600 I = 1,N 1 125 XMXPR(I) = PP(1,I) 1 126 ICLUS(I) = 1 1 127 DO 1600 IC = 2,K 2 128 IF ( PP(IC,I) .GT. XMXPR(I) ) GO TO 1500 2 129 GO TO 1600 2 130 1500 XMXPR(I) = PP(IC,I) 2 131 ICLUS(I) = IC 2 132 1600 CONTINUE C 133 IF (N .GE. 31) GO TO 1650 C C WRITE NEW LABELS: 134 WRITE (6,16000) 135 WRITE (6,18000) (IOTA(I), I=1,N) 136 WRITE (6,20000) (ICLUS(I), I=1,N) 137 1650 CONTINUE C C UPDATE CLUSTER PRIOR PROBABILITIES P(IC), MEANS XMEAN(IC) AND C VARIANCES VAR(IC): 138 WGSS = 0.0 139 DO 1800 IC = 1,K 1 140 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 141 SUM(IC) = 0.0 1 142 SS(IC) = 0.0 1 143 DO 1700 I = 1,N 2 144 XNC(IC) = XNC(IC) + PP(IC,I) 2 145 SUM(IC) = SUM(IC) + PP(IC,I)*X(I) 2 146 SS(IC) = SS(IC) + PP(IC,I)*X(I)*X(I) 2 147 1700 CONTINUE 1 148 IF ( XNC(IC) .EQ. 0.0 ) GO TO 2700 1 149 XMEAN(IC) = SUM(IC)/XNC(IC) 1 150 SSD(IC) = SS(IC) - SUM(IC)*SUM(IC)/XNC(IC) 1 151 VAR(IC) = SSD(IC)/(XNC(IC)) C IF ( VAR(IC) .LE. 0.0 ) VAR(IC) = 0.0001 1 152 SD(IC) = SQRT(VAR(IC)) 1 153 P(IC) = XNC(IC)/XN C P(IC) = estimate of the mixing prob of component IC. C Update the parms of the beta: 1 154 RPLUSS(IC) = XMEAN(IC)*(1-XMEAN(IC))/VAR(IC) - 1 1 155 R(IC) =XMEAN(IC)*RPLUSS(IC) 1 156 S(IC) = (1-XMEAN(IC) )*RPLUSS(IC) 1 157 RECIPBET(IC) = GAMMA(R(IC))*GAMMA(S(IC))/GAMMA(R(IC)+S(IC)) 1 158 WGSS = WGSS + SSD(IC) 1 159 1800 CONTINUE C C COUNT NUMBERS IN CLUSTERS: 160 DO 1900 IC = 1,K 1 161 NC(IC) = 0 1 162 1900 CONTINUE 163 DO 2000 I = 1,N 1 164 IGROUP = ICLUS(I) 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 7 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 1 165 NC(IGROUP) = NC(IGROUP) + 1 1 166 2000 CONTINUE C 167 VARHAT = WGSS/XN C 168 SUMLNF = 0.0 169 DO 2200 I=1,N 1 170 SUMPXF = 0.0 1 171 DO 2100 IC=1,K 2 172 SUMPXF = SUMPXF + P(IC)*F(I,IC) 2 173 2100 CONTINUE 1 174 SUMLNF = SUMLNF + DLOG(SUMPXF) 1 175 2200 CONTINUE C 176 XMN2LL = -2.0*SUMLNF C C 177 WGMS = WGSS/(N-K) 178 WRITE (6,40000) WGSS, XMN2LL, WGMS 179 STDERR = SQRT(WGMS) 180 WRITE (6,74000) STDERR C 181 KMINUS1 = K-1 182 DO 2300 IC=1,KMINUS1 1 183 ICPLUS1 = IC+1 C B(IC) IS approximate BOUNDARY BETWEEN classes c and c+1. 1 184 W1 = P(IC) 1 185 W2 = P(ICPLUS1) 1 186 B(IC) = P(IC)*XMEAN(IC) + P(ICPLUS1)*XMEAN(ICPLUS1) 1 187 B(IC) = B(IC)/( P(IC)+P(ICPLUS1) ) C BDYADJ = DLOG(P(ICPLUS1)) - DLOG(P(IC)) C BDYADJ = BDYADJ/(XMEAN(ICPLUS1) - XMEAN(IC)) C BDYADJ = VARHAT*BDYADJ C B(IC) = B(IC) - BDYADJ 1 188 2300 CONTINUE 189 WRITE (6,58000) ITER 190 WRITE (6,42000) (B(IC),IC=1,KMINUS1) 191 WRITE (6,38000) (XMEAN(IC),IC=1,K) C 192 WRITE (6,38100) ( P(IC),IC=1,K) C 193 WRITE (6,38300) (NC(IC),IC=1,K) C 194 IF (ITER .EQ. 1) GO TO 2500 195 DO 2400 I = 1,N C If any case changes cluster, another iteration will be performed. 1 196 IF (ICLUS(I) .EQ. ICLSOL(I)) GO TO 2400 1 197 GO TO 2500 1 198 2400 CONTINUE 199 WRITE (6,79000) C Even if no case changes cluster, another iteration will be C performed if the parameter estimates have changed much. 200 DO 2450 IC = 1,K 1 201 TEST=(XMEAN(IC)-XMEANOL(IC))/SQRT( XMEAN(IC)*(1-XMEAN(IC)) ) C TEST=(XMEAN(IC)-XMEANOL(IC))/XMEANOL(IC) C 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 8 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 1 202 TOLERANS = 0.05 C Another iteration will be performed iff. C (1) the max number of iterations has not been reached C and (2) the new means are not close to the old means. C or if (3) some case had changed clusters. C 1 203 IF (ABS(TEST) .LE. TOLERANS) GO TO 2450 1 204 GO TO 2500 1 205 2450 CONTINUE C If changes in all means are small, stop and write results. 206 GO TO 2900 207 2500 CONTINUE 208 ITER = ITER + 1 209 IF (ITER.GE. 5) GO TO 2600 210 GO TO 700 C 211 2600 WRITE (6,80000) 212 2700 CONTINUE 213 WRITE (6,88000) C C Next the results will be written: C 214 2900 CONTINUE C 215 WRITE (6,60000) ( P(IC), IC = 1,K ) 216 WRITE (6,62000) (NC(IC),IC=1,K) 217 WRITE (6,72000) (VAR(IC),IC=1,K) 218 WRITE (6,68000) (SD(IC), IC=1,K) 219 WRITE (6,68100) (R(IC), IC=1,K) 220 WRITE (6,68200) (S(IC), IC=1,K) C VARHAT IS MLE OF VARIANCE. 221 VARHAT = WGSS/N 222 WRITE (6,30000) VARHAT C C 223 IF (N .GE. 31) GO TO 3050 224 WRITE (6,78000) 225 DO 3000 I=1,N 1 226 WRITE (6,76000) X(I),ICLUS(I),(PP(IC,I),IC=1,K) 1 227 3000 CONTINUE C 228 3050 CONTINUE C C COMPUTE MODEL-SELECTION CRITERIA: C NO. PARAMETERS = K MEANS + K VARIANCES + (K-1) PROBS. 229 NOPARM = K + K + (K-1) C 230 WRITE (6,86000) NOPARM 231 AIC = XMN2LL + 2.0*NOPARM 232 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. 233 TEMP = 2*VARHAT**3 234 IF ( TEMP .LE. 0.0 ) TEMP = 0.0001 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 9 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 236 XKASH = SCH - DLOG(TEMP) C 237 WRITE (6,82000) AIC 238 WRITE (6,84000) SCH 239 WRITE (6,48000) XKASH C C C 240 STOP 241 10000 FORMAT(2X,I1) 242 12000 FORMAT(2X,I4) 243 14000 FORMAT(5X,F3.2,2X,F8.2,2X,F8.2) 244 16000 FORMAT(1X,'CLUSTERING') 245 18000 FORMAT(1X,'POINT: '/, (1X,40I3)) 246 20000 FORMAT(1X,'CLUSTER: '/, (1X,40I3)) 247 24000 FORMAT(1X,'DATA:'/) 248 26000 FORMAT(/1X,'INITIAL PRIOR PROBS, MEANS AND VARIANCES:') 249 28000 FORMAT(1X,I2,2X,3F11.2) 250 28100 FORMAT(3X,I2,2X,'R', F11.2, ' S', F11.2) 251 30000 FORMAT(1X, 'M.L. ESTIMATE OF COMMON VARIANCE = ',F14.5/) 252 32000 FORMAT(1X,'Program MIX1BETA CLUSPAC') 253 34000 FORMAT(1X,'COPYRIGHT (C) 1998 STANLEY LOUIS SCLOVE'/) 254 36000 FORMAT('1',1X,'K = ',I1,' CLUSTERS') 255 38000 FORMAT(1X,'MEANS: ',9F10.3) 256 38100 FORMAT(1X,'PROBS: ',9F10.3) C 257 38300 FORMAT(1X,'Ns: ',9F10.3) 258 40000 FORMAT(/1X,'WGSS = ',F14.4,' MINUS 2 LOG LIKELIHOOD = ', XF14.4, ' WGMS = ',F14.4/) 259 42000 FORMAT(1X,'BOUNDARIES:', 8X, 9F10.3) 260 44000 FORMAT(1X,'MINIMUM OF SAMPLE: ') 261 46000 FORMAT(1X,'MAXIMUM OF SAMPLE: ') 262 48000 FORMAT( 1X,'KASHYAP CRITERION = ', F14.4) 263 50000 FORMAT(/,1X,'MEAN = ', F14.4) 264 52000 FORMAT(1X, 'M.L. ESTIMATE OF VARIANCE = ',F14.5/) 265 54000 FORMAT(/1X,'SSDEVS = ',F14.4,' MINUS 2 LOG LIKELIHOOD = ', XF14.4, ' STDDEV = ',F14.4/) 266 56000 FORMAT('1','.................................................', X/,1X,'PROGRAM MIX1BETA CLUSPAC'/1X,'MIXTURE MODEL CLUSTERING'/ X,1X,'FOR UNIVARIATE DATA (DATA ON THE LINE)'/ X,1X,'WITH BETA DISTRIBUTIONS '/ X1X,'DEVELOPED AND PROGRAMMED BY DR. STANLEY L. SCLOVE' X/,1X,'VERSION 2.0 13-Dec-98 '/) 267 58000 FORMAT(1X,'ITERATION ', I2) 268 60000 FORMAT(1X,'PROBS:',5X,9F10.3) 269 62000 FORMAT(1X,'NUMBERS:',3X,9(I10,3X)) 270 64000 FORMAT(18A4) 271 66000 FORMAT(1X,'N = ',I3/) 272 68000 FORMAT(1X,'STD.DEVS.: ',9F10.3) 273 68100 FORMAT(1X,'Rs: ',9F10.3) 274 68200 FORMAT(1X,'Ss: ',9F10.3) 275 70000 FORMAT(1X,18A4/) 276 72000 FORMAT(1X,'VARIANCES: ',9F10.3) 277 74000 FORMAT(1X,'STD.ERROR=SQRT(WGMS) = ',F13.4/) 278 76000 FORMAT(1X,F8.2,2X,I2,3X,(8F6.3/)) 279 78000 FORMAT(1X/,1X,'DATA, LABELS AND PROBS OF CLUSTER MEMBERSHIP'/) 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 10 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 280 79000 FORMAT(/,1X,'NO CASE CHANGED CLUSTERS IN THIS ITERATION.') 281 80000 FORMAT(1X,'PROGRAM HAS NOT CONVERGED IN ', X' the prescribed number of ITERATIONS. STOP') 282 82000 FORMAT(1X,'AIC = ', F14.4 ) 283 84000 FORMAT(1X,'SCHWARZ CRITERION = ', F14.4 ) 284 86000 FORMAT(/1X,'NUMBER OF PARAMETERS = ',I4 ) 285 88000 FORMAT(1X,'NO OBSERVATIONS IN GROUP ',I3,'. STOP') 286 END 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 11 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 203 AIC R*4 71S 72F 231S 237F ALOG R*4 I 63 63 66 73 232 B R*4 A 4 186S 187F 187S 190F C R*4 A 2 UNREFERENCED DENOM R*8 AT 14 16 114S 116F 116S 120F 121S 122F DLOG R*8 I 77 174 236 DSQ R*4 A 2 UNREFERENCED F R*8 AT 10 15 109S 110F 110S 116F 122F 172F FMT R*4 A 5 26S 33 47F 50F 52F GAMMA GI 66 66 66 107 107 107 157 157 157 I I*4 32S 33 34F 35F 36F 37F 38F 39F 39F 40F 41F 42F 43B 47S 47S 62S 63F 63B 99S 100F 100B 108S 109F 109F 109F 110F 110B 113S 114F 116F 116F 116B 118S 120F 121F 122F 122F 122B 124S 125F 125F 126F 128F 128F 130F 130F 131B 135S 135S 136S 136S 143S 144F 145F 145F 146F 146F 146B 163S 164B 169S 172B 195S 196F 196B 225S 226F 226F 226B IC I*4 81S 82 82 82 83F 83F 83F 83F 84F 84F 84F 85F 85F 85B 88S 89F 89F 89F 89F 90F 90F 90B 102S 103F 103B 106S 107F 107F 107F 107F 107F 109F 109F 109F 110F 110F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 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 110B 119S 122F 122F 122B 127S 128F 130F 131B 139S 140F 141F 142F 144F 144F 144F 145F 145F 145F 146F 146F 146F 148F 149F 149F 149F 150F 150F 150F 150F 150F 151F 151F 151F 152F 152F 153F 153F 154F 154F 154F 154F 155F 155F 155F 156F 156F 156F 157F 157F 157F 157F 157F 158B 160S 161B 171S 172F 172B 182S 183F 184F 186F 186F 186F 187F 187F 187B 190S 190S 191S 191S 192S 192S 193S 193S 200S 201F 201F 201F 201B 215S 215S 216S 216S 217S 217S 218S 218S 219S 219S 220S 220S 226S 226S ICLSOL I*4 A 9 100S 196F ICLUS I*4 A 1 100F 126S 131S 136F 164F 196F 226F ICPLUS1 I*4 183S 185F 186F 186F 187F IGROUP I*4 164S 165F 165F IH I*4 115S 116F 116B INTGER I*4 92S 93F 93B IOTA I*4 A 1 93S 135F ITER I*4 95S 97F 189F 194F 208F 208S 209F K I*4 79S 80F 81F 88F 102F 106F 115F 119F 127F 139F 160F 171F 177F 181F 191F 192F 193F 200F 215F 216F 217F 218F 219F 220F 226F 229F 229F 229F KMINUS1 I*4 181S 182F 190F N I*4 23S 24F 25F 32F 46F 47F 62F 92F 99F 108F 113F 118F 124F 133F 135F 136F 143F 163F 169F 177F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 13 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 195F 221F 223F 225F NC I*4 A 4 161S 165F 165S 193F 216F NOPARM I*4 70S 71F 73F 229S 230F 231F 232F P R*8 AT 11 15 82S 89F 116F 122F 153S 172F 184F 185F 186F 186F 187F 187F 192F 215F PP R*8 AT 12 15 122S 125F 128F 130F 144F 145F 146F 226F R R*4 A 2 84S 90F 107F 107F 109F 155S 157F 157F 219F RECIPBET R*4 A 4 107S 110F 157S RPLUSS R*4 A 2 83S 84F 85F 154S 155F 156F S R*4 A 2 85S 90F 107F 107F 109F 156S 157F 157F 220F SCH R*4 73S 77F 232S 236F 238F SD R*4 A 7 152S 218F SQRT GI 68 152 179 201 SS R*8 AT 6 15 142S 146F 146S 150F SSD R*8 AT 6 15 150S 151F 158F SSDEVS R*8 T 17 29S 54S 55F 69F STDDEV R*4 68S 69F STDERR R*4 179S 180F SUM R*8 AT 2 15 141S 145F 145S 149F 150F 150F SUMLNF R*8 T 16 168S 174F 174S 176F SUMPXF R*8 T 16 170S 172F 172S 174F SUMSQS R*8 T 15 28S 39F 39S 54F TEMP R*8 T 16 63S 64F 74S 75F 76S 77F 233S 234F 235S 236F TEST R*4 201S 203F TITLE R*4 A 3 18S 22F TOLERANS R*4 202S 203F TOTAL R*8 T 15 27S 38F 38S 53F 54F 54F VAR R*8 AT 8 15 82S 83F 89F 151S 152F 154F 217F VARHAT R*4 55S 57F 58F 68F 74F 167S 221S 222F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 14 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 233F WGMS R*4 177S 178F 179F WGSS R*8 T 17 138S 158F 158S 167F 177F 178F 221F W1 R*4 184S W2 R*4 185S X R*4 A 1 33S 34F 35S 36F 37S 38F 39F 39F 40F 41F 42F 43F 47F 63F 63F 109F 109F 145F 146F 146F 226F XBAR R*4 53S 56F 58F 58F 59F 60F XC R*4 58S 59F 60F XKASH R*4 77S 78F 236S 239F XLL R*8 T 16 61S 64F 64S 66F 66S 67F XMAX R*4 30S 42F 43S 52F XMEAN R*8 AT 4 17 82S 83F 83F 84F 85F 89F 103F 149S 154F 154F 155F 156F 186F 186F 191F 201F 201F 201F XMEANOL R*8 AT 4 17 103S 201F XMIN R*4 31S 40F 41S 50F XMNDSQ R*4 A 1 UNREFERENCED XMN2LL R*8 T 17 67S 69F 71F 73F 98F 176S 178F 231F 232F XMN2LLOL R*8 T 17 98S XMXPR R*4 A 13 125S 128F 130S XN R*4 24S 53F 54F 55F 73F 153F 167F 232F XNC R*8 AT 11 15 140S 144F 144S 148F 149F 150F 151F 153F XR R*4 59S 63F 66F 66F XS R*4 60S 63F 66F 66F 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 15 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 105 65 62 300 44 32 350 B 48 46 400 86 81 500 91 88 600 94 92 700 B 96 210 800 101 99 850 104 102 900 B 105 97 1100 111 106 108 1200 112 UNREFERENCED 1300 117 113 115 1400 123 118 119 1500 B 130 128 1600 B 132 124 127 129 1650 B 137 133 1700 147 143 1800 159 139 1900 162 160 2000 166 163 2100 173 171 2200 175 169 2300 188 182 2400 B 198 195 196 2450 B 205 200 203 2500 B 207 194 197 204 2600 B 211 209 2700 B 212 148 2900 B 214 206 3000 227 225 3050 B 228 223 10000 NF 241 79 12000 NF 242 23 14000 NF 243 82 16000 NF 244 134 18000 NF 245 135 20000 NF 246 136 24000 NF 247 45 26000 NF 248 87 28000 NF 249 89 28100 NF 250 90 30000 NF 251 222 32000 NF 252 20 34000 NF 253 21 36000 NF 254 80 38000 NF 255 191 38100 NF 256 192 38300 NF 257 193 40000 NF 258 178 1LEVEL 2.5.0 (JUNE 1991) VS FORTRAN DEC 16, 1998 07:53:08 NAME:MAIN# PAGE: 16 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 42000 NF 259 190 44000 NF 260 49 46000 NF 261 51 48000 NF 262 78 239 50000 NF 263 56 52000 NF 264 57 54000 NF 265 69 56000 NF 266 19 58000 NF 267 189 60000 NF 268 215 62000 NF 269 216 64000 NF 270 18 26 66000 NF 271 25 68000 NF 272 218 68100 NF 273 219 68200 NF 274 220 70000 NF 275 22 72000 NF 276 217 74000 NF 277 180 76000 NF 278 226 78000 NF 279 224 79000 NF 280 199 80000 NF 281 211 82000 NF 282 72 237 84000 NF 283 238 86000 NF 284 230 88000 NF 285 213 0*STATISTICS* SOURCE STATEMENTS: 279, PROGRAM SIZE: 215040 BYTES, PROGRAM NAME: MAIN#, PAGE: 1 *STATISTICS* NO DIAGNOSTICS GENERATED. **MAIN#** END OF COMPILATION 1 ****** TIME STAMP: 98.35007.53.08