C 10000 FORMAT('1','**********************************************'/// X' MIXPDTA CLUSPAC - MIXTURE-MODEL CLUSTERING OF CASES'// X' VARYING COVARIANCE MATRICES '/ X' AUTOMATIC SETTING OF INITIAL PARAMETER ESTIMATES '// X' Copyright (C) 1991 Stanley Louis Sclove '// X' Developed and programmed by: '// X19X,'Prof. Stanley L. Sclove, Ph.D. Phone 312-996-2681'/ X19X,'Information & Decision Sciences Dept. M/C 294 '/ X19X,'University of Illinois at Chicago '/ X19X,'601 S. Morgan St. '/ X19X,'Chicago IL 60607-7124 '// X' Program Version: 2.8 01-Jan-90 '/ X' (VM/CMS) ') C 11097 FORMAT(/' MIXPCMA CLUSPAC - MIXTURE-MODEL CLUSTERING OF CASES'/) 10999 FORMAT( ' PROBLEM TITLE IS' ) 11096 FORMAT(/' NUMBER OF CLUSTERS TO REPORT. . . . . . . . 1 TO', XI2/) 21000 FORMAT(/1X,'MAXIMUM FOR EACH VARIABLE: ',/) 10050 FORMAT(/' ................................................', X'................................'//) 11111 FORMAT(/' FIRST FOUR DATA POINTS') 11112 FORMAT(' ----- ---- ---- ------'/) 19000 FORMAT( 1X,'MINIMUM FOR EACH VARIABLE: ',/) 10100 FORMAT(' MEAN VECTOR: ',(9F8.2/)/) 10200 FORMAT(/' COVARIANCE MATRIX: '/) 10500 FORMAT(//' STATISTICS FOR WHOLE (UNCLUSTERED) SAMPLE: '/) 11001 FORMAT(1X,18A4/) 11002 FORMAT(/1X,18A4/) 11100 FORMAT(1X,I2,'/',I2,'/',I4,3X,'AT ',I2,':',I2,':',I2/) 44444 FORMAT(' CLUSTER ',I2) 44445 FORMAT(1X,A6) 10900 FORMAT(/' MINUS 2 LOG LIKELIHOOD = ',F13.1) 11000 FORMAT(//1X,18A4 ) 11500 FORMAT(' DET = ',E15.7, ' IDET = ',I2/) 12000 FORMAT(2X,I4) 13000 FORMAT( ' NUMBER OF OBSERVATIONS (SAMPLE SIZE), N . . ',I3/) 14000 FORMAT(//1X,'CLUSTERING:') 15000 FORMAT(18A4) 16000 FORMAT(/,1X,'CASES AND LABELS:--'/) 17000 FORMAT(' PREC.MX. OF CLUSTER ',I2,' HAS BEEN SET TO POOLED ', X'PREC.MX.') 18000 FORMAT(1X, (7(1X,A6,': ',I2)/) ) 21999 FORMAT(' '//) 22000 FORMAT(1X,'MEAN VECTOR FOR CLUSTER ',I2,': ') 26000 FORMAT(/, 1X,'K = ',I2,' CLUSTERS'/) 30000 FORMAT(1X,'MEAN VECTOR FOR CLUSTER ',I2,': '/,(9F13.3/)) 32000 FORMAT(/,1X,'ITERATION ', I2) C WRITE(6,46000) K,XM2LL(K),IPARM(K),AICVEC(K),AICPPH(K), C XSCHVEC(K),SCHPPH(K) 14560 FORMAT(1X,'THE DENOMINATOR OF THE POSTERIOR PROBS FOR ', X1X,A6, ' IS 0, SO THE POSTERIOR PROBS PP(IC,I), IC=1,...,K,'/ X' HAVE BEEN SET EQUAL TO THOSE FROM THE ', X' PRECEDING ITERATION SO THAT THE I-TH CASE WILL STILL '/ X' HAVE WEIGHT IN SUCCEEDING COMPUTATIONS.') 36000 FORMAT(/,1X,'NUMBERS IN CLUSTERS:',3X,9(I5,3X)//) 40000 FORMAT(/,1X,'MIXING PROBABILITIES:',3X,9(F4.2,3X)/) 42000 FORMAT(//,1X,'COMMON COVARIANCE MATRIX (MLE):',/ ) 34000 FORMAT(1X,'MODEL SELECTION CRITERIA'// X' MODEL-SELECTION CRITERIA '/ X' NUMBER MINUS 2 NUMBER (PPH=POST.PROB.OF MODEL):'/ X' OF TIMES OF AKAIKE SCHWARZ '/ X' CLUSTERS, LOG OF PARAM- ------------ ------------ '/ X' K LIKELIHOOD ETERS VALUE PPH VALUE PPH '/ X' --------- ---------- ----- ------- --- ------- ---'//) C 9 XXXX.X 143 XXXX.X .XX XXXX.X .XX 46000 FORMAT(5X,I1,6X,F8.1, 4X,I4,3X,F7.1,1X,F4.2,3X,F7.1,1X,F4.2) 48000 FORMAT(1X,5F13.3) 50000 FORMAT(1X,A6,': ',I2,(8F13.3/)) 52000 FORMAT(//,1X,'CLUSTER MEMBERSHIP OF EACH CASE:'/ X1X,'CASE LABEL DATA'/ ) 54000 FORMAT(3X,I2) 56000 FORMAT(/' NUMBER OF VARIABLES . . . . . . . . . . . . ',I3 ) 58000 FORMAT(/1X,'CONVERGENCE: NO CASE CHANGED CLUSTERS AFTER ', X'ITERATION ',I2,'.'//) 60000 FORMAT(/,1X,'JFLG = ',I2,'. IF JFLG=0, COMPUTATION OF DET', X' WENT WELL; OTHERWISE, THERE WAS TROUBLE OR MATRIX WAS ', X'ILL-CONDITIONED.'//) 62000 FORMAT(/1X,'DET = ',E13.5,' IDET = ',I3,5X, X'ACTUAL DET. = DET*10**IDET',//) 70000 FORMAT(1X,' AIC = ', F13.1) 78000 FORMAT(1X,' SCHWARZ CRITERION = ', F13.1/) 72000 FORMAT( ' NUMBER OF PARAMETERS = ',I4) 74000 FORMAT(/1X,'EXPECTED NO. OF OBSERVATIONS IN GROUP ',I3,'IS ', X' LESS THAN 0.5. ', X' END COMPUTATION FOR THIS NUMBER OF CLUSTERS. '/) 76000 FORMAT(/' POSTERIOR PROB OF GROUP MEMBERSHIP FOR 1ST 4 CASES:') 82000 FORMAT(1X,29F5.2) 88000 FORMAT(' ALGORITHM HAS NOT CONVERGED IN 20 ITERATIONS. STOP') 96000 FORMAT(/' NGFLAG=1: FOR SOME K, AT LEAST ONE CLUSTER HAS '/ X' AN EXP.NO. OF OBSERVATIONS LESS TH AN 0.5 .') 98000 FORMAT(//1X,'PROGRAM ENDED NORMALLY.') END C