1 J E S 2 J O B L O G -- S Y S T E M U I C C -- N O D E U I C M V S - 18.58.20 JOB 369 IEF677I WARNING MESSAGE(S) FOR JOB CMS37331 ISSUED 18.58.20 JOB 369 $HASP373 CMS37331 STARTED - INIT 1 - CLASS A - SYS UICC 18.58.42 JOB 369 +AFB240I VABEX : ABEND CODE IS: SYSTEM 0C4, USER 0. 18.58.42 JOB 369 +AFB240I VABEX : PSW=078D000400150032 ENTRY POINT=117010. 18.58.42 JOB 369 +AFB240I VABEX : REGS 0 - 3 00000032 00117B0C 001179B2 C1497905 18.58.42 JOB 369 +AFB240I VABEX : REGS 4 - 7 587A0E55 0013BD02 00000008 00118D08 18.58.42 JOB 369 +AFB240I VABEX : REGS 8 -11 00000000 8014FC76 00114FE0 001442B0 18.58.42 JOB 369 +AFB240I VABEX : REGS 12-15 8014FBB0 00144F40 6013BD5C 8014FC5E 18.58.42 JOB 369 +AFB240I VABEX : FRGS 0 & 2 42320000 00000000 3E96226B 15979FD5 18.58.42 JOB 369 +AFB240I VABEX : FRGS 4 & 6 B7104D7D 10C00000 3E6F80D6 950CEC80 18.58.42 JOB 369 +IEW1991 ERROR - USER PROGRAM HAS ABNORMALLY TERMINATED 18.58.42 JOB 369 IEF450I CMS37331 GO - ABEND S000 U0240 18.58.43 JOB 369 $HASP395 CMS37331 ENDED 0------ JES2 JOB STATISTICS ------ - 08 FEB 90 JOB EXECUTION DATE - 911 CARDS READ - 1,464 SYSOUT PRINT RECORDS - 0 SYSOUT PUNCH RECORDS - 0.38 MINUTES EXECUTION TIME 0 1 //CMS37331 JOB JOB 369 ***ACCOUNT U37331 *CMSWYL ***KEY CARD IGNORED - DUPLICATE ***JOBPARM Q=F *CMSWYL ***ROUTE PRINT BSB *CMSWYL ***JOBPARM T=(0,40),L=3,R=2048 1. ***ROUTE PRINT UICVM.U37331 2. 2 // EXEC FORTVCLG 3. *********************************************************************// 0000000 *** *// 0000000 *** VS FORTRAN 2.3 (LANGLVL(77)) *// 0000000 *** *// 0000000 *** USES IMSL LIBRARY EDITION 10 *// 0000000 *** *// 0000000 *** NEW DEFAULTS: XREF CROSS REFERENCE LISTING *// 0000000 *** SXM XREF FORMATTED FOR TERMINAL DISPLAY *// 0000000 *** USE // EXEC FORTVCXX,XREF=NOXREF *// 0000000 *** OR // EXEC FORTVCXX,SXM=NOSXM TO CANCEL DEFAULTS *// 0000000 *** *// 0000000 *** NOTE: RELEASE 2 LOADS SUBROUTINES AT RUN TIME. *// 0000000 *** FOR LINKING SUBROUTINES SEE INFORM DOCUMENT 904 AND USE: *// 0000000 *** // EXEC FORTVXX,LIB='SYS1.VSFORT.VLNKMLIB.VER23' *// 0000000 *** *// 0000000 *********************************************************************// 0000000 3 XXFORTVCLG PROC FVPGM=FORTVS,FVREGN=800K,FVPDECK=NODECK, +0000000 XX FVPOLST=NOLIST,FVPOPT=0,FVTERM='DUMMY', +0000000 XX FVLNSPC='3200,(25,6)', +0000000 XX GOF5DD='DDNAME=SYSIN', +0000000 XX GOF6DD='SYSOUT=A', +0000000 XX GOF7DD='SYSOUT=B',GOREGN=100K, +0000000 XX LIB='UICCAPP.DUMMY',VER=23,LANGLVL='LANGLVL(77)', +0000000 XX IMSL=ED10,XREF=XREF,SXM=SXM,INTERF='UICCAPP.DUMMY' 0000000 *** 0000000 *** PARAMETER DEFAULT-VALUE USAGE 0000000 *** 0000000 *** FVPGM FORTVS COMPILER NAME 0000000 *** FVPDECK NODECK COMPILER DECK OPTION 0000000 *** FVPOLST NOLIST COMPILER LIST OPTION 0000000 *** FVPOPT 0 COMPILER OPTIMIZATION 0000000 *** FVTERM DUMMY FORT.SYSTERM OPERAND 0000000 *** FVLNSPC 3200,(25,6) FORT.SYSLIN SPACE 0000000 *** GOF5DD DDNAME=SYSIN GO.FT05F001 OPERAND 0000000 *** GOF6DD SYSOUT=A GO.FT06F001 OPERAND 0000000 *** GOF7DD SYSOUT=B GO.FT07F001 OPERAND 0000000 *** 0000000 4 XXFORT EXEC PGM=&FVPGM,REGION=&FVREGN, 0000000 XX PARM='&XREF,&SXM,&FVPDECK,&FVPOLST,OPT(&FVPOPT),&LANGLVL' 0000000 5 XXSTEPLIB DD DSN=SYS1.VSFORT.FORTVS.VER&VER,DISP=SHR 0000000 6 XXSYSPRINT DD SYSOUT=A,DCB=BLKSIZE=3429 0000000 7 XXSYSTERM DD &FVTERM 0000000 8 XXSYSPUNCH DD SYSOUT=B,DCB=BLKSIZE=3440 0000000 9 XXSYSLIN DD DSN=&&UICOBJ,DISP=(MOD,PASS),UNIT=SCRATCH, + 0000000 XX SPACE=(&FVLNSPC),DCB=BLKSIZE=3200 0000000 10 //FORT.SYSIN DD * 4. 11 XXGO EXEC PGM=LOADER,REGION=&GOREGN,COND=(7,LT,FORT), +0000000 XX PARM='LET,NORES,EP=MAIN' 0000000 12 XXSYSLIN DD DSN=&&UICOBJ,DISP=(OLD,DELETE) 0000000 13 XX DD DDNAME=OBJMODS 0000000 14 XXSYSLOUT DD SYSOUT=A 0000000 15 XXSYSLIB DD DSN=&LIB,DISP=SHR 0000000 16 XX DD DSN=SYS1.VSFORT.VFORTLIB.VER&VER,DISP=SHR 0000000 17 XX DD DSN=&INTERF,DISP=SHR 0000000 18 XX DD DSN=UICCAPP.IMSLIB.&IMSL,DISP=SHR 0000000 19 XX DD DSN=UICCAPP.EISPACK,DISP=SHR 0000000 20 XX DD DSN=UICCAPP.MINPACK,DISP=SHR 0000000 21 XX DD DSN=UICCAPP.LINPACK,DISP=SHR 0000000 22 XX DD DSN=UICCAPP.FORTLIB2,DISP=SHR 0000000 23 XX DD DSN=UICCAPP.HARLIB,DISP=SHR 0000000 24 XXSTEPLIB DD DSN=SYS1.VSFORT.VFLODLIB.VER&VER,DISP=SHR 0000000 25 XXFT05F001 DD &GOF5DD 0000000 26 XXFT06F001 DD &GOF6DD 0000000 27 XXFT07F001 DD &GOF7DD 0000000 28 //GO.SYSIN DD * 867. 29 //SYSIN DD * GENERATED STATEMENT 30 //SYSIN DD * GENERATED STATEMENT 0 STMT NO. MESSAGE - 4 IEF653I SUBSTITUTION JCL - PGM=FORTVS,REGION=800K, 4 IEF653I SUBSTITUTION JCL - PARM='XREF,SXM,NODECK,NOLIST,OPT(0),LANGLVL(77)' 5 IEF653I SUBSTITUTION JCL - DSN=SYS1.VSFORT.FORTVS.VER23,DISP=SHR 7 IEF653I SUBSTITUTION JCL - DUMMY 9 IEF653I SUBSTITUTION JCL - SPACE=(3200,(25,6)),DCB=BLKSIZE=3200 11 IEF653I SUBSTITUTION JCL - PGM=LOADER,REGION=100K,COND=(7,LT,FORT), 15 IEF653I SUBSTITUTION JCL - DSN=UICCAPP.DUMMY,DISP=SHR 16 IEF653I SUBSTITUTION JCL - DSN=SYS1.VSFORT.VFORTLIB.VER23,DISP=SHR 17 IEF653I SUBSTITUTION JCL - DSN=UICCAPP.DUMMY,DISP=SHR 18 IEF653I SUBSTITUTION JCL - DSN=UICCAPP.IMSLIB.ED10,DISP=SHR 24 IEF653I SUBSTITUTION JCL - DSN=SYS1.VSFORT.VFLODLIB.VER23,DISP=SHR 25 IEF653I SUBSTITUTION JCL - DDNAME=SYSIN 26 IEF653I SUBSTITUTION JCL - SYSOUT=A 27 IEF653I SUBSTITUTION JCL - SYSOUT=B 30 IEF686I DDNAME REFERRED TO ON DDNAME KEYWORD IN PRIOR STEP WAS NOT RESOLVED IEF142I CMS37331 FORT - STEP WAS EXECUTED - COND CODE 0000 IEF373I STEP /FORT / START 90039.1858 IEF374I STEP /FORT / STOP 90039.1858 CPU 0MIN 00.70SEC SRB 0MIN 00.12SEC VIRT 1816K SYS 308K CPU TIME(SEC)=.70; REGION REQUESTED(K)=2,048; SRB TIME(SEC)=.12; REGION USED(K)=1,816; DISK EXCPS=297; REAL TIME(SEC)=13.81; PAGES-IN= 00; SWAPS= 00; PAGES-OUT= 00; AFB240I VABEX : ABEND CODE IS: SYSTEM 0C4, USER 0. AFB240I VABEX : PSW=078D000400150032 ENTRY POINT=117010. AFB240I VABEX : REGS 0 - 3 00000032 00117B0C 001179B2 C1497905 AFB240I VABEX : REGS 4 - 7 587A0E55 0013BD02 00000008 00118D08 AFB240I VABEX : REGS 8 -11 00000000 8014FC76 00114FE0 001442B0 AFB240I VABEX : REGS 12-15 8014FBB0 00144F40 6013BD5C 8014FC5E AFB240I VABEX : FRGS 0 & 2 42320000 00000000 3E96226B 15979FD5 AFB240I VABEX : FRGS 4 & 6 B7104D7D 10C00000 3E6F80D6 950CEC80 IEW1991 ERROR - USER PROGRAM HAS ABNORMALLY TERMINATED IEF472I CMS37331 GO - COMPLETION CODE - SYSTEM=000 USER=0240 IEF242I ALLOC. FOR CMS37331 GO AT ABEND IEF237I 549 ALLOCATED TO SYSLIN IEF237I DMY ALLOCATED TO IEF237I JES2 ALLOCATED TO SYSLOUT IEF237I 357 ALLOCATED TO SYSLIB IEF237I 352 ALLOCATED TO IEF237I 357 ALLOCATED TO IEF237I 352 ALLOCATED TO IEF237I C43 ALLOCATED TO IEF237I 352 ALLOCATED TO IEF237I 352 ALLOCATED TO IEF237I 352 ALLOCATED TO IEF237I 352 ALLOCATED TO IEF237I 352 ALLOCATED TO STEPLIB IEF237I JES2 ALLOCATED TO FT05F001 IEF237I JES2 ALLOCATED TO FT06F001 IEF237I JES2 ALLOCATED TO FT07F001 IEF237I JES2 ALLOCATED TO SYSIN IEF237I JES2 ALLOCATED TO SYSIN IEF285I SYS90039.T185820.RA000.CMS37331.UICOBJ DELETED IEF285I VOL SER NOS= SCR05 . IEF285I JES2.JOB00369.SO0107 SYSOUT IEF285I UICCAPP.DUMMY KEPT IEF285I VOL SER NOS= ARCH32. IEF285I SYS1.VSFORT.VFORTLIB.VER23 KEPT IEF285I VOL SER NOS= MVS004. IEF285I UICCAPP.DUMMY KEPT IEF285I VOL SER NOS= ARCH32. IEF285I UICCAPP.IMSLIB.ED10 KEPT IEF285I VOL SER NOS= MVS004. IEF285I UICCAPP.EISPACK KEPT IEF285I VOL SER NOS= MVS005. IEF285I UICCAPP.MINPACK KEPT IEF285I VOL SER NOS= MVS004. IEF285I UICCAPP.LINPACK KEPT IEF285I VOL SER NOS= MVS004. IEF285I UICCAPP.FORTLIB2 KEPT IEF285I VOL SER NOS= MVS004. IEF285I UICCAPP.HARLIB KEPT IEF285I VOL SER NOS= MVS004. IEF285I SYS1.VSFORT.VFLODLIB.VER23 KEPT IEF285I VOL SER NOS= MVS004. IEF285I JES2.JOB00369.SI0102 SYSIN IEF285I JES2.JOB00369.SO0108 SYSOUT IEF285I JES2.JOB00369.SO0109 SYSOUT IEF285I JES2.JOB00369.SI0103 SYSIN IEF285I JES2.JOB00369.SI0104 SYSIN IEF373I STEP /GO / START 90039.1858 IEF374I STEP /GO / STOP 90039.1858 CPU 0MIN 00.31SEC SRB 0MIN 00.08SEC VIRT 2048K SYS 312K CPU TIME(SEC)=.31; REGION REQUESTED(K)=2,048; SRB TIME(SEC)=.08; REGION USED(K)=2,048; DISK EXCPS=202; REAL TIME(SEC)=8.98; PAGES-IN= 00; SWAPS= 00; PAGES-OUT= 00; IEF375I JOB /CMS37331/ START 90039.1858 IEF376I JOB /CMS37331/ STOP 90039.1858 CPU 0MIN 01.01SEC SRB 0MIN 00.20SEC 0 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 PAGE: 1 0REQUESTED OPTIONS (EXECUTE): XREF,SXM,NODECK,NOLIST,OPT(0),LANGLVL(77) 0OPTIONS IN EFFECT: NOLIST NOMAP XREF NOGOSTMT NODECK SOURCE TERM OBJECT FIXED NOTRMFLG SRCFLG NOSYM NORENT SDUMP(ISN) SXM NOVECTOR IL(DIM) NOTEST NODC NOICA NODIRECTIVE NODBCS NOSAA OPT(0) LANGLVL(77) NOFIPS FLAG(I) AUTODBL(NONE) NAME(MAIN) LINECOUNT(60) CHARLEN(500) 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C 5. C 6. C ..............................................................C 7.1 C C 8. C C 9. C CMS FILE: EVALUATE CLUSPAC C 10. C VERSION 1.4 30-JAN-90 C 11. C C 12. C PROGRAMMED BY: C 13. C C 14. C DR. STANLEY L. SCLOVE 312/996-2681 C 15. C INFORMATION AND DECISION SCIENCE DEPT. M/C 294 C 17. C COLLEGE OF BUSINESS ADMINISTRATION C 18. C UNIVERSITY OF ILLINOIS AT CHICAGO C 19. C BOX 4348, CHICAGO, IL 60680 C 20. C C 21. C C 22. C COPYRIGHT (C) 1990 STANLEY L. SCLOVE. ALL RIGHTS RESERVED. C 23. C C 24. C ..............................................................C 24.01 C C 24.1 C C 24.2 C Program EVALUATE CLUSPAC evaluates the likelihood C 25. C of classified (labeled) data under the assumption of C 25.1 C Gaussian class-conditional distributions. The evaluation C 26. C is made both for a common covariance matrix and varying C 27. C covariance matrices. Values of model-selection criteria C 28. C are computed. These may be compared with those obtained C 29. C by clustering. C 30. C C 31. C ..............................................................C 32. C C 40. C Input: C 41. C ----- C 42. C Number of groups (K), data and labels. C 43. C C 49. C Program restrictions (can be modified): C 50. C -------------------------------------- C 51. C N, sample sizes, at most 1000; C 52. C IP, number of variables, at most 20; C 53. C K, number of groups, at most 29; C 54. C C 56. C C 57. C Subroutines called: C 58. C MATEQ, which calls MATDT C 59. C C 60. C IV is a work array for subroutine MATEQ. C 61. C C 62. C C 63. C ..............................................................C 65. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 2 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C C 66. C C 67. C C 68. C CONTROL CARDS: C 69. C C 70. C //GO.SYSIN DD * 70.1 C FISHER IRIS DATA 70.2 C K=03 70.21 C IP=04 70.3 C (4(1X,F3.1),1X,I1) 70.4 C N=0050 70.6 C C 70.7 C THESE ARE: 70.8 C JCL "GO"-STEP CARD 70.9 C DATASET TITLE C 71. C K, NUMBER OF GROUPS, IN FORMAT (2X,I2) C 71.1 C NO. OF VARIABLES, IP, IN FORMAT (3X,I2) 73. C DATFMT, IN FORMAT (18A4), E.G., (4F4.1) C 73.1 C "DATFMT" WILL ALSO BE USED FOR OUTPUT, SO ALLOW AT LEAST C 73.2 C ONE BLANK AT THE BEGINNING FOR CARRIAGE CONTROL. C 73.3 C N FOR FIRST SAMPLE, IN FORMAT (2X,I4) 73.4 C DATA FOR FIRST SAMPLE 73.5 C 73.6 C (DATA, ONE CASE AT A TIME, IN FORMAT SPECIFIED BY DATFMT) 73.7 C 73.8 C N FOR SECOND SAMPLE, IN FORMAT (2X,I4) 73.9 C DATA FOR SECOND SAMPLE 74. C . 74.1 C . 74.2 C . 74.3 C 74.4 C N FOR K-TH SAMPLE, IN FORMAT (2X,I4) 74.5 C DATA FOR K-TH SAMPLE 74.6 C 74.7 C C 74.9 C C 75. C C 75.1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 75.2 C 75.4 C 75.5 C 75.6 1 DIMENSION X(1000,20),SUM(29,20),XMEAN(20) 93. C ELEMENTS OF X ARE X(I,JV), I=1,2,...,N, JV=1,2,...,IP. 93.1 2 DIMENSION TITLE(18) 98. 3 DIMENSION DATFMT(18) 98.1 4 DIMENSION NG(29),TOTAL(29),SUMSQS(20,20) 101. 5 DIMENSION SSD(29,20,20),SIGMA(29,20,20),SSDEVS(20,20) 103. 6 DIMENSION WGSS(20,20) 104. 7 DIMENSION VARHAT(20,20),COMSIG(20,20),COVMX(20,20) 105. 8 DIMENSION XMIN(20),XMAX(20) 106. 9 DIMENSION IV(20,20) 107. 10 DIMENSION XLLIKE(29) 108.2 C 110. 11 DOUBLE PRECISION SUM 111. 12 DOUBLE PRECISION WGSS,SSD 112. 13 DOUBLE PRECISION VARHAT,COVMX,TOTAL,SUMSQS,SSDEVS 113. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 3 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 14 DOUBLE PRECISION P 114. 15 DOUBLE PRECISION DET,TRUDET 115. 16 DOUBLE PRECISION DSQ 116. 17 DOUBLE PRECISION DEVV,DEVW 118. 18 DOUBLE PRECISION F 119. 19 DOUBLE PRECISION XLLIKE 119.1 C 120. C 121. C 122. C FLOW OF PROGRAM: 123. C 124. C 125. C FOR EACH SAMPLE:-- 125.1 C READ DATA. 126. C COMPUTE MEAN VECTOR AND COVARIANCE MATRIX. 127. C EVALUATE LIKELIHOOD. 127.1 C L = EXP(-NP/2)/( (2*PI)**NP/2*DET(COV MX)**N/2 ) 127.2 C EVALUATE LIKELIHOOD WITH VARYING COVARIANCE MATRICES BY 128. C MULTIPLYING LIKELIHOODS FOR THE SEPARATE SAMPLES. 129. C COMPUTE MODEL-SELECTION CRITERIA. 130. C 130.1 C POOL TO OBTAIN COMMON COVARIANCE MATRIX. 131. C EVALUATE LIKELIHOOD WITH COMMON COVARIANCE MATRIX. 132. C COMPUTE MODEL-SELECTION CRITERIA. 140. C PRINT RESULTS. 141. C 142. C 142.1 C 142.2 20 PI = .314159265358979324D+01 143.1 C 144. 21 READ (5,36000) TITLE 145. C 146. C WRITE PROGRAM INFORMATION. 147. 22 WRITE (6,10000) 148. 23 WRITE (6,10050) 148.1 24 WRITE (6,40000) TITLE 149. C 150. C READ K, NUMBER OF GROUPS. 151. 25 READ (5,11000) K 152. 26 WRITE (6,26000) K 153. C READ NUMBER OF VARIABLES, IP. 155. 27 READ (5,54000) IP 156. 28 WRITE (6,13000) IP 157. 29 WRITE (6,10050) 157.1 C 158. C READ DATA FORMAT. 159. 30 READ (5,36000) DATFMT 160. C 160.1 C 160.2 31 CALL XUFLOW(0) 160.3 C 161. C DO FOR EACH SAMPLE: 161.01 C 161.02 C 161.021 32 DO 4100 ISAMPL = 1,K 161.03 1 33 WRITE (6,70000) ISAMPL 161.04 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 4 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C READ SAMPLE SIZE, N. 161.1 1 34 READ (5,12000) N 161.2 1 35 XN = N 161.3 1 36 NG(ISAMPL) = N 161.31 1 37 WRITE (6,15000) N 161.4 C READ DATA. 162. 1 38 DO 500 I = 1,N 163. 2 39 READ (5,DATFMT) (X(I,JV), JV = 1,IP) 164. 2 40 IF (I .EQ. 1) GO TO 100 165. 2 41 GO TO 300 166. 2 42 100 CONTINUE 167. C COMPUTE MINIMA AND MAXIMA: 168. 2 43 DO 200 JV = 1,IP 169. 3 44 XMAX(JV) = X(1,JV) 170. 3 45 XMIN(JV) = X(1,JV) 171. 3 46 200 CONTINUE 172. 2 47 300 CONTINUE 173. 2 48 DO 400 JV = 1,IP 174. 3 49 IF (X(I,JV) .LT. XMIN(JV)) XMIN(JV)=X(I,JV) 175. 3 51 IF (X(I,JV) .GT. XMAX(JV)) XMAX(JV)=X(I,JV) 176. 3 53 400 CONTINUE 177. 2 54 500 CONTINUE 178. C WRITE MINIMA AND MAXIMA: 179. 1 55 WRITE (6,64000) 180. 1 56 WRITE (6,DATFMT) (XMIN(JV),JV=1,IP) 181. 1 57 WRITE (6,66000) 182. 1 58 WRITE (6,DATFMT) (XMAX(JV),JV=1,IP) 183. C 184. C 188. C 216. C COMPUTE SUM VECTOR AND SUM-OF-PRODUCTS MATRIX: 216.1 1 59 DO 105 JV = 1,IP 216.2 2 60 TOTAL(JV) = 0.0D+00 216.3 2 61 DO 105 JW = 1,IP 216.4 3 62 SUMSQS(JV,JW) = 0.0D+00 216.5 3 63 SSDEVS(JV,JW) = 0.0D+00 216.51 3 64 105 CONTINUE 216.6 1 65 DO 110 JV = 1,IP 216.7 2 66 DO 110 I = 1,N 216.8 3 67 TOTAL(JV) = TOTAL(JV) + X(I,JV) 216.9 3 68 110 CONTINUE 217. C 217.01 1 69 DO 130 JV = 1,IP 217.02 2 70 XMEAN(JV) = TOTAL(JV)/XN 217.03 2 71 130 CONTINUE 217.04 1 72 DO 120 JV = 1,IP 217.1 2 73 DO 120 JW = 1,IP 217.2 3 74 DO 120 I = 1,N 217.3 4 75 DEVV=X(I,JV)-XMEAN(JV) 217.31 4 76 DEVW=X(I,JW)-XMEAN(JW) 217.32 4 77 SSDEVS(JV,JW) = SSDEVS(JV,JW) + DEVV*DEVW 217.4 4 78 120 CONTINUE 217.5 1 79 DO 140 JV = 1,IP 217.9 2 80 DO 140 JW = 1,IP 218. 3 81 VARHAT(JV,JW) = SSDEVS(JV,JW)/XN 218.2 3 82 SIGMA(ISAMPL,JV,JW) = SSDEVS(JV,JW) 218.21 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 5 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 3 83 140 CONTINUE 218.3 C WRITE SUMMARY STATISTICS FOR THIS SAMPLE: 218.4 1 84 WRITE(6,21000) 218.6 1 85 WRITE(6,22000) 218.7 1 86 WRITE(6,10100) ( XMEAN(JV), JV = 1,IP ) 218.8 1 87 WRITE(6,10200) 218.9 1 88 DO 150 JV = 1,IP 219. 2 89 WRITE(6,48000) ( VARHAT(JV,JW), JW=1,IP ) 219.1 2 90 150 CONTINUE 219.2 C 219.3 C 219.4 C CALL SUBROUTINE TO COMPUTE DETERMINANT OF SAMPLE COVARIANCE 219.5 C MATRIX: 219.6 C 219.7 1 91 IDET = 1 219.8 1 92 NRS1 = 0 219.9 1 93 DO 2404 JV = 1,IP 220. 2 94 DO 2404 JW = 1,IP 220.1 3 95 COVMX(JV,JW) = VARHAT(JV,JW) 220.2 3 96 2404 CONTINUE 220.3 1 97 WRITE (6,11500) 220.31 1 98 WRITE (6,11700) COVMX(1,1) 220.32 1 99 CALL MATEQ(COVMX,IP,20,JFLG,DET,IDET,IV,NRS1,P,20) 220.4 1 100 WRITE (6,11600) 220.41 C 220.5 C GENERAL FORM OF CALL IS: 220.6 C CALL MATEQ(A,M,N,JFLG,DET,IDET,IV,NRS1,P,LL) 220.7 C SEE SUBROUTINE LISTING FOR FULLER EXPLANATION. 220.8 C DET(A) = DET*10.0**IDET 220.9 C 221. 1 101 WRITE (6,60000) JFLG 222. C 223. C 224. 1 102 XIDET = IDET 225. C XLGDET = DLOG(DET) + XIDET*ALOG(10.0) 226. 1 103 TRUDET = DET*(10.0**IDET) 227. 1 104 WRITE (6,62000) DET,IDET 228. C 229. C COMPUTE LIKELIHOOD ACCORDING TO 230. C MAX LIKELIHOOD = 231. C (2*PI)**(-N*P/2)*(DET OF COV MX)**(-N/2)*EXP(-N*P/2) 232. C COMPUTE LOG OF MAX LIKELIHOOD: 233. 1 105 XLL=-(XN*IP/2.0)*DLOG(2.0D+00*PI)-(XN/2.0)*DLOG(TRUDET)-XN*IP/2.0 234. 1 106 XLLIKE(ISAMPL) = XLL 235. C 236. C 237. C 238. C GO ON TO NEXT SAMPLE (I.E., NEXT VALUE OF ISAMPL): 464. C 465. 1 107 4100 CONTINUE 466. C 467. C EVALUATE LIKELIHOOD WITH VARYING COVARIANCE MATRICES BY 468. C MULTIPLYING LIKELIHOODS FOR THE SEPARATE SAMPLES, 469. C I.E., BY ADDING THEIR LOGARITHMS: 469.1 108 XLOGL = 0.0 470. 109 DO 4150 ISAMPL=1,K 471. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 6 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 1 110 XLOGL = XLOGL + XLLIKE(ISAMPL) 472. 1 111 4150 CONTINUE 473. C COMPUTE MODEL-SELECTION CRITERIA: 474. 112 WRITE (6,69000) 475. C COMPUTE MODEL-SELECTION CRITERIA FOR MODEL WITH 476. C VARYING COVARIANCE MATRICES:-- 477. 113 WRITE (6,69100) 478. C 479. C PARAMETERS: 480. C K MEAN VECTORS OF DIMENSION P AND K P-BY-P COVARIANCE MATRICES, 481. C WHERE P IS THE NUMBER OF VARIABLES:-- 482. C 483. 114 NOPARM = K*IP + K*IP*(IP+1)/2 484. C 484.1 115 WRITE (6,72000) NOPARM 485. 116 XMN2LL = -2.0*XLOGL 486. 117 WRITE (6,76000) XMN2LL 486.1 118 AIC = XMN2LL + 2.0*NOPARM 487. 119 WRITE (6,78000) AIC 488. 120 SCH = XMN2LL + ALOG(XN)*NOPARM 489. 121 WRITE (6,80000) SCH 490. C COMPUTE MODEL-SELECTION CRITERIA FOR MODEL WITH 491. C COMMON COVARIANCE MATRIX:-- 492. C 493. C BEGIN BY COMPUTING POOLED ESTIMATE OF COVARIANCE MATRIX: 494. C 494.1 C POOL: 494.2 C 494.3 122 DO 1950 JV = 1,IP 494.4 1 123 DO 1950 JW = 1,IP 494.5 2 124 WGSS(JV,JW) = 0.0D+00 494.6 2 125 1950 CONTINUE 494.7 C 494.8 126 DO 2300 JV = 1,IP 494.9 1 127 DO 2300 JW = 1,IP 495. 2 128 DO 2300 IC = 1,K 495.1 3 129 WGSS(JV,JW) = WGSS(JV,JW) + SIGMA(IC,JV,JW) 495.2 3 130 2300 CONTINUE 495.3 C 495.301 131 N = 0 495.31 132 DO 4160 ISAMPL=1,K 495.32 1 133 N = N + NG(ISAMPL) 495.33 1 134 4160 CONTINUE 495.34 135 XN = N 495.35 C 495.36 C COMPUTE COMSIG, MLE OF COMMON COVARIANCE MATRIX: 495.4 136 DO 2400 JV = 1,IP 495.5 1 137 DO 2400 JW = 1,IP 495.6 2 138 COMSIG(JV,JW) = WGSS(JV,JW)/XN 495.7 2 139 2400 CONTINUE 495.8 C 495.9 140 WRITE (6,48500) 495.901 141 DO 2410 JV = 1,IP 495.91 1 142 WRITE (6,48000) ( COMSIG(JV,JW),JW=1,IP ) 495.92 1 143 2410 CONTINUE 495.94 C 495.95 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 7 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C 496. C 496.1 C CALL SUBROUTINE TO COMPUTE DETERMINANT OF SAMPLE COVARIANCE 496.2 C MATRIX: 496.3 C 496.4 144 IDET = 1 496.5 145 NRS1 = 0 496.6 146 DO 2405 JV = 1,IP 496.7 1 147 DO 2405 JW = 1,IP 496.8 2 148 COVMX(JV,JW) = COMSIG(JV,JW) 496.9 2 149 2405 CONTINUE 497. 150 CALL MATEQ(COVMX, IP,20,JFLG,DET,IDET,IV,NRS1,P,20) 497.1 C 497.2 C GENERAL FORM OF CALL IS: 497.3 C CALL MATEQ(A,M,N,JFLG,DET,IDET,IV,NRS1,P,LL) 497.4 C SEE SUBROUTINE LISTING FOR FULLER EXPLANATION. 497.5 C DET(VARHAT) = DET*10.0**IDET 497.6 C 497.7 151 WRITE (6,60000) JFLG 497.8 C 497.9 C 498. 152 XIDET = IDET 498.1 C XLGDET = DLOG(DET) + XIDET*ALOG(10.0) 498.2 153 TRUDET = DET*(10.0**IDET) 498.3 154 IF (JFLG .GT. 0) WRITE (6,62000) DET,IDET 498.4 C 498.5 C COMPUTE LIKELIHOOD OF MODEL WITH COMMON COVARIANCE MATRIX 498.6 C ACCORDING TO 498.7 C MAX LIKELIHOOD = 498.8 C (2*PI)**(-N*P/2)*(DET OF COV MX)**(-N/2)*EXP(-N*P/2) 498.9 C COMPUTE LOG OF MAX LIKELIHOOD: 499. 156 XLL=-(XN*IP/2.0)*DLOG(2.0D+00*PI)-(XN/2.0)*DLOG(TRUDET)-XN*IP/2.0 499.1 C 499.3 157 XMN2LL = -2.0*XLL 499.4 C 499.41 158 WRITE (6,68000) 499.5 159 NOPARM = K*IP + IP*(IP+1)/2 499.6 160 WRITE (6,72000) NOPARM 499.7 161 WRITE (6,76000) XMN2LL 499.71 162 AIC = XMN2LL + 2.0*NOPARM 499.8 163 WRITE (6,78000) AIC 499.9 164 SCH = XMN2LL + ALOG(XN)*NOPARM 501. 165 WRITE (6,80000) SCH 502. C 503. C 518. 166 WRITE (6,84000) 519. C 520. C 521. 167 STOP 522. C 523. C 524. C 525. 168 10000 FORMAT('1','**********************************************'/// 526. X' EVALUATION OF LIKELIHOOD IN ONE-WAY CLASSIFICATION '/ 527. X' WITH COMMON AND VARYING COVARIANCE MATRICES '// 528. X' CMS FILE: EVALUATE CLUSPAC '/// 529. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 8 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 X' DEVELOPED AND PROGRAMMED BY: '// 530. X' DR. STANLEY L. SCLOVE '/ 531. X' INFORMATION & DECISION SCIENCES DEPT. M/C 294 '/ 532. X' COLLEGE OF BUSINESS ADMINISTRATION '/ 533. X' UNIVERSITY OF ILLINOIS AT CHICAGO '/ 534. X' BOX 4348, CHICAGO, IL 60680 '// 535. X' 312/996-2681 '// 536. X' EVALUATE CLUSPAC VERSION 1.4 30-JAN-90 '// 537. X' COPYRIGHT (C) 1990 STANLEY L. SCLOVE. ', 538. X' ALL RIGHTS RESERVED.'//) 539. 169 10050 FORMAT(/' ................................................', 539.7 X'................................'//) 539.8 170 11000 FORMAT(2X,I2) 540. 171 11500 FORMAT(2X,'ENTER SUBROUTINES.') 540.1 172 11700 FORMAT(2X,'(1,1) ELEMENT OF INPUT MATRIX = ',E15.7) 540.2 173 11600 FORMAT(2X,'LEAVE SUBROUTINES.') 540.3 174 12000 FORMAT(2X,I4) 541. 175 13000 FORMAT(1X,'NUMBER OF VARIABLES USED .................',I3/) 542. 176 26000 FORMAT(1X,'NUMBER OF GROUPS, K ......................',I3/) 542.1 177 15000 FORMAT(1X,'NUMBER OF OBSERVATIONS (SAMPLE SIZE), N...',I3/) 543. 178 18000 FORMAT(15(I5,I3)) 546. 179 20000 FORMAT(/1X,' MEANS: '/) 547. 180 25000 FORMAT(5X, F3.2) 550. 181 21000 FORMAT(/ ' STATISTICS FOR THIS SAMPLE: ' ) 552. 182 22000 FORMAT( ' ---------- --- ---- ------ ' /) 553. 183 10100 FORMAT(' MEAN VECTOR: ',(9F8.2/)//) 553.1 184 10200 FORMAT(/' COVARIANCE MATRIX: '/) 553.2 185 28000 FORMAT(1X,'MEAN VECTOR FOR GROUP ',I2,': ',(8F13.5/)) 553.3 186 30000 FORMAT(/1X,' MINUS 2 LOG LIKELIHOOD = ', F13.5/) 554. 187 36000 FORMAT(18A4) 558. 188 40000 FORMAT(//,1X,18A4//) 559. 189 42000 FORMAT(/,' COVARIANCE MATRICES: ',/) 560. 190 42500 FORMAT(1X,' GROUP ',I2/) 561. 191 48000 FORMAT(1X,8F13.5/) 562. 192 48500 FORMAT(1X,'MLE OF COMMON COVARIANCE MATRIX: '/) 562.1 193 50000 FORMAT(1X,I4,1X,I2) 563. 194 54000 FORMAT(3X,I2) 565. 195 60000 FORMAT(/,1X,'JFLG = ',I2,'. IF JFLG=0, COMPUTATION OF DET', 568. X' WENT WELL; OTHERWISE, THERE WAS TROUBLE OR MATRIX WAS ', 569. X'ILL-CONDITIONED.'//) 570. 196 62000 FORMAT(/1X,'DET = ',E18.7,' IDET = ',I3,5X, 571. X'ACTUAL DET. = DET*10**IDET',//) 572. 197 64000 FORMAT(/1X,'MINIMUM FOR EACH VARIABLE: ',/) 573. 198 66000 FORMAT(/1X,'MAXIMUM FOR EACH VARIABLE: ',/) 574. 199 69000 FORMAT(1X,'VALUES OF MODEL SELECTION CRITERIA: '/) 576. 200 69100 FORMAT(1X,'FOR MODEL WITH VARYING COVARIANCE MATRICES:--') 577. 201 68000 FORMAT(1X,'FOR MODEL WITH COMMON COVARIANCE MATRIX:--') 578. 202 70000 FORMAT(1X,'RESULTS FOR SAMPLE ',I2/) 579. 203 72000 FORMAT(1X,'NUMBER OF PARAMETERS = ',5X,I4/) 582. 204 76000 FORMAT(1X,' MINUS 2 LOG L = ', F15.5 /) 582.1 205 78000 FORMAT(1X,' AIC = ', F15.5 /) 582.3 206 80000 FORMAT(1X,' SCHWARZ CRITERION = ', F15.5 /) 582.4 207 84000 FORMAT(1X,'PROGRAM ENDED NORMALLY.') 584. 208 END 585. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 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 0NAME TYPE TAG DECLARED REFERENCED +______ ____ ______ ________ ________________________________________ 0AIC R*4 118 119 162 163 ALOG R*4 I 120 164 COMSIG R*4 A 7 138 142 148 COVMX R*8 AT 7 13 95 98 99 148 150 DATFMT R*4 A 3 30 39 56 58 DET R*8 T 15 99 103 104 150 153 155 DEVV R*8 T 17 75 77 DEVW R*8 T 17 76 77 DLOG R*8 I 105 105 156 156 DSQ R*8 T 16 UNREFERENCED F R*8 T 18 UNREFERENCED I I*4 38 39 40 49 50 51 52 66 67 74 75 76 IC I*4 128 129 IDET I*4 91 99 102 103 104 144 150 152 153 155 IP I*4 27 28 39 43 48 56 58 59 61 65 69 72 73 79 80 86 88 89 93 94 99 105 105 114 114 114 122 123 126 127 136 137 141 142 146 147 150 156 156 159 159 159 ISAMPL I*4 32 33 36 82 106 109 110 132 133 IV I*4 A 9 99 150 JFLG I*4 99 101 150 151 154 JV I*4 39 39 43 44 44 45 45 48 49 49 50 50 51 51 52 52 56 56 58 58 59 60 62 63 65 67 67 67 69 70 70 72 75 75 77 77 79 81 81 82 82 86 86 88 89 93 95 95 122 124 126 129 129 129 136 138 138 141 142 146 148 148 JW I*4 61 62 63 73 76 76 77 77 80 81 81 82 82 89 89 94 95 95 123 124 127 129 129 129 137 138 138 142 142 147 148 148 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 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 0NAME TYPE TAG DECLARED REFERENCED +______ ____ ______ ________ ________________________________________ 0K I*4 25 26 32 109 114 114 128 132 159 MATEQ X 99 150 N I*4 34 35 36 37 38 66 74 131 133 133 135 NG I*4 A 4 36 133 NOPARM I*4 114 115 118 120 159 160 162 164 NRS1 I*4 92 99 145 150 P R*8 T 14 99 150 PI R*4 20 105 156 SCH R*4 120 121 164 165 SIGMA R*4 A 5 82 129 SSD R*8 AT 5 12 SSDEVS R*8 AT 5 13 63 77 77 81 82 SUM R*8 AT 1 11 SUMSQS R*8 AT 4 13 62 TITLE R*4 A 2 21 24 TOTAL R*8 AT 4 13 60 67 67 70 TRUDET R*8 T 15 103 105 153 156 VARHAT R*8 AT 7 13 81 89 95 WGSS R*8 AT 6 12 124 129 129 138 X R*4 A 1 39 44 45 49 50 51 52 67 75 76 XIDET R*4 102 152 XLL R*4 105 106 156 157 XLLIKE R*8 AT 10 19 106 110 XLOGL R*4 108 110 110 116 XMAX R*4 A 8 44 51 52 58 XMEAN R*4 A 1 70 75 76 86 XMIN R*4 A 8 45 49 50 56 XMN2LL R*4 116 117 118 120 157 161 162 164 XN R*4 35 70 81 105 105 105 120 135 138 156 156 156 164 XUFLOW X 31 0 0VARIABLES REFERENCED BUT NOT SET. (* POSSIBLY SET AS ARGUMENT.) 0DET* IV* JFLG* P* 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 11 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 42 40 105 64 59 61 110 68 65 66 120 78 72 73 74 130 71 69 140 83 79 80 150 90 88 200 46 43 300 B 47 41 400 53 48 500 54 38 1950 125 122 123 2300 130 126 127 128 2400 139 136 137 2404 96 93 94 2405 149 146 147 2410 143 141 4100 107 32 4150 111 109 4160 134 132 10000 NF 168 22 10050 NF 169 23 29 10100 NF 183 86 10200 NF 184 87 11000 NF 170 25 11500 NF 171 97 11600 NF 173 100 11700 NF 172 98 12000 NF 174 34 13000 NF 175 28 15000 NF 177 37 18000 NF 178 UNREFERENCED 20000 NF 179 UNREFERENCED 21000 NF 181 84 22000 NF 182 85 25000 NF 180 UNREFERENCED 26000 NF 176 26 28000 NF 185 UNREFERENCED 30000 NF 186 UNREFERENCED 36000 NF 187 21 30 40000 NF 188 24 42000 NF 189 UNREFERENCED 42500 NF 190 UNREFERENCED 48000 NF 191 89 142 48500 NF 192 140 50000 NF 193 UNREFERENCED 54000 NF 194 27 60000 NF 195 101 151 62000 NF 196 104 155 64000 NF 197 55 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:31 NAME:MAIN PAGE: 12 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 66000 NF 198 57 68000 NF 201 158 69000 NF 199 112 69100 NF 200 113 70000 NF 202 33 72000 NF 203 115 160 76000 NF 204 117 161 78000 NF 205 119 163 80000 NF 206 121 165 84000 NF 207 166 0*STATISTICS* SOURCE STATEMENTS = 205, PROGRAM SIZE = 156836 BYTES, PROGRAM NAME = MAIN PAGE: 1. 0*STATISTICS* NO DIAGNOSTICS GENERATED. 0**MAIN** END OF COMPILATION 1 ****** TIME STAMP: 90.03918.58.31 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 PAGE: 13 0OPTIONS IN EFFECT: NOLIST NOMAP XREF NOGOSTMT NODECK SOURCE TERM OBJECT FIXED NOTRMFLG SRCFLG NOSYM NORENT SDUMP(ISN) SXM NOVECTOR IL(DIM) NOTEST NODC NOICA NODIRECTIVE NODBCS NOSAA OPT(0) LANGLVL(77) NOFIPS FLAG(I) AUTODBL(NONE) NAME(MAIN) LINECOUNT(60) CHARLEN(500) 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C 586. C 587. 1 SUBROUTINE MATEQ(A,M,N,JFLG,DET,IDET,IV,NRS1,P,LL) 588. C SUBROUTINE MATEQ IS DMATEQ FROM THE UICC SUBROUTINE LIBRARY. 589. C PAGE 1 590. C 591. C SUBROUTINE DMATEQ 592. C ***************** 593. C THIS ROUTINE WILL SOLVE A REAL*8 SYSTEM OF LINEAR EQUATIONS,COMPUTE 594. C THE DETERMINANT, WITHOUT UNDERFLOW OR OVERFLOW, OF A REAL*8 MATRIX, 595. C AND/OR INVERT A REAL*8 MATRIX. 596. C CALLING SEQUENCE: 597. C CALL DMATEQ(A,N,IA,JFLG,DET,IDET,IV,NRS,P,IP) WHERE; 598. C A (INPUT) - IS THE REAL*8 MATRIX ON WHICH THE ROUTINE IS 599. C TO WORK. IN THE PROCESS OF COMPUTATION THE 600. C CONTENTS OF THIS MATRIX ARE DESTROYED. 601. C N (INPUT) - IS AN INTEGER*4 VARIABLE WHICH SPECIFIES THE 602. C ORDER OF THE A MATRIX. 603. C IA (INPUT) - IS AN INTEGER*4 VARIABLE WHICH SPECIFIES THE 604. C ACTUAL ROW DIMENSION OF A AS DIMENSIONED IN 605. C THE CALLING PROGRAM. IA MUST BE GREATER THAN 606. C OR EQUAL TO N. 607. C JFLG (OUTPUT) - IS AN INTEGER*4 RETURN CODE VARIABLE. UPON 608. C RETURN FROM DMATEQ IF; 609. C JFLG=0, ALL WENT WELL. 610. C JFLG=1, THE A MATRIX WAS SINGULAR OR NEAR 611. C SINGULAR AND THE COMPUTATIONS COULD NOT BE 612. C COMPLETED. THE CONTENTS OF THE VARIABLES 613. C A, DET, IDET AND P ARE MEANINGLESS. 614. C DET (OUTPUT) - IS A REAL*8 VARIABLE WHICH CONTAINS THE 615. C DETERMINANT OF A. (SEE IDET) 616. C IDET (INPUT) - IS AN INTEGER*4 VARIABLE. ON INPUT IF; 617. C IDET=0, NO DETERMINANT IS CALCULATED. 618. C IDET NOT 0, THE DETERMINANT OF A IS COMPUTED. 619. C ON OUTPUT IDET CONTAINS THE POWER OF 10 620. C THAT DET SHOULD BE MULTIPLIED BY TO GIVE THE 621. C CORRECT VALUE OF THE DETERMINANT. I.E. 622. C DET(A)=DET*10.0D0**IDET. 623. C IF DET(A) CAN BE COMPUTED WITHOUT UNDER OR 624. C OVERFLOW, THEN IDET=0 OTHERWISE IDET IS SET 625. C TO THE PROPER VALUE SO THAT NO UNDER OR OVER- 626. C FLOW WILL OCCUR IN COMPUTING DET. 627. C IV (INPUT) - IS AN INTEGER*4 WORK ARRAY WHICH SHOULD BE 628. C DIMENSIONED AT LEAST IV(N). 629. C PAGE 2 630. C 631. C NRS (INPUT) - IS AN INTEGER*4 VARIABLE WITH THE FOLLOWING 632. C INTERPRETATION: 633. C NRS>0, SOLVE A SYSTEM OF LINEAR EQUATIONS 634. C WITH NRS RIGHT HAND SIDES. 635. C NRS=0, INVERT THE A MATRIX. 636. C NRS<0, ONLY COMPUTE THE DETERMINANT OF A. 637. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MAIN PAGE: 14 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C IN THIS CASE IDET MUST BE DIFFERENT FROM 0. 638. C P (INPUT) - IS A REAL*8 ARRAY WITH THE FOLLOWING INTER- 639. C PRETATION: 640. C IF NRS>0, THEN P CONTAINS THE NRS RIGHT HAND 641. C SIDES STORED BY COLUMNS. IN THIS CASE P MUST 642. C BE DIMENSIONED AT LEAST P(N,NRS). ON RETURN 643. C THE COLUMNS OF P ARE REPLACED BY THE RESPEC- 644. C TIVE SOLUTIONS. 645. C IF NRS=0, THEN P MUST BE DIMENSIONED AT LEAST 646. C P(N,N). ON RETURN P WILL CONTAIN THE INVERSE 647. C OF A. 648. C IF NRS<0,THEN P NEED ONLY BE A DUMMY VARIABLE 649. C IN THIS CASE P IS NEVER ACCESSED BY DMATEQ. 650. C IP (INPUT) - IS AN INTEGER*4 VARIABLE WHICH CONTAINS THE 651. C ACTUAL ROW DIMENSION OF P AS DIMENSIONED IN 652. C THE CALLING PROGRAM. IP MUST BE GREATER THAN 653. C OR EQUAL TO N. 654. C NOTE: IMMEDIATELY ON RETURN FROM DMATEQ THE CONDITION CODE FLAG, 655. C JFLG, SHOULD BE INTERROGATED. IF JFLG=1, THEN THE ROUTINE 656. C COULD NOT COMPUTE A SOLUTION. 657. C METHOD - THE ALGORITHM USED IS GAUSSIAN ELIMINATION WITH PARTIAL 658. C -1 659. C PIVOTING. IN ESSENCE THE ROUTINE GENERATES A MATRIX L SUCH 660. C -1 661. C THAT L *A = U, WHERE U IS AN UPPER TRIANGULAR MATRIX. THEN IT 662. C SOLVES THE SYSTEM A*X = P BY MEANS OF THE EQUIVALENT SYSTEM 663. C -1 -1 664. C U*X = L *A*X = L *P BY BACK SUBSTITUTION. 665. C -1 666. C THE L MATRIX CAN BE WRITTEN AS A PRODUCT OF THE FORM 667. C -1 668. C L = L *P *....*L *P WHERE EACH P IS A PERMUTATION 669. C N-1 N-1 1 1 K 670. C MATRIX OBTAINED BY INTERCHANGING AT MOST TWO ROWS OF THE 671. C IDENTITY MATRIX. ( THIS REPRESENTS THE INTERCHANGING OF TWO 672. C ROWS). THE L MATRICES ARE ELIMINATION MATRICES WHICH ARE 673. C K 674. C CHOSEN TO INTRODUCE ZEROS IN THE LAST N-K ENTRIES OF THE K-TH 675. C COLUMN OF THE MATRIX. 676. C PAGE 3 677. C 678. C -1 -1 679. C THE CALCULATIONS OF L *A AND L *P ARE DONE BY PERFORMING 680. C THE PERMUTATIONS ON A AND P RESPECTIVELY. THE ACTUAL L AND P 681. C K K 682. C ARE NOT COMPUTED. 683. C SUBROUTINES CALLED: DMATDT 684. C REFERENCE: 685. C G. W. STEWART, INTRODUCTION TO MATRIX COMPUTATIONS, 686. C ACADEMIC PRESS, 1973. 687. 2 REAL*8 A(N,1),DET,P(LL,1) 688. 3 REAL*8 DNORM,DEN,DMULT,DSUM,DISIGN 689. 4 DIMENSION IV(1) 690. 5 NRS=NRS1 691. 6 IF (NRS.EQ.0) IDET=1 692. 8 DISIGN=1.0D+00 693. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATEQ PAGE: 15 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 9 DET=0.0D+00 694. 10 JFLG=0 695. C 696. C JFLG IS A TROUBLE FLAG.UPON EXIT IF JFLG=0 THEN THE MATRIX WAS PROCESS 697. C WITHOUT TROUBLE.IF JFLG=1 EITHER THE MATRIX IS SINGULAR OR TROUBLE 698. C OCCURED.ISIGN=-ISIGN EVERY TIME A ROW IS INTERCHANGED.THIS IS USED TO 699. C INSURE THAT THE DETERMINANT HAS THE PROPER SIGN. 700. C 701. 11 M1=M-1 702. 12 DO 100 I=1,M 703. 1 13 100 IV(I)=I 704. 14 IF (NRS) 500,200,500 705. 15 200 DO 300 I=1,M 706. 1 16 DO 300 J=1,M 707. 2 17 300 P(I,J)=0.0D+00 708. 18 DO 400 I=1,M 709. 1 19 400 P(I,I)=1.0D+00 710. 20 NRS=M 711. C 712. C INSTEAD OF ACTUALLY INTERCHANGING ROWS A POINTER ARRAY IS USED TO KEEP 713. C TRACK OF THE ROW POSITIONS. 714. C 715. C BEGIN ELIMINATION LOOP. 716. C 717. 21 500 DO 1200 K=1,M1 718. 1 22 ICOL=K 719. 1 23 IPCOL=K 720. C 721. C SEARCHING FOR LARGEST ELEMENT IN ABSOLUTE VALUE IN COLUMN K. 722. C 723. 1 24 DNORM=A(IV(K),K) 724. 1 25 IFLG=0 725. 1 26 KK=K+1 726. 1 27 DO 600 J=KK,M 727. 2 28 IF (DABS(A(IV(J),K)).LE.DABS(DNORM)) GO TO 600 728. 2 29 IFLG=1 729. 2 30 IPCOL=IV(J) 730. 2 31 DNORM=A(IPCOL,K) 731. 2 32 600 CONTINUE 732. C 733. C IF IFLG=0 NO ROW INTERCHANGE TOOK PLACE.IF IFLG=1 A ROW INTERCHANGE 734. C TOOK PLACE AND THE POINTER ARRAY IV MUST BE UPDATED. 735. C 736. 1 33 IF (IFLG.EQ.0) GO TO 800 737. 1 34 ISAVE=IV(ICOL) 738. 1 35 IV(ICOL)=IPCOL 739. 1 36 ICOL1=ICOL+1 740. 1 37 DO 700 L=ICOL1,M 741. 2 38 IF (IV(L).EQ.IPCOL) IV(L)=ISAVE 742. 2 40 700 CONTINUE 743. 1 41 DISIGN=-DISIGN 744. 1 42 800 IF (DNORM.EQ.0.0D+00) GO TO 1900 745. C 746. C BEGIN ELIMINATION OF ROW BELOW IV(K).DEN IS THE PIVOT ELEMENT. 747. C 748. 1 43 K1=K+1 749. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATEQ PAGE: 16 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 1 44 DO 1100 IM=K1,M 750. C 751. C BEFORE ACTUALLY ELIMINATING WE CHECK TO SEE IF A(IV(IM),K) HAS 752. C ALREADY BEEN ANIHILATED. 753. C 754. 2 45 IF (A(IV(IM),K).EQ.0.0D+00) GO TO 1100 755. C 756. C CACULATE ELIMINATION FACTOR. 757. C 758. 2 46 DMULT=-A(IV(IM),K) 759. C 760. C WE NOW CALCULATE VALUE OF OTHER ELEMENTS IN ROW IV(IM). 761. C 762. 2 47 DO 900 NN=K1,M 763. 3 48 900 A(IV(IM),NN)=(DMULT*A(IV(K),NN))/DNORM+A(IV(IM),NN) 764. 2 49 IF (NRS.LE.0) GO TO 1100 765. 2 50 DO 1000 IN=1,NRS 766. 3 51 1000 P(IV(IM),IN)=(DMULT*P(IV(K),IN))/DNORM+P(IV(IM),IN) 767. 2 52 1100 CONTINUE 768. 1 53 1200 CONTINUE 769. C 770. C CALCULATE VALUE OF DETERMINANT. 771. C 772. 54 IF (A(IV(M),M).EQ.0.0D0) GO TO 1900 773. 55 DET=DISIGN 774. 56 IF (IDET.NE.0) CALL DMATDT(A,N,M,DET,IV,IDET) 775. 58 IF (DET.EQ.0.0D+00) GO TO 1900 776. 59 IF (NRS.LE.0) GO TO 2000 777. C 778. C WE START SOLVING RIGHT HAND SIDES.THE SOLUTION REPLACES THE RIGHT HAND 779. C VECTOR. 780. C 781. 60 1300 N1=M-1 782. 61 DO 1600 JJ=1,NRS 783. C 784. C BEGIN BACK SUBSTITUTION. 785. C 786. 1 62 P(IV(M),JJ)=P(IV(M),JJ)/A(IV(M),M) 787. 1 63 DO 1500 I=1,N1 788. 2 64 DSUM=0.0D+00 789. 2 65 DO 1400 J=1,I 790. 3 66 1400 DSUM=DSUM-A(IV(M-I),M-J+1)*P(IV(M-J+1),JJ) 791. 2 67 1500 P(IV(M-I),JJ)=(P(IV(M-I),JJ)+DSUM)/A(IV(M-I),M-I) 792. 1 68 1600 CONTINUE 793. 69 DO 1800 JJ=1,NRS 794. 1 70 DO 1700 IND=1,M 795. 2 71 1700 A(IND,1)=P(IV(IND),JJ) 796. 1 72 DO 1800 IND=1,M 797. 2 73 1800 P(IND,JJ)=A(IND,1) 798. 74 RETURN 799. 75 1900 JFLG=1 800. 76 IDET=0 801. 77 2000 RETURN 802. 78 END 803. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATEQ PAGE: 17 0SYMBOL CROSS REFERENCE DICTIONARY 0PROGRAM NAME: MATEQ 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 0NAME TYPE TAG DECLARED REFERENCED +______ ____ ______ ________ ________________________________________ 0A R*8 ADT 1 2 24 28 31 45 46 48 48 48 54 57 62 66 67 71 73 DABS R*8 I 28 28 DEN R*8 T 3 UNREFERENCED DET R*8 DT 1 2 9 55 57 58 DISIGN R*8 T 3 8 41 41 55 DMATDT X 57 DMULT R*8 T 3 46 48 51 DNORM R*8 T 3 24 28 31 42 48 51 DSUM R*8 T 3 64 66 66 67 I I*4 12 13 13 15 17 18 19 19 63 65 66 67 67 67 67 ICOL I*4 22 34 35 36 ICOL1 I*4 36 37 IDET I*4 D 1 7 56 57 76 IFLG I*4 25 29 33 IM I*4 44 45 46 48 48 51 51 IN I*4 50 51 51 51 IND I*4 70 71 71 72 73 73 IPCOL I*4 23 30 31 35 38 ISAVE I*4 34 39 IV I*4 AD 1 4 13 24 28 30 34 35 38 39 45 46 48 48 48 51 51 51 54 57 62 62 62 66 66 67 67 67 71 J I*4 16 17 27 28 30 65 66 66 JFLG I*4 D 1 10 75 JJ I*4 61 62 62 66 67 67 69 71 73 K I*4 21 22 23 24 24 26 28 31 43 45 46 48 51 KK I*4 26 27 K1 I*4 43 44 47 L I*4 37 38 39 LL I*4 D 1 2 M I*4 D 1 11 12 15 16 18 20 27 37 44 47 54 54 57 60 62 62 62 62 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATEQ PAGE: 18 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 0NAME TYPE TAG DECLARED REFERENCED +______ ____ ______ ________ ________________________________________ 0 66 66 66 67 67 67 67 70 72 MATEQ R 1 M1 I*4 11 21 N I*4 D 1 2 57 NN I*4 47 48 48 48 NRS I*4 5 6 14 20 49 50 59 61 69 NRS1 I*4 D 1 5 N1 I*4 60 63 P R*8 ADT 1 2 17 19 51 51 51 62 62 66 67 67 71 73 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 13 12 200 B 15 14 300 17 15 16 400 19 18 500 B 21 14 14 600 B 32 27 28 700 40 37 800 B 42 33 900 48 47 1000 51 50 1100 B 52 44 45 49 1200 53 21 1300 60 UNREFERENCED 1400 66 65 1500 67 63 1600 68 61 1700 71 70 1800 73 69 72 1900 B 75 42 54 58 2000 B 77 59 0*STATISTICS* SOURCE STATEMENTS = 75, PROGRAM SIZE = 4792 BYTES, PROGRAM NAME = MATEQ PAGE: 13. 0*STATISTICS* NO DIAGNOSTICS GENERATED. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATEQ PAGE: 19 0**MATEQ** END OF COMPILATION 2 ****** TIME STAMP: 90.03918.58.33 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 PAGE: 20 0OPTIONS IN EFFECT: NOLIST NOMAP XREF NOGOSTMT NODECK SOURCE TERM OBJECT FIXED NOTRMFLG SRCFLG NOSYM NORENT SDUMP(ISN) SXM NOVECTOR IL(DIM) NOTEST NODC NOICA NODIRECTIVE NODBCS NOSAA OPT(0) LANGLVL(77) NOFIPS FLAG(I) AUTODBL(NONE) NAME(MAIN) LINECOUNT(60) CHARLEN(500) 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 1 SUBROUTINE MATDT(A,IA,N,DET,IV,IDET) 804. C SUBROUTINE MATDT IS DMATDT FROM THE UICC SUBROUTINE LIBRARY. 805. 2 REAL*8 A(IA,1),DET,B,LOG16 806. 3 INTEGER*4 IV(1),K 807. 4 EQUIVALENCE (B,K) 808. 5 NUM=16777216 809. 6 LOG16=.120411998265592457D+01 810. 7 IF (A(IV(N),N).EQ.0.0D+00) GO TO 300 811. 8 L=0 812. 9 DO 100 I=1,N 813. 1 10 B=DABS(A(IV(I),I)) 814. 1 11 K=K/NUM-64 815. 1 12 L=L+K 816. 1 13 100 DET=DET*(A(IV(I),I)/16.0D+00**K) 817. 14 B=DABS(DET) 818. 15 K=K/NUM-64 819. 16 IW=L+K 820. 17 IF ((IW.LT.-64).OR.(IW.GT.63)) GO TO 200 821. 18 DET=DET*16.0D+00**L 822. 19 IDET=0 823. 20 GO TO 400 824. 21 200 DET=DET*16.0D+00**(-K) 825. 22 IDET=L+K 826. 23 B=IDET*LOG16 827. 24 IDET=B 828. 25 B=B-DFLOAT(IDET) 829. 26 DET=DET*1.0D+01**B 830. 27 GO TO 400 831. 28 300 DET=0.0D+00 832. 29 IDET=0 833. 30 400 RETURN 834. C CONTROL CARDS FOR PROGRAM EVALUATE CLUSPAC: C 835. C C 836. C //GO.SYSIN DD * 837. C FISHER IRIS DATA 838. C K=03 839. C IP=04 840. C (4(1X,F3.1),1X,I1) 841. C N=0050 842. C C 843. C THESE ARE: 844. C JCL "GO"-STEP CARD 845. C DATASET TITLE C 846. C K, NUMBER OF GROUPS, IN FORMAT (2X,I2) C 847. C NO. OF VARIABLES, IP, IN FORMAT (3X,I2) 848. C DATFMT, IN FORMAT (18A4), E.G., (4F4.1) C 849. C "DATFMT" WILL ALSO BE USED FOR OUTPUT, SO ALLOW AT LEAST C 850. C ONE BLANK AT THE BEGINNING FOR CARRIAGE CONTROL. C 851. C N FOR FIRST SAMPLE, IN FORMAT (2X,I4) 852. C DATA FOR FIRST SAMPLE 853. C 854. C (DATA, ONE CASE AT A TIME, IN FORMAT SPECIFIED BY DATFMT) 855. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATDT PAGE: 21 0 IF DO ISN *....*...1.........2.........3.........4.........5.........6.........7.*.......8 0 C 856. C N FOR SECOND SAMPLE, IN FORMAT (2X,I4) 857. C DATA FOR SECOND SAMPLE 858. C . 859. C . 860. C . 861. C 862. C N FOR K-TH SAMPLE, IN FORMAT (2X,I4) 863. C DATA FOR K-TH SAMPLE 864. C 865. 31 END 866. 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:33 NAME:MATDT PAGE: 22 0SYMBOL CROSS REFERENCE DICTIONARY 0PROGRAM NAME: MATDT 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 0NAME TYPE TAG DECLARED REFERENCED +______ ____ ______ ________ ________________________________________ 0A R*8 ADT 1 2 7 10 13 B R*8 ET 2 4 10 14 23 24 25 25 26 DABS R*8 I 10 14 DET R*8 DT 1 2 13 13 14 18 18 21 21 26 26 28 DFLOAT R*8 I 25 I I*4 9 10 10 13 13 IA I*4 D 1 2 IDET I*4 D 1 19 22 23 24 25 29 IV I*4 ADT 1 3 7 10 13 IW I*4 16 17 17 K I*4 ET 3 4 11 11 12 13 15 15 16 21 22 L I*4 8 12 12 16 18 22 LOG16 R*8 T 2 6 23 MATDT R 1 N I*4 D 1 7 7 9 NUM I*4 5 11 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 100 13 9 200 B 21 17 300 B 28 7 400 B 30 20 27 0*STATISTICS* SOURCE STATEMENTS = 31, PROGRAM SIZE = 1820 BYTES, PROGRAM NAME = MATDT PAGE: 20. 0*STATISTICS* NO DIAGNOSTICS GENERATED. 0**MATDT** END OF COMPILATION 3 ****** TIME STAMP: 90.03918.58.33 1LEVEL 2.3.0 (MAR 1988) VS FORTRAN FEB 08, 1990 18:58:34 PAGE: 23 0SUMMARY OF MESSAGES AND STATISTICS FOR ALL COMPILATIONS 0*STATISTICS* SOURCE STATEMENTS = 205, PROGRAM SIZE = 156836 BYTES, PROGRAM NAME = MAIN PAGE: 1. 0*STATISTICS* NO DIAGNOSTICS GENERATED. 0**MAIN** END OF COMPILATION 1 ****** TIME STAMP: 90.03918.58.33 0*STATISTICS* SOURCE STATEMENTS = 75, PROGRAM SIZE = 4792 BYTES, PROGRAM NAME = MATEQ PAGE: 13. 0*STATISTICS* NO DIAGNOSTICS GENERATED. 0**MATEQ** END OF COMPILATION 2 ****** TIME STAMP: 90.03918.58.33 0*STATISTICS* SOURCE STATEMENTS = 31, PROGRAM SIZE = 1820 BYTES, PROGRAM NAME = MATDT PAGE: 20. 0*STATISTICS* NO DIAGNOSTICS GENERATED. 0**MATDT** END OF COMPILATION 3 ****** TIME STAMP: 90.03918.58.33 0******* SUMMARY STATISTICS ******* 0 DIAGNOSTICS GENERATED. HIGHEST SEVERITY CODE IS 0. 0 1 MVS/3.7 LOADER (UICC VERSION OF 4/20/77) 0OPTIONS USED - PRINT,MAP,LET,CALL,NORES,NOTERM,SIZE=2084864,NAME=**GO EP=MAIN - NAME TYPE ADDR NAME TYPE ADDR NAME TYPE ADDR NAME TYPE ADDR NAME TYPE ADDR 0 MAIN SD 117010 MATEQ SD 13D4B8 MATDT SD 13E770 AFBVIOFM*1 SD 13EE90 IFYVIOFM* LR 13EE90 VFRSF# *1 LR 13EEBE VFWSF# *1 LR 13EEE8 VFCSF# * LR 13EF12 VFDSF# * LR 13EF3C VFESF# * LR 13EF66 VFIXF# *1 LR 13EF90 VFSXF# *1 LR 13EFBA VFFXF# *1 LR 13EFE4 AFBFRXPI*1 SD 13F0C0 IFYFRXPI* LR 13F0C0 FRXPI# *1 LR 13F0FA AFBLLGN *1 SD 13F2B8 IH$DLOG * LR 13F2F2 IFYLLGN * LR 13F2F2 D#LOG * LR 13F2F2 DLOG *1 LR 13F2F2 AFBSLGNC*1 SD 13F7A0 IH$ALOG0* LR 13F7DA IFYSLGC * LR 13F7DA A#LOG1 * LR 13F7DA L#LOG1 * LR 13F7DA ALOG10 * LR 13F7DA LOG10 * LR 13F7DA IFYSLGN * LR 13F804 IH$ALOG * LR 13F804 A#LOG * LR 13F804 ALOG *1 LR 13F804 L#OG * LR 13F804 LOG * LR 13F804 AFBVINTE*1 SD 13FB68 IFYVINTE* LR 13FB68 VFEIM# *1 LR 13FB92 VFEIN# * LR 13FBC6 VFEIL# * LR 13FC0C VFEE# *1 LR 13FC52 VFES# * LR 13FC9C VFEP# * LR 13FCD8 VFELC# * LR 13FD14 AFBVXMSK* SD 13FDF8 IFYVXMSK* LR 13FDF8 XUFLOW *1 LR 13FE32 DMATDT *7 SD 13FF08 AFBFDXPI*1 SD 140348 IFYFDXPI* LR 140348 FDXPI# *1 LR 140382 AFBFDXPD*1 SD 140558 IFYFDXPD* LR 140592 FDXPD# *1 LR 140592 AFBVCLSI*1 SD 140FD8 AFBVCLMI*1 SD 141368 AFBVCLMA* LR 1414A8 AFBVEMGI* SD 141530 VFIEMGN#*1 LR 14155E AFBTLGNC* SD 1415E0 LOGTBL# *1 LR 1415E0 AFBVLCI0*1 SD 142600 AFBCLCI0*1 LR 142600 AFBDLCI0*1 LR 142600 AFBVLCI1*1 LR 142600 AFBVLCIN*1 LR 142600 AFBTDXPD* SD 142D28 XYTBL# *1 LR 142D28 0CONCATENATION NUMBER DICTIONARY 0 1 SYS1.VSFORT.VFORTLIB.VER23 MVS004 MEMBER(S) USED - AFBCLCI0 AFBDLCI0 AFBFDXPD AFBFDXPI AFBFRXPI AFBLLGN AFBSLGNC AFBVCLMI AFBVCLSI AFBVINTE AFBVIOFM AFBVLCIN AFBVLCI0 AFBVLCI1 ALOG DLOG FDXPD# FDXPI# FRXPI# LOGTBL# VFEE# VFEIM# VFFXF# VFIEMGN# VFIXF# VFRSF# VFSXF# VFWSF# XUFLOW XYTBL# 0 7 UICCAPP.FORTLIB2 MVS004 MEMBER(S) USED - DMATDT 0 TOTAL LENGTH 2CD18 ENTRY ADDRESS 117010 0 1********************************************** EVALUATION OF LIKELIHOOD IN ONE-WAY CLASSIFICATION WITH COMMON AND VARYING COVARIANCE MATRICES CMS FILE: EVALUATE CLUSPAC DEVELOPED AND PROGRAMMED BY: DR. STANLEY L. SCLOVE INFORMATION & DECISION SCIENCES DEPT. M/C 294 COLLEGE OF BUSINESS ADMINISTRATION UNIVERSITY OF ILLINOIS AT CHICAGO BOX 4348, CHICAGO, IL 60680 312/996-2681 EVALUATE CLUSPAC VERSION 1.4 30-JAN-90 COPYRIGHT (C) 1990 STANLEY L. SCLOVE. ALL RIGHTS RESERVED. ................................................................................ FISHER IRIS DATA NUMBER OF GROUPS, K ...................... 3 NUMBER OF VARIABLES USED ................. 4 ................................................................................ RESULTS FOR SAMPLE 1 NUMBER OF OBSERVATIONS (SAMPLE SIZE), N... 50 MINIMUM FOR EACH VARIABLE: 4.3 2.3 1.0 0.1 MAXIMUM FOR EACH VARIABLE: 5.8 4.4 1.9 0.6 STATISTICS FOR THIS SAMPLE: ---------- --- ---- ------ MEAN VECTOR: 5.01 3.43 1.46 0.25 COVARIANCE MATRIX: 0.12176 0.09723 0.01603 0.01012 0.09723 0.14082 0.01146 0.00911 0.01603 0.01146 0.02956 0.00595 0.01012 0.00911 0.00595 0.01088 ENTER SUBROUTINES. (1,1) ELEMENT OF INPUT MATRIX = 0.1217640E+00 LEAVE SUBROUTINES. JFLG = 0. IF JFLG=0, COMPUTATION OF DET WENT WELL; OTHERWISE, THERE WAS TROUBLE OR MATRIX WAS ILL-CONDITIONED. 0AFB207I VFNTH : PROGRAM INTERRUPT - FLOATING-POINT OVERFLOW EXCEPTION VFNTH : PSW 078D000C4013F14C, REGISTER CONTAINS 2B24EE8C23C040000000000000000000 VFNTH : LAST EXECUTED FORTRAN STATEMENT IN PROGRAM MAIN (OFFSET 02813C). 0TRACEBACK OF CALLING ROUTINES; MODULE ENTRY ADDRESS = 117010. ----------------------------------------------------------------------------- MAIN (117010) CALLED BY OPERATING SYSTEM. ----------------------------------------------------------------------------- 0STANDARD CORRECTIVE ACTION TAKEN. EXECUTION CONTINUING. 0AFB207I VFNTH : PROGRAM INTERRUPT - FLOATING-POINT OVERFLOW EXCEPTION VFNTH : PSW 078D000C4013F156, REGISTER CONTAINS 2C10B46AEF4B95000000000000000000 VFNTH : LAST EXECUTED FORTRAN STATEMENT IN PROGRAM MAIN (OFFSET 028146). 0TRACEBACK OF CALLING ROUTINES; MODULE ENTRY ADDRESS = 117010. ----------------------------------------------------------------------------- MAIN (117010) CALLED BY OPERATING SYSTEM. ----------------------------------------------------------------------------- 0STANDARD CORRECTIVE ACTION TAKEN. EXECUTION CONTINUING. 0AFB207I VFNTH : PROGRAM INTERRUPT - FLOATING-POINT OVERFLOW EXCEPTION VFNTH : PSW 078D000C4013F14C, REGISTER CONTAINS 3EFFFFFE000001000000000000000000 VFNTH : LAST EXECUTED FORTRAN STATEMENT IN PROGRAM MAIN (OFFSET 02813C). 0TRACEBACK OF CALLING ROUTINES; MODULE ENTRY ADDRESS = 117010. ----------------------------------------------------------------------------- MAIN (117010) CALLED BY OPERATING SYSTEM. ----------------------------------------------------------------------------- 0STANDARD CORRECTIVE ACTION TAKEN. EXECUTION CONTINUING. 0AFB207I VFNTH : PROGRAM INTERRUPT - FLOATING-POINT OVERFLOW EXCEPTION VFNTH : PSW 078D000C4013F156, REGISTER CONTAINS 3EFFFFFE000001000000000000000000 VFNTH : LAST EXECUTED FORTRAN STATEMENT IN PROGRAM MAIN (OFFSET 028146). 0TRACEBACK OF CALLING ROUTINES; MODULE ENTRY ADDRESS = 117010. ----------------------------------------------------------------------------- MAIN (117010) CALLED BY OPERATING SYSTEM. ----------------------------------------------------------------------------- 0STANDARD CORRECTIVE ACTION TAKEN. EXECUTION CONTINUING. 0AFB207I VFNTH : PROGRAM INTERRUPT - FLOATING-POINT OVERFLOW EXCEPTION VFNTH : PSW 078D000C4013F14C, REGISTER CONTAINS 3EFFFFFE000001000000000000000000 VFNTH : LAST EXECUTED FORTRAN STATEMENT IN PROGRAM MAIN (OFFSET 02813C). 0TRACEBACK OF CALLING ROUTINES; MODULE ENTRY ADDRESS = 117010. ----------------------------------------------------------------------------- MAIN (117010) CALLED BY OPERATING SYSTEM. ----------------------------------------------------------------------------- 0STANDARD CORRECTIVE ACTION TAKEN. EXECUTION CONTINUING. DET = 0.1949043E-05 IDET = *** ACTUAL DET. = DET*10**IDET RESULTS FOR SAMPLE 2 0AFB210I VFNTH : PROGRAM INTERRUPT - PROTECTION EXCEPTION VFNTH : PSW 078D000400150032 VFNTH : LAST EXECUTED FORTRAN STATEMENT IN PROGRAM MAIN AT ISN 37 (OFFSET 024D4C). 0TRACEBACK OF CALLING ROUTINES; MODULE ENTRY ADDRESS = 117010. ----------------------------------------------------------------------------- VFWSF# (14FC5E) CALLED BY MAIN (117010) AT ISN 37 AT OFFSET 024D4C. ----------------------------------------------------------------------------- MAIN (117010) CALLED BY OPERATING SYSTEM. ----------------------------------------------------------------------------- 0STANDARD CORRECTIVE ACTION TAKEN. EXECUTION CONTINUING. 0AFB240I VABEX : ABEND CODE IS: SYSTEM 0C4, USER 0. VABEX : PSW=078D000400150032 ENTRY POINT=117010. VABEX : REGS 0 - 3 00000032 00117B0C 001179B2 C1497905 VABEX : REGS 4 - 7 587A0E55 0013BD02 00000008 00118D08 VABEX : REGS 8 -11 00000000 8014FC76 00114FE0 001442B0 VABEX : REGS 12-15 8014FBB0 00144F40 6013BD5C 8014FC5E VABEX : FRGS 0 & 2 42320000 00000000 3E96226B 15979FD5 VABEX : FRGS 4 & 6 B7104D7D 10C00000 3E6F80D6 950CEC80 VABEX : ABEND IN MODULE VFWSF# (OFFSET 000003D4).