ŧY/*JOBPARM REGION=1024,LINES=19}// EXEC FORTGCLG //FORT.SYSIN DD *ȫCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCChC0000~CC0000~CC0xPROGRAMIMSGPCMA0JC C0VERSION2.219-SEP-820~CC0000~CC0000~C8CpPROGRAMIMSGPCMA(IMageSeGmentation,P-channel data,CoMmon\C!covariancematrix,Automatic mode)ISFORSEGMENTINGMULTI-CHANNELHCDIGITALIMAGES.(FORSINGLE-CHANNELDATATHE"IMSG1"ACCPROGRAMSMAYBEUSED.)00JCqCmTHISPROGRAMUSESDISTANCEINTHEMETRICOFTHEPOOLED,/CC)WITHIN-GROUPSCOVARIANCEMATRIX.00CC0000~C&CrAUTOMATIC MODE: TRIESDIFFERENTNUMBERSOFCLASSES,jCCWITHAUTOMATICSETTINGOFINITIALMEANS.USEPROGRAMAC+CIMSGPCMTOTRYA FIXEDNUMBEROFCLASSES.0jCC0000~C0C0000~CQCPROGRAMMEDBY:000C){CDR.STANLEYL.SCLOVE0312/996-26810C/yC*QUANTITATIVEMETHODSDEPARTMENT%312/996-26760C݌CCOLLEGEOFBUSINESSADMINISTRATION00CC!UNIVERSITYOFILLINOISATCHICAGO00C_-CBOX 4348,CHICAGO,IL 6068000CC0000~CC0000~C` CqRESEARCHSUPPORTEDINPARTBY:00C'C0000~CYCONRCONTRACTN00014-80-C-0408,TASKNR042-4430Cd} CAROCONTRACTDAAG29-82-K-015500/ChC0000~CC0000~ChCxPROGRAMRESTRICTIONS(CANBEMODIFIED):0/CmCNR,NUMBEROFROWSOFARRAY,ATMOST256;0jCqCNC,NUMBEROFCOLUMNSOFARRAY,ATMOST256;0/C CIP,NUMBEROFCHANNELS,ATMOST20;01CCKU,MAXIMUMNUMBEROFCLASSES,ATMOST29;0aCuCITER,MAXIMUMNUMBEROFITERATIONS,20.0Cx8C0000~C\C[SUBROUTINE(S)CALLED:00jCHCMATEQ, WHICH CALLS MATDT00CC0000~CqC0000~CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCĝC0000~CC&CGDIMENSIONX(56,56,10)CCHANGENO.OFCHANNELSFROM4TO20CANDNO.OFROWSANDCOLUMNSFROM56TO256K0CINDIMENSIONSTATEMENTSOQCWHENALLOWABLEREGIONISINCREASED.{DIMENSION IOTA(56)yDIMENSIONJOTA(29)SDIMENSIOND(29),ICLUS(56,56)WDIMENSIONTITLE(18)-DIMENSIONNG(29),XMEAN(29,4)DIMENSIONFMT(18)2DIMENSIONSUM(29,4),SS(29,4,4)6DIMENSIONSSD(29,4,4)DIMENSION WGSS(4,4)Y DIMENSIONVARHAT(20,20),WGMS(4,4)}DIMENSIONICLSOL(56,56) DIMENSION XMIN(4),XMAX(4)DIMENSIONIV(20,20) hDIMENSIONP(20,20)DIMENSIONNT(29,29,29),IRSUM(29,29),TP(29,29,29) DIMENSIONPROB(29) CDOUBLEPRECISIONSS,SUMDOUBLEPRECISIONWGSS,SSD8DOUBLEPRECISIONVARHAT\DOUBLEPRECISIONPH DOUBLEPRECISIONDETDOUBLEPRECISIONDqDOUBLEPRECISION XMEANDOUBLEPRECISIONTEMPIV,TEMPJVDOUBLEPRECISIONF &DOUBLEPRECISIONCFDOUBLEPRECISIONTP,PROBC<C|0CIVISAWORK ARRAYFORSUBROUTINEMATEQ.QC{CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC$yC0000~C^CCONTROL"CARDS":000CC0000~C- CDATASET TITLE000C,CNUMBEROF ROWS,NR,INFORMAT(3X,I3)0C?CNUMBEROFCOLUMNS,NC,INFORMAT(3X,I3)0CCNUMBEROFCHANNELS,IP,INFORMAT(3X,I2)0jCCFMT,INFORMAT(18A4), E.G.,(1X,F4.1)0AC:YCe"FMT"WILLALSOBEUSEDFOROUTPUT,SO ALLOWAT LEASTONE BLANC=}CbATTHEBEGINNINGFORCARRIAGECONTROL.0jCbCDATA,ONENUMBERATA TIME,INFORMATSPECIFIEDBYFMTCfCdDATAISINDEXEDBYCHANNEL,ROW,ANDCOLUMN.0ChCfCOLUMNCHANGESFIRST,THENROW,THENCHANNEL.0CC0000~Ck CSMALLESTNUMBEROFCLASSESTOBETRIED,KL,INFORMAT 2X,I2/Co CLARGESTNUMBEROFCLASSESTOBETRIED,KU,INFORMAT 2X,I2CC(2<=KL<=KU<=29)00JCC0000~C~8C0000~Cv\CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCzHC[Cq READ(5,38000) TITLECCWRITEPROGRAMINFORMATION.]&WRITE(6,22000)AWRITE(6,44000) TITLEECICREADNUMBERSOFROWSANDCOLUMNS,NRANDNC.0READ(5,12000)NR}QREAD(5,12000)NCM{WRITE(6,40000)NRQyWRITE(6,42000)NC CREADNUMBEROFCHANNELS,IP.\CbINSTATISTICALTERMINOLOGY,A"CHANNEL"ISA"VARIABLE."U-READ(5,64000)IPYWRITE(6,66000)IPC0CREADDATAFORMAT.4READ(5,38000)FMT7YC} CREAD DATA.DO300 ICHAN=1,IP0DO300I=1,NRh0DO300J=1,NC 0READ(5,FMT)X(I,J,ICHAN)0IF(I.EQ.1 .AND.J.EQ.1)GOTO100 0GOTO2001000CONTINUE0DO200IVAR=1,IP80XMAX(IVAR)=X(1,1,IVAR)\0XMIN(IVAR)=X(1,1,IVAR)H2000CONTINUE0DO300IVAR=1,IPq0IF(X(I,J,IVAR).LT.XMIN(IVAR))XMIN(IVAR)=X0X(I,J,IVAR)0IF(X(I,J,IVAR).GT.XMAX(IVAR))XMAX(IVAR)=&X0X(I,J,IVAR)3000CONTINUEWRITE(6,74000)WRITE(6,FMT)(XMIN(IVAR),IVAR=1,IP)0WRITE(6,76000)(QWRITE(6,FMT)(XMAX(IVAR),IVAR=1,IP){CREADKLANDKU,THEMINANDMAXNUMBEROFCLASSESTOBETRIED.yREAD(5,10000)KLREAD(5,10000)KU)WRITE(6,24000) KL,KU/-CݫCDO5800K= KL,KU_CCOMPUTEINITIALMEANS.0DO400IG=1,KY0DO400IVAR=1,IP}0>XMEAN(IG,IVAR)=XMIN(IVAR)+IG*(XMAX(IVAR)-@X0^XMIN(IVAR))/(K+1)4000CONTINUEdhChWRITE(6,26000)K0DO500IG=1,K 0 WRITE(6,28000)IG,(XMEAN(IG,IVAR),IVAR=1,IP)m 500CONTINUEqC8C\CtHCSETCONSTANTS.xCqN= NR*NCXN=N0DO600INTEG=1,NC& 0IOTA(INTEG)= INTEG 600CONTINUE0DO700J=1,KC 0JOTA(J)=JF0 700CONTINUEQC{PI=3.1415927KyC;DISTRIBUTIONALPARAMETERS:OCKMEANVECTORSOFDIMENSIONIPANDANIP-BY-IPCOVARIANCEMATRIX, CWHEREIPISTHENUMBEROFCHANNELS-CNO.OFINDEPENDENTTRANSITIONPROB.'S=(K**2)*(K-1)SCWCNO.OFPARAMETERS=NO.OFDISTRIBUTIONALPARAMETERSC0A+NO.OFINDEPENDENTTRANSITIONPROB.'SNOPARM=K*IP+IP*(IP+1)/2+(K**2)*(K-1)1YC!INITIALIZEVARIABLES.5}C9CFOR FIRSTITERATION,MARGINALDISTRIBUTIONOFLABELS, PROB,CIS TAKENTOBEUNIFORM.h0DO800IG=1,K 800PROB(IG)= 1.0/KC CITER=1  900CONTINUE8 WRITE(6,34000)ITER\IF (ITER.EQ.1)GOTO1100H0DO1000I=1,NR0DO1000J=1,NCq10000ICLSOL(I,J)=ICLUS(I,J)CSAVEPREVIOUS VALUEOF-2LOGMAXLIKELIHOOD.XMN2OL=XMN2LL&1100CONTINUECCCOMMENCEDISTANCECOMPUTATIONS. C00DO2300I=1,NRQ0DO2300J=1,NC.{C!INITIALIZEDISTANCES(TOBEACCUMULATED)AT ZERO.|y0DO1200L=1,K12000D(L)=0.0CFOR FIRSTITERATION,EUCLIDEANDISTANCEISUSEDBECAUSE$-CNOCOVARIANCEMATRIXISYETAVAILABLE. AFTERTHE FIRST^C!ITERATION,DISTANCEWILLBE TAKENINTHEMETRICOFTHE C!COVARIANCEMATRIX.0IF (ITER.GT.1)GOTO1400,0DO1300 L=1,K>Y00DO1300IVAR=1,IP}13000D(L)=D(L)+(XMEAN(L,IVAR)-X(I,J,X05IVAR))**2: 0GOTO1600=h14000CONTINUEbC(AFTER FIRSTITERATION)COMPUTEDISTANCESINMETRICfCOFCOVARIANCEMATRIX 0DO1500L=1,K00DO1500IVAR=1,IPk00TEMPIV=XMEAN(L,IVAR)-X(I,J,IVAR)n800DO1500JV=1,IPr\00TEMPJV=XMEAN(L,JV)-X(I,J,JV)H150000D(L)=D(L)+TEMPIV*P(IVAR,JV)*~ X00TEMPJVvq16000CONTINUEz0DO2100L=1,K[CFOR FIRSTITERATION,CLASSIFICATIONISSIMPLYBY&CMINIMUMDISTANCE.OTHERWISE,THETRANSITIONC)PROBABILITIESENTER.0IF (ITER.EQ.1)GOTO2100]0ARG=-D(L)/2.0{00IF (ARG.LT.-180.2)GOTO1700DQ 0GOTO1800H{17000D(L)=0.0y 0GOTO1900}18000CONTINUEMCQ-0D(L)=EXP(ARG)19000CONTINUE\CU0IF(I.EQ.1.OR.J.EQ.1)GOTOYX0u2000Y0IG1=ICLUS(I,J-1)}0IG2=ICLUS(I-1,J)30D(L)=-TP(IG1,IG2,L)*D(L)7 0GOTO2100h20000CONTINUECCLASSIFYBOUNDARYOBSERVATIONS0D(L)=-D(L)*PROB(L) 21000CONTINUE C0/F=D(1)80ICLUS(I,J)=1\0DO2300L=2,KH0IF(D(L)-F)2200,2300,230022000/F=D(L)q0ICLUS(I,J)=L23000CONTINUEC &WRITE(6,14000)WRITE(6,16000)(IOTA(J),J=1,NC)0DO2400I=1,NR0 WRITE(6,18000)I,(ICLUS(I,J),J=1,NC)02400CONTINUEQC{0DO2500IG=1,Ky 0NG(IG)=0( 0DO2500IVAR=1,IP0%SUM(IG,IVAR)=0.0-0DO2500JV=1,IP0SSD(IG,IVAR,JV)=0.0)25000>SS(IG,IVAR,JV)=0.0/0DO2700I=1,NR0DO2700J=1,NCY 0IGROUP=ICLUS(I,J)%}0NG(IGROUP)=NG(IGROUP)+1׼ 0DO2600IVAR=1,IP0SUM(IGROUP,IVAR)=SUM(IGROUP,IVAR)+X(I,J,hX0vIVAR)@00DO2600JV=1,IP26000SS(IGROUP,IVAR,JV)=SS(IGROUP,IVAR,JV)+d X00X(I,J,IVAR)*X(I,J,JV)h27000CONTINUE 0DO2800IVAR=1,IP80DO2800JV=1,IPl\ 0_WGSS(IVAR,JV)=0.0pH28000CONTINUECqCt0DO3100IG=1,Kx0IF(NG(IG).EQ.0)GOTO2900& 0GOTO30002900WRITE(6,88000)IG 0GOTO53003000CONTINUE0 0DO3100IVAR=1,IPQCCOMPUTEMEANVECTORS.B{0>XMEAN(IG,IVAR)=SUM(IG,IVAR)/NG(IG)FyCCOMPUTESUM-OF-PRODUCTSMATRICES.ˌ0DO3100JV=1,IP0CF=SUM(IG,IVAR)*SUM(IG,JV)/NG(IG)K-0SSD(IG,IVAR,JV)=SS(IG,IVAR,JV)-CFO31000CONTINUECCSCVYCZ}CPOOL:C10DO3200IG=1,K5h 0DO3200IVAR=1,IP90DO3200JV=1,IP0]WGSS(IVAR,JV)=WGSS(IVAR,JV)+SSD(IG,IVAR,  X0tJV) 32000CONTINUE C 8C \CCOMPUTEVARHAT,MLEOFCOMMONCOVARIANCEMATRIX: HC  0DO3300IVAR=1,IP q0DO3300JV=1,IP 0?VARHAT(IVAR,JV)=WGSS(IVAR,JV)/N 33000CONTINUE &CPRINTVARHAT: WRITE(6,54000) 0DO3400IVAR=1,IP 0 WRITE(6,62000)(VARHAT(IVAR,JV),JV=1,IP) 03400CONTINUE QC {IDET=1 yNRS1=0 NVRHT=20 .NP=20 |-CALLMATEQ(VARHAT,IP,NVRHT,JFLG,DET,IDET,IV,NRS1,P,NP) C CGENERALFORMOFCALLIS: $CCALLMATEQ(A,M,N,JFLG,DET,IDET,IV,NRS1,P,LL) ^CSEELISTINGOFSUBROUTINEFORFULLEREXPLANATION. YC$DET(VARHAT)=DET*10.0**IDET }C IF (JFLG.GT.0) WRITE(6,70000)JFLG >WRITE(6,72000)DET,IDET hC WRITE(6,58000) :0DO3500IVAR=1,IP = 0 WRITE(6,62000)(P(IVAR,JV),JV= 1,IP) b3500CONTINUE fC i8 XIDET=IDET \XLGDET=DLOG(DET)+XIDET*ALOG(10.0) jHXMN2LL=N*(IP*ALOG(2.0*PI)+IP+XLGDET) nC rqWRITE(6,32000)XMN2LL C ~0DO3600IG=1,K v&0 WRITE(6,30000)IG,(XMEAN(IG,IVAR) ,IVAR=1,IP) z3600CONTINUE [IF (ITER.EQ.1)GOTO3800 0DO3700I=1,NR 00DO3700J=1,NC Q0IF(ICLUS(I,J).EQ.ICLSOL(I,J))GOTO3700 { 0GOTO3800 {y37000CONTINUE D GOTO5400 H3800CONTINUE -C }CCOMPUTETRANSITIONPROBABILITYMATRIX: MC Q0DO3900I1=1,K 0DO3900I2=1,K Y0DO3900J=1,K T}0NT(I1,I2,J)=0 X39000CONTINUE 0DO4000I=2,NR h0DO4000J=2,NC 3 0IM1=I-1 7 0JM1=J-1 0INORTH=ICLUS(IM1,J) 0IWEST=ICLUS(I,JM1)0IY=ICLUS(I,J)80?NT(IWEST,INORTH,IY)=NT(IWEST,INORTH,IY)+1\40000CONTINUEH0DO4100I1=1,K0DO4100I2=1,Kq0%IRSUM(I1,I2)=041000CONTINUE0DO4200I1=1,K&0DO4200I2=1,K0DO4200J=1,K0IRSUM(I1,I2)=IRSUM(I1,I2)+NT(I1,I2,J)42000CONTINUE00DO4500I1=1,KQ0DO4500I2=1,K{0?XDENOM=IRSUM(I1,I2)y0IF(XDENOM.EQ.0.0)GOTO4300 0GOTO440043000XDENOM=K-44000CONTINUE(C0DO4500J=1,K0XNUM=NT(I1,I2,J)CIF THEREARENOTRANSITIONSFROM(I1,I2),THENTP(I1,I2,J)*YCISSET EQUALTO ZERO,FORALLJ=1,2,...,K.-}C45000,TP(I1,I2,J)=XNUM/XDENOM4600CONTINUE%hCCCOMPUTEMARGINALDISTRIBUTIONOFLABELS:0DO4700IG=1,K 0PROB(IG)=NG(IG)@0PROB(IG)=PROB(IG)/N4700CONTINUEc8Cg\ TRANS=0.0H0DO4800I1=1,Kŵ0DO4800I2=1,Klq0DO4800 J=1,Kp 0ITEST=NT(I1,I2,J)0IF(ITEST.EQ.0)GOTO4800& 0IF(TP(I1,I2,J).EQ.0.)GOTO4800t0TRANS= TRANS+NT(I1,I2,J)*DLOG(TP(I1,I2,J))x48000CONTINUEC N0TRANS=-2.0*TRANSQC{CACCOUNTFORLABELSOFBORDEROBSERVATIONS:y FIRST=0.0،DO516J=1,NCBLABEL1=ICLUS(1,J)F-PROBAB=PROB(LABEL1)˫FIRST= FIRST+ALOG(PROBAB) 516CONTINUEKDO517I=2,NROLABEL1=ICLUS(I,1)RYPROBAB=PROB(LABEL1)}FIRST= FIRST+ALOG(PROBAB) 517CONTINUEVFIRST=-2.0*FIRSTZhCCCOMPUTE MODELSELECTIONCRITERIA:1CEXCEPTFOR FIRSTITERATION,COMPUTEAIC BASEDON5 CNEWCLUSTERINGANDOLDCOVARIANCEMATRIX.9IF (ITER.EQ.1)GOTO4900AICOLD=XMN2OL+2.0*NOPARM8SCHOLD=XMN2OL+ALOG(XN)*NOPARM\WRITE(6,92000)AICOLDHWRITE(6,92005)SCHOLD 4900CONTINUEqCwXMN2LL=XMN2LL+ TRANS+ FIRSTC&AIC=XMN2LL+2.0*NOPARMWRITE(6,94000)AICSCH=XMN2LL+ALOG(XN)*NOPARMWRITE(6,94005)SCH0CQC{CWRITETRANSITIONPROBABILITIES:yWRITE(6,46000) WRITE(6,50000)(JOTA(JAY),JAY=1,K)0DO5000I1=1,K-0DO5000I2=1,K0WRITE(6,48000)I1,I2,(NT(I1,I2,J),J=1,K).50000CONTINUE|WRITE(6,52000) WRITE(6,50000)(JOTA(JAY),JAY=1,K)Y0DO5100I1=1,K!}0DO5100I2=1,K;0WRITE(6,56000)I1,I2,(TP(I1,I2,J),J=1,K)51000CONTINUEhWRITE(6,90000)(PROB(IG),IG=1,K)WRITE(6,80000) TRANS>ITER=ITER+1 IF(ITER.GE.21)GOTO5200GOTO900: 5200WRITE(6,78000)'85300STOPa\CeHCi5400CONTINUEqCj WRITE(6,68000)ITERnWRITE(6,36000)(NG(IG),IG=1,K)r&WRITE(6,60000) 0DO5500IVAR=1,IP~0DO5500JV=1,IPv0_WGMS(IVAR,JV)=WGSS(IVAR,JV)/(N-K)y055000CONTINUEQ0DO5600IVAR=1,IP{0 WRITE(6,62000)(WGMS(IVAR,JV),JV=1,IP)y5600CONTINUECC{-WRITE(6,16000)(IOTA(J),J=1,NC)D0DO5700I=1,NRH0 WRITE(6,18000)I,(ICLUS(I,J),J=1,NC)5700CONTINUE}CLYWRITE(6,84000)NOPARMP}AIC=XMN2LL+2.0*NOPARMWRITE(6,82000)AICCThCX WRITE(6,86000)ITERC 5800CONTINUE3C7 STOP810000 FORMAT(2X,I2)\12000 FORMAT(3X,I3)H14000 FORMAT(//1X,'SEGMENTATION:'/) 16000FORMAT(/,1X,'ROW:gCOLUMN:'/,4X,(40I3/))q18000 FORMAT(1X, I3, (40I3/) )20000 FORMAT(1X, (8E13.5/)//)22000 FORMAT('1','#######################################*',&;X//,1X,'PROGAMIMSGPCMA'/X,1X,'FOR IMAGESEGMENTATION'/$X,1X,'USINGDISTANCEINTHEMETRICOFTHECOVARIANCEMATRIX'/^X//,1X,'DEVELOPEDANDPROGRAMMEDBYDR.STANLEYL.SCLOVE' 0^X//,1X,'VERSION2.219-SEP-82'//)Q24000 FORMAT('1',//,1X,'MIN AND MAX NUMBER OF CLASSES TO BE TRIED ARE ',{XI2,'AND',I2//)y26000 FORMAT(1H1,'K = ',I2,' CLASSES'//)28000 FORMAT(/1X,'INITIAL MEAN VECTOR FOR CLASS ',I2,': ',(4E14.5/)) 30000FORMAT(1X,'MEANVECTORFOR CLASS',I2,':',(8E13.5/))- 32000FORMAT(/1X,'MINUS2LOGLIKELIHOOD=',E13.5//)34000 FORMAT(/,1X,'ITERATION ', I2,//)36000 FORMAT(/,1X,'NUMBERS:',(9I13/)/)38000 FORMAT(18A4)(40000 FORMAT(/1X,'NUMBER OF ROWS = ',I3/)&Y42000 FORMAT(/1X,'NUMBER OF COLUMNS = ',I3/)} 44000 FORMAT(1X,18A4)46000 FORMAT(/1X,'TRANSITIONS'/)*48000 FORMAT(1X, 2I4, (9I7/))-h50000 FORMAT(/9X,(9I7/))52000 FORMAT(/1X,'TRANSITION PROBABILITIES'/)54000 FORMAT(///,1X,'COMMON COVARIANCE MATRIX (MLE):',//)% 56000 FORMAT(1X,2I4,3X,(9F7.4/))58000 FORMAT(///,1X,'INVERSE COVARIANCE MATRIX:',//)60000 FORMAT(///,1X,'COMMON COVARIANCE MATRIX (UNBIASED ESTIMATE):',//)862000 FORMAT(1X,(8E13.5/))#\64000 FORMAT(3X,I2)"H66000 FORMAT(/1X,'NUMBER OF CHANNELS = ',I2/)c 68000 FORMAT(/1X,'CONVERGENCE: NO CASE CHANGED CLASSES AFTER ',gq$X'ITERATION',I2,'.RESULTSAREPRINTEDBELOW.'//)70000 FORMAT(/,1X,'JFLG = ',I2,'. IF JFLG=0, COMPUTATION OF DET',ŝX'WENT WELL;OTHERWISE, THEREWASTROUBLEORMATRIXWAS',l&^X'ILL-CONDITIONED.'//)p 72000FORMAT(/1X,'DET=',E13.5,'IDET=',I3,5X,X'ACTUALDET.=DET*10**IDET',//) 074000 FORMAT(/1X,'MINIMUM FOR EACH CHANNEL: ',/)s076000 FORMAT(/1X,'MAXIMUM FOR EACH CHANNEL: ',/)wQ78000 FORMAT(1X,'ROUTINE HAS NOT CONVERGED IN 20 ITERATIONS. STOP'){80000 FORMAT(//,1X,'CONTRIBUTION OF TRANS. PROBS. TO ',yX'LOGLIKELIHOOD=',E15.5/)82000 FORMAT(1X,'AIC = ', E15.5///)84000 FORMAT(//1X,'NUMBER OF PARAMETERS = ',I4//)-86000 FORMAT(1X,'RESULTS DID NOT CHANGE IN ITERATION NUMBER',I3,'.',ثX/,1X,'PROGRAM ENDEDSUCCESSFULLY.')B 88000FORMAT(1X,'NOOBSERVATIONSIN GROUP',I3,'.STOP')F90000 FORMAT(/1X,'MARGINAL PROB. VECTOR:',(9F11.4/))92000 FORMAT(/1X,'AIC BASED ON NEW LABELS AND OLD DISTRIBUTIONAL ',Y)X'PARAMETERS:',E15.5/)J}92005 FORMAT(/1X,'SCHWARZ CRITERION BASED ON NEW LABELS AND OLD ',N^X'DISTRIBUTIONALPARAMETERS:',E15.5/)RȀ94000 FORMAT(/1X,'AIC BASED ON NEW LABELS AND NEW DISTRIBUTIONAL ',h)X'PARAMETERS:',E15.5/)׀94005 FORMAT(/1X,'SCHWARZ CRITERION BASED ON NEW LABELS AND NEW ',DV^X'DISTRIBUTIONALPARAMETERS:',E15.5/)Z ENDSUBROUTINEMATEQ(A,M,N,JFLG,DET,IDET,IV,NRS1,P,LL)1C48C!SUBROUTINE MATEQISDMATEQFROMTHEUICSUBROUTINELIBRARY.8\CH C0SUBROUTINEDMATEQC0*****************qCTHISROUTINEWILL SOLVEAREAL*8SYSTEMOFLINEAREQUATIONS,COMPUTECTHEDETERMINANT,WITHOUTUNDERFLOWOROVERFLOW,OFAREAL*8MATRIX, CAND/ORINVERTAREAL*8MATRIX.&CCALLINGSEQUENCE:CCALLDMATEQ(A,N,IA,JFLG,DET,IDET,IV,NRS,P,IP)WHERE;CjA(INPUT)-ISTHEREAL*8MATRIXON WHICHTHEROUTINEISC0TO WORK.INTHEPROCESSOFCOMPUTATIONTHE0C0CONTENTSOFTHISMATRIXAREDESTROYED.QCjN(INPUT)-ISANINTEGER*4VARIABLE WHICHSPECIFIESTHE{C0ORDEROFTHEAMATRIX.yCkIA(INPUT)-ISANINTEGER*4VARIABLE WHICHSPECIFIESTHEC0ACTUALROWDIMENSIONOFAASDIMENSIONEDINC0THECALLINGPROGRAM.IAMUSTBEGREATERTHAN-C0OR EQUALTON.CmJFLG(OUTPUT)-ISANINTEGER*4RETURNCODEVARIABLE.UPONC0RETURNFROMDMATEQIF;C0JFLG=0,ALLWENT WELL.C0JFLG=1,THEAMATRIXWASSINGULARORNEARYC0SINGULARANDTHECOMPUTATIONS COULDNOTBE+}C0COMPLETED.THECONTENTSOFTHEVARIABLESC0A,DET,IDETANDPAREMEANINGLESS.ClDET(OUTPUT)-ISAREAL*8VARIABLE WHICHCONTAINSTHE!hC0DETERMINANTOFA.(SEE IDET);CmIDET(INPUT)-ISANINTEGER*4VARIABLE.ON INPUTIF;ߨC0IDET=0,NODETERMINANTISCALCULATED. C0IDETNOT0,THEDETERMINANTOFAISCOMPUTED.C0KONOUTPUTIDETCONTAINSTHE POWEROF10>C0THATDETSHOULDBEMULTIPLIEDBYTOGIVETHE8C0CORRECT VALUEOFTHEDETERMINANT.I.E.\C0DET(A)=DET*10.0D0**IDET.`HC0KIFDET(A)CANBECOMPUTEDWITHOUT UNDEROR'C0OVERFLOW,THENIDET=0OTHERWISEIDETISSETaqC0TOTHEPROPER VALUESOTHATNO UNDEROR OVER-e C0FLOWWILL OCCURINCOMPUTINGDET.iCkIV(INPUT)-ISANINTEGER*4WORK ARRAY WHICHSHOULDBE&C0DIMENSIONEDAT LEASTIV(N).jCnClNRS(INPUT)-ISANINTEGER*4VARIABLEWITHTHEFOLLOWINGrC0INTERPRETATION:0C0NRS>0, SOLVEASYSTEMOFLINEAREQUATIONSQC0WITHNRS RIGHTHANDSIDES.u{C0NRS=0,INVERTTHEAMATRIX.yyC0NRS<0,ONLYCOMPUTETHEDETERMINANTOFA.ڌC0INTHISCASEIDETMUSTBEDIFFERENTFROM0.CjP(INPUT)-ISAREAL*8 ARRAYWITHTHEFOLLOWINGINTER--C0PRETATION:C0IFNRS>0,THENPCONTAINSTHENRS RIGHTHANDC0SIDESSTOREDBYCOLUMNS.INTHISCASEPMUST{C0BEDIMENSIONEDAT LEASTP(N,NRS).ONRETURNDC0THECOLUMNSOFPAREREPLACEDBYTHERESPEC-GYC0TIVESOLUTIONS.}C0IFNRS=0,THENPMUSTBEDIMENSIONEDAT LEASTC0P(N,N).ONRETURNPWILLCONTAINTHEINVERSELC0OFA.PhC0IFNRS<0,THENPNEEDONLYBEA DUMMYVARIABLEC0INTHISCASEPIS NEVERACCESSEDBYDMATEQ.CkIP(INPUT)-ISANINTEGER*4VARIABLE WHICHCONTAINSTHET C0ACTUALROWDIMENSIONOFPASDIMENSIONEDINXC0THECALLINGPROGRAM.IPMUSTBEGREATERTHANC0OR EQUALTON. 68CNOTE:IMMEDIATELYONRETURNFROMDMATEQTHECONDITIONCODE FLAG,2\CvJFLG,SHOULDBEINTERROGATED.IFJFLG=1,THENTHEROUTINE6HCvCOULDNOTCOMPUTEASOLUTION.CMETHOD-THEALGORITHMUSEDISGAUSSIANELIMINATIONWITHPARTIALqC0000-1CrPIVOTING.INESSENCETHEROUTINEGENERATESAMATRIXLSUCH C2-1& CmTHATL*A=U, WHEREUISAN UPPERTRIANGULARMATRIX.THENITCoSOLVESTHESYSTEMA*X=PBY MEANSOFTHEEQUIVALENTSYSTEMC0-1b-1ClU*X=L*A*X=L*PBYBACKSUBSTITUTION.0 C2-1QCtTHELMATRIXCANBEWRITTENASAPRODUCTOFTHEFORM{ Ck-1 yCaL=L*P*....*L*P WHEREEACHPISAPERMUTATIONC3N-1N-1a11JKCfMATRIXOBTAINEDBYINTERCHANGINGATMOSTTWOROWSOFTHE-ChIDENTITYMATRIX.(THISREPRESENTSTHEINTERCHANGINGOFTWOCfROWS).THELMATRICESAREELIMINATIONMATRICES WHICHARE C0K CfCHOSENTOINTRODUCE ZEROSINTHELASTN-KENTRIESOFTHEK-THCfCOLUMNOFTHEMATRIX.YC}C02-1b-1<ClTHECALCULATIONSOFL*AANDL*PAREDONEBYPERFORMING&CTHEPERMUTATIONSONAANDPRESPECTIVELY.THEACTUALLANDPhC0000K/KCARENOTCOMPUTED.*CSUBROUTINESCALLED:DMATDT- CREFERENCE:C G. W. STEWART, INTRODUCTION TO MATRIX COMPUTATIONS,C ACADEMIC PRESS, 1973.,8REAL*8A(N,1),DET,P(LL,1)?\REAL*8DNORM,DEN,DMULT ,DSUM,DISIGNHDIMENSION IV(1)ֵNRS=NRS1#q IF(NRS.EQ.0)IDET=1">DISIGN=1.0D+00c ,DET=0.0D+00g&JFLG=0CC JFLG IS A TROUBLE FLAG.UPON EXIT IF JFLG=0 THEN THE MATRIX WAS PROCESSlC WITHOUT TROUBLE.IF JFLG=1 EITHER THE MATRIX IS SINGULAR OR TROUBLEo0C OCCURED.ISIGN=-ISIGN EVERY TIME A ROW IS INTERCHANGED.THIS IS USED TOQC INSURE THAT THE DETERMINANT HAS THE PROPER SIGN.{CsyM1=M-1w DO100 I=1,M 100IV(I)=I- IF (NRS)500,200,500200DO300 I=1,M0DO300 J=1,M300P(I,J)=0.0D+00 DO400 I=1,MAY400>P(I,I)=1.0D+00E}NRS=MICȫC INSTEAD OF ACTUALLY INTERCHANGING ROWS A POINTER ARRAY IS USED TO KEEPJhC TRACK OF THE ROW POSITIONS.NCRC BEGIN ELIMINATION LOOP. C500DO1200K=1,M1VICOL=KY8IPCOL=K\C0HC SEARCHING FOR LARGEST ELEMENT IN ABSOLUTE VALUE IN COLUMN K.4C8qDNORM=A(IV(K),K)IFLG=0 KK=K+1 &0DO600J=KK,M 0IF(DABS(A(IV(J),K)).LE.DABS(DNORM))GOTO600 0IFLG=1  0IPCOL=IV(J) 00DNORM=A(IPCOL,K) Q 600CONTINUE {C yC IF IFLG=0 NO ROW INTERCHANGE TOOK PLACE.IF IFLG=1 A ROW INTERCHANGE C TOOK PLACE AND THE POINTER ARRAY IV MUST BE UPDATED. C -IF(IFLG.EQ.0)GOTO800 ISAVE=IV(ICOL) IV(ICOL)=IPCOL  ICOL1=ICOL+1 0DO700L=ICOL1,M Y0IF(IV(L).EQ.IPCOL)IV(L)=ISAVE } 700CONTINUE DISIGN=-DISIGN 800IF(DNORM.EQ.0.0D+00)GOTO1900 +hC C BEGIN ELIMINATION OF ROW BELOW IV(K).DEN IS THE PIVOT ELEMENT. C ! K1=K+1 ;0DO1100IM=K1,M C 8C BEFORE ACTUALLY ELIMINTING WE CHECK TO SEE IF A(IV(IM),K) HAS ALREADY \C BEEN ANIHALATED. _HC 0IF(A(IV(IM),K).EQ.0.0D+00)GOTO1100 qC `C CACULATE ELIMINATION FACTOR. 'C a& 0DMULT=-A(IV(IM),K) eC iC WE NOW CALCULATE VALUE OF OTHER ELEMENTS IN ROW IV(IM). C 00DO900NN=K1,M mQ9000A(IV(IM),NN)=(DMULT*A(IV(K),NN))/DNORM+A(IV(IM),NN) q{0IF(NRS.LE.0)GOTO1100 y0DO1000IN=1,NRS Ȍ10000P(IV(IM),IN)=(DMULT*P(IV(K),IN))/DNORM+P(IV(IM),IN) u1100CONTINUE y-1200CONTINUE ګC C CALCULATE VALUE OF DETERMINANT. C IF(A(IV(M),M).EQ.0.0D0)GOTO1900 Y DET=DISIGN }IF(IDET.NE.0)CALLDMATDT(A,N,M,DET,IV,IDET) CIF(DET.EQ.0.0D+00)GOTO1900 GIF(NRS.LE.0)GOTO2000 hC ׫C WE START SOLVING RIGHT HAND SIDES.THE SOLUTION REPLACES THE RIGHT HAND LC VECTOR. P C  1300N1=M-1 DO1600JJ=1,NRS S8C W\C BEGIN BACK SUBSTITUTION. HC P(IV(M),JJ)=P(IV(M),JJ)/A(IV(M),M) 2q0DO1500I=1,N1 60DSUM=0.0D+00 0DO1400 J=1,I &14000DSUM=DSUM-A(IV(M-I),M-J+1)*P(I V(M-J+1),JJ) 4 1500P(IV(M-I),JJ)=(P(IV(M-I),JJ)+DSUM)/A(IV(M-I),M-I) 1600CONTINUE DO1800JJ=1,NRS 00DO1700IND=1,M Q1700A(IND,1)=P(IV(IND),JJ) {0DO1800IND=1,M y1800P(IND,JJ)=A(IND,1) RETURN  1900JFLG=1 -IDET=0  2000RETURN END SUBROUTINEMATDT(A,IA,N,DET,IV,IDET) C!SUBROUTINE MATDTISDMATDTFROMTHEUICSUBROUTINELIBRARY. YREAL*8A(IA,1),DET,B,LOG16 }INTEGER*4IV(1),K ,EQUIVALENCE (B,K) %NUM=16777216 h?LOG16=.120411998265592457D+01 <IF(A(IV(N),N).EQ.0.0D+00)GOTO300 &L=0 DO100 I=1,N B=DABS(A(IV(I),I)) * K=K/NUM-64 ^8L=L+K \100?DET=DET*(A(IV(I),I)/16.0D+00**K) H ,B=DABS(DET) , K=K/NUM-64 ?qIW=L+K IF((IW.LT.-64).OR .(IW.GT.63))GOTO200 ֝ ?DET=DET*16.0D+00**L #&IDET=0 "GOTO400 c200DET=DET*16.0D+00**(-K) gIDET=L+K 0 %B=IDET*LOG16 QIDET=B k{?B=B-DFLOAT(IDET) oy?DET=DET*1.0D+01**B GOTO400 300DET=0.0D+00 s-IDET=0 w 400RETURN END //GO.SYSIN DD * FISHER IRIS DATA YNR=015 }NC=010 ] IP=04 A(1X,F3.1) Eh 5.1VARIABLE1A1001 I 4.9VARIABLE1A1002 4.7VARIABLE1A1003 J 4.6VARIABLE1A1004 N 5.0VARIABLE1A1005 R 5.4VARIABLE1A1006 8 4.6VARIABLE1A1007 \\ 5.0VARIABLE1A1008 UH 4.4VARIABLE1A1009 Y 4.9VARIABLE1A1010 q 5.4VARIABLE1A1011 0 4.8VARIABLE1A1012 4 4.8VARIABLE1A1013 8& 4.3VARIABLE1A1014  5.8VARIABLE1A1015  5.7VARIABLE1A1016  5.4VARIABLE1A1017 0 5.1VARIABLE1A1018 Q 5.7VARIABLE1A1019 { 5.1VARIABLE1A1020 y 5.4VARIABLE1A1021  5.1VARIABLE1A1022  4.6VARIABLE1A1023 - 5.1VARIABLE1A1024  4.8VARIABLE1A1025  5.0VARIABLE1A1026  5.0VARIABLE1A1027  5.2VARIABLE1A1028 Y 5.2VARIABLE1A1029 } 4.7VARIABLE1A1030  4.8VARIABLE1A1031  5.4VARIABLE1A1032 h 5.2VARIABLE1A1033 5.5VARIABLE1A1034 4.9VARIABLE1A1035 + 5.0VARIABLE1A1036  5.5VARIABLE1A1037  4.9VARIABLE1A1038 8 4.4VARIABLE1A1039 )\ 5.1VARIABLE1A1040 /H 5.0VARIABLE1A1041 ݵ 4.5VARIABLE1A1042 q 4.4VARIABLE1A1043 _ 5.0VARIABLE1A1044 5.1VARIABLE1A1045 & 4.8VARIABLE1A1046 ` 5.1VARIABLE1A1047 ' 4.6VARIABLE1A1048 a 5.3VARIABLE1A1049 d0 5.0VARIABLE1A1050 hQ 7.0VARIABLE1A2051 { 6.4VARIABLE1A2052 y 6.9VARIABLE1A2053 m 5.5VARIABLE1A2054 q 6.5VARIABLE1A2055 - 5.7VARIABLE1A2056 ȫ 6.3VARIABLE1A2057 u 4.9VARIABLE1A2058 y 6.6VARIABLE1A2059 5.2VARIABLE1A2060 Y 5.0VARIABLE1A2061 } 5.9VARIABLE1A2062 6.0VARIABLE1A2063 6.1VARIABLE1A2064 h 5.6VARIABLE1A2065 C 6.7VARIABLE1A2066 G 5.6VARIABLE1A2067 5.8VARIABLE1A2068  6.2VARIABLE1A2069 L 5.6VARIABLE1A2070 O8 5.9VARIABLE1A2071 \ 6.1VARIABLE1A2072 H 6.3VARIABLE1A2073 S 6.1VARIABLE1A2074 Wq 6.4VARIABLE1A2075 6.6VARIABLE1A2076 6.8VARIABLE1A2077 2& 6.7VARIABLE1A2078 6 6.0VARIABLE1A2079  5.7VARIABLE1A2080  5.5VARIABLE1A20810 5.5VARIABLE1A2082 Q 5.8VARIABLE1A2083{ 6.0VARIABLE1A2084 y 5.4VARIABLE1A2085 6.0VARIABLE1A2086 6.7VARIABLE1A2087- 6.3VARIABLE1A2088 5.6VARIABLE1A2089 5.5VARIABLE1A2090  5.5VARIABLE1A2091 6.1VARIABLE1A2092Y 5.8VARIABLE1A2093} 5.0VARIABLE1A2094 5.6VARIABLE1A2095 5.7VARIABLE1A2096h 5.7VARIABLE1A2097 6.2VARIABLE1A2098 5.1VARIABLE1A2099 5.7VARIABLE1A2100< 6.3VARIABLE1A3101& 5.8VARIABLE1A31028 7.1VARIABLE1A3103\ 6.3VARIABLE1A3104$H 6.5VARIABLE1A3105^ 7.6VARIABLE1A3106q 4.9VARIABLE1A3107 7.3VARIABLE1A3108, 6.7VARIABLE1A3109?& 7.2VARIABLE1A3110 6.5VARIABLE1A3111 6.4VARIABLE1A3112# 6.8VARIABLE1A3113=0 5.7VARIABLE1A3114 4bQ 5.8VARIABLE1A3115f{ 6.4VARIABLE1A3116y 6.5VARIABLE1A3117 7.7VARIABLE1A3118k 7.7VARIABLE1A3119o- 6.0VARIABLE1A3120 6.9VARIABLE1A3121 5.6VARIABLE1A3122s 7.7VARIABLE1A3123w 6.3VARIABLE1A3124zY 6.7VARIABLE1A3125[} 7.2VARIABLE1A3126 6.2VARIABLE1A3127 6.1VARIABLE1A3128h 6.4VARIABLE1A3129] 7.2VARIABLE1A3130A 7.4VARIABLE1A3131E 7.9VARIABLE1A3132I 6.4VARIABLE1A3133 6.3VARIABLE1A3134}8 6.1VARIABLE1A3135M\ 7.7VARIABLE1A3136QH 6.3VARIABLE1A3137 6.4VARIABLE1A3138\q 6.0VARIABLE1A3139U 6.9VARIABLE1A3140Y 6.7VARIABLE1A3141& 6.9VARIABLE1A31420 5.8VARIABLE1A31434 6.8VARIABLE1A31448 6.7VARIABLE1A31450 6.7VARIABLE1A3146Q 6.3VARIABLE1A3147{ 6.5VARIABLE1A3148y 6.2VARIABLE1A3149 5.9VARIABLE1A3150 3.5VARIABLE21001- 3.0VARIABLE21002 3.2VARIABLE21003 3.1VARIABLE21004 3.6VARIABLE21005 3.9VARIABLE21006Y 3.4VARIABLE21007} 3.4VARIABLE21008 2.9VARIABLE21009 3.1VARIABLE21010h 3.7VARIABLE21011 3.4VARIABLE21012 3.0VARIABLE21013 3.0VARIABLE21014 4.0VARIABLE21015 4.4VARIABLE21016(8 3.9VARIABLE21017\ 3.5VARIABLE21018H 3.8VARIABLE21019 3.8VARIABLE21020)q 3.4VARIABLE21021/ 3.7VARIABLE21022ݝ 3.6VARIABLE21023& 3.3VARIABLE21024_ 3.4VARIABLE21025 3.0VARIABLE21026 3.4VARIABLE210270 3.5VARIABLE21028@Q 3.4VARIABLE21029{ 3.2VARIABLE21030dy 3.1VARIABLE21031h 3.4VARIABLE21032 4.1VARIABLE21033- 4.2VARIABLE21034m 3.1VARIABLE21035q 3.2VARIABLE21036 3.5VARIABLE21037 3.6VARIABLE21038tY 3.0VARIABLE21039x} 3.4VARIABLE21040 3.5VARIABLE21041 2.3VARIABLE21042h 3.2VARIABLE21043 3.5VARIABLE21044٨ 3.8VARIABLE21045 3.0VARIABLE21046C 3.8VARIABLE21047G 3.2VARIABLE210488 3.7VARIABLE21049\ 3.3VARIABLE21050KH 3.2VARIABLE22051O 3.2VARIABLE22052q 3.1VARIABLE22053 2.3VARIABLE22054S 2.8VARIABLE22055W& 2.8VARIABLE22056 3.3VARIABLE22057 2.4VARIABLE220582 2.9VARIABLE2205950 2.7VARIABLE220609Q 2.0VARIABLE22061{ 3.0VARIABLE22062y 2.2VARIABLE22063 2.9VARIABLE22064 2.9VARIABLE22065 - 3.1VARIABLE22066 3.0VARIABLE22067  2.7VARIABLE22068 2.2VARIABLE22069 2.5VARIABLE22070Y 3.2VARIABLE22071} 2.8VARIABLE22072 2.5VARIABLE22073 2.8VARIABLE22074h 2.9VARIABLE22075 3.0VARIABLE22076 2.8VARIABLE22077 3.0VARIABLE22078  2.9VARIABLE22079 2.6VARIABLE220808 2.4VARIABLE22081.\ 2.4VARIABLE22082|H 2.7VARIABLE22083 2.7VARIABLE22084q 3.0VARIABLE22085$ 3.4VARIABLE22086^ 3.1VARIABLE22087& 2.3VARIABLE22088 3.0VARIABLE22089, 2.5VARIABLE22090? 2.6VARIABLE220910 3.0VARIABLE22092Q 2.6VARIABLE22093:{ 2.3VARIABLE22094=y 2.7VARIABLE22095b 3.0VARIABLE22096f 2.9VARIABLE22097- 2.9VARIABLE22098 2.5VARIABLE22099k 2.8VARIABLE22100o 3.3VARIABLE23101 2.7VARIABLE23102Y 3.0VARIABLE23103~} 2.9VARIABLE23104v 3.0VARIABLE23105z 3.0VARIABLE23106[h 2.5VARIABLE23107 2.9VARIABLE23108 2.5VARIABLE23109 3.6VARIABLE23110] 3.2VARIABLE23111A 2.7VARIABLE23112D8 3.0VARIABLE23113H\ 2.5VARIABLE23114H 2.8VARIABLE23115} 3.2VARIABLE23116Mq 3.0VARIABLE23117Q 3.8VARIABLE23118 2.6VARIABLE23119\& 2.2VARIABLE23120U 3.2VARIABLE23121 4Y 2.8VARIABLE23122 2.8VARIABLE231230 2.7VARIABLE231243Q 3.3VARIABLE231257{ 3.2VARIABLE23126y 2.8VARIABLE23127 3.0VARIABLE23128 2.8VARIABLE23129- 3.0VARIABLE23130 2.8VARIABLE23131 3.8VARIABLE23132 2.8VARIABLE23133 2.8VARIABLE23134Y 2.6VARIABLE23135} 3.0VARIABLE23136 3.4VARIABLE23137 3.1VARIABLE23138h 3.0VARIABLE23139 3.1VARIABLE23140 3.1VARIABLE23141 3.1VARIABLE23142 2.7VARIABLE23143 3.2VARIABLE231448 3.3VARIABLE23145\ 3.0VARIABLE23146H 2.5VARIABLE23147( 3.0VARIABLE23148q 3.4VARIABLE23149 3.0VARIABLE23150 1.4VARIABLE3/1001)& 1.4VARIABLE3/1002/ 1.3VARIABLE3/1003 1.5VARIABLE3/1004 1.4VARIABLE3/1005%0 1.7VARIABLE3/1006Q 1.4VARIABLE3/1007{ 1.5VARIABLE3/1008y 1.4VARIABLE3/1009@ 1.5VARIABLE3/1010 1.5VARIABLE3/1011d- 1.6VARIABLE3/1012h 1.4VARIABLE3/1013 1.1VARIABLE3/1014 1.2VARIABLE3/1015m 1.5VARIABLE3/1016pY 1.3VARIABLE3/1017} 1.4VARIABLE3/1018 1.7VARIABLE3/1019t 1.5VARIABLE3/1020xh 1.7VARIABLE3/1021 1.5VARIABLE3/1022 1.0VARIABLE3/1023 1.7VARIABLE3/1024 1.9VARIABLE3/1025 1.6VARIABLE3/10268 1.6VARIABLE3/1027B\ 1.5VARIABLE3/1028FH 1.4VARIABLE3/1029˵ 1.6VARIABLE3/1030q 1.6VARIABLE3/1031K 1.5VARIABLE3/1032O 1.5VARIABLE3/1033& 1.4VARIABLE3/1034 1.5VARIABLE3/1035S 1.2VARIABLE3/1036W 1.3VARIABLE3/1037Z0 1.4VARIABLE3/1038Q 1.3VARIABLE3/10391{ 1.5VARIABLE3/10405y 1.3VARIABLE3/10419 1.3VARIABLE3/1042 1.3VARIABLE3/1043- 1.6VARIABLE3/1044 1.9VARIABLE3/1045 1.4VARIABLE3/1046  1.6VARIABLE3/1047 1.4VARIABLE3/1048Y 1.5VARIABLE3/1049} 1.4VARIABLE3/1050 4.7VARIABLE3/2051 4.5VARIABLE3/2052h 4.9VARIABLE3/2053 4.0VARIABLE3/2054 4.6VARIABLE3/2055 4.5VARIABLE3/2056 4.7VARIABLE3/2057 3.3VARIABLE3/20588 4.6VARIABLE3/2059\ 3.9VARIABLE3/2060H 3.5VARIABLE3/2061 4.2VARIABLE3/2062.q 4.0VARIABLE3/2063| 4.7VARIABLE3/2064 3.6VARIABLE3/2065& 4.4VARIABLE3/2066$ 4.5VARIABLE3/2067^ 4.1VARIABLE3/2068 4.5VARIABLE3/20690 3.9VARIABLE3/2070Q 4.8VARIABLE3/2071>{ 4.0VARIABLE3/2072y 4.9VARIABLE3/2073 4.7VARIABLE3/2074: 4.3VARIABLE3/2075=- 4.4VARIABLE3/2076b 4.8VARIABLE3/2077f 5.0VARIABLE3/2078 4.5VARIABLE3/2079 3.5VARIABLE3/2080jY 3.8VARIABLE3/2081n} 3.7VARIABLE3/2082r 3.9VARIABLE3/2083 5.1VARIABLE3/2084~h 4.5VARIABLE3/2085v 4.5VARIABLE3/2086z 4.7VARIABLE3/2087[ 4.4VARIABLE3/2088 4.1VARIABLE3/2089 4.0VARIABLE3/20908 4.4VARIABLE3/2091\ 4.6VARIABLE3/2092{H 4.0VARIABLE3/2093D 3.3VARIABLE3/2094Hq 4.2VARIABLE3/2095 4.2VARIABLE3/2096} 4.2VARIABLE3/2097M& 4.3VARIABLE3/2098Q 3.0VARIABLE3/2099 4.1VARIABLE3/2100\ 6.0VARIABLE3/3101T0 5.1VARIABLE3/3102XQ 5.9VARIABLE3/3103{ 5.6VARIABLE3/3104y 5.8VARIABLE3/31053 6.6VARIABLE3/31067 4.5VARIABLE3/3107- 6.3VARIABLE3/3108 5.8VARIABLE3/3109 6.1VARIABLE3/3110 5.1VARIABLE3/3111 5.3VARIABLE3/3112Y 5.5VARIABLE3/3113} 5.0VARIABLE3/3114 5.1VARIABLE3/3115 5.3VARIABLE3/3116h 5.5VARIABLE3/3117 6.7VARIABLE3/3118 6.9VARIABLE3/3119 5.0VARIABLE3/3120 5.7VARIABLE3/3121 4.9VARIABLE3/31228 6.7VARIABLE3/3123\ 4.9VARIABLE3/3124H 5.7VARIABLE3/3125 6.0VARIABLE3/3126q 4.8VARIABLE3/3127 4.9VARIABLE3/3128 4( 5.6VARIABLE3/3129& 5.8VARIABLE3/3130 6.1VARIABLE3/3131 6.4VARIABLE3/3132) 5.6VARIABLE3/3133-0 5.1VARIABLE3/3134Q 5.6VARIABLE3/3135{ 6.1VARIABLE3/3136%y 5.6VARIABLE3/3137׌ 5.5VARIABLE3/3138 4.8VARIABLE3/3139- 5.4VARIABLE3/3140@ 5.6VARIABLE3/3141 5.1VARIABLE3/3142d 5.1VARIABLE3/3143h 5.9VARIABLE3/3144Y 5.7VARIABLE3/3145} 5.2VARIABLE3/3146l 5.0VARIABLE3/3147p 5.2VARIABLE3/3148h 5.4VARIABLE3/3149 5.1VARIABLE3/3150t 0.2VARIABLE4a1001x 0.2VARIABLE4a1002 0.2VARIABLE4a1003 0.2VARIABLE4a10048 0.2VARIABLE4a1005\ 0.4VARIABLE4a1006H 0.3VARIABLE4a1007ص 0.2VARIABLE4a1008Bq 0.2VARIABLE4a1009F 0.1VARIABLE4a1010˝ 0.2VARIABLE4a1011& 0.2VARIABLE4a1012K 0.1VARIABLE4a1013O 0.1VARIABLE4a1014 0.2VARIABLE4a10150 0.4VARIABLE4a1016Q 0.4VARIABLE4a1017V{ 0.3VARIABLE4a1018Zy 0.3VARIABLE4a1019 0.3VARIABLE4a10201 0.2VARIABLE4a10215- 0.4VARIABLE4a10229 0.2VARIABLE4a1023 0.5VARIABLE4a1024 0.2VARIABLE4a1025 0.2VARIABLE4a1026Y 0.4VARIABLE4a1027 } 0.2VARIABLE4a1028 0.2VARIABLE4a1029 0.2VARIABLE4a1030h 0.2VARIABLE4a1031 0.4VARIABLE4a1032 0.1VARIABLE4a1033 0.2VARIABLE4a1034 0.2VARIABLE4a1035 0.2VARIABLE4a10368 0.2VARIABLE4a1037\ 0.1VARIABLE4a1038H 0.2VARIABLE4a1039 0.2VARIABLE4a1040q 0.3VARIABLE4a1041 0.3VARIABLE4a1042 0.2VARIABLE4a1043.& 0.6VARIABLE4a1044| 0.4VARIABLE4a1045 0.3VARIABLE4a1046 0.2VARIABLE4a1047!0 0.2VARIABLE4a1048;Q 0.2VARIABLE4a1049{ 0.2VARIABLE4a1050y 1.4VARIABLE4a2051 1.5VARIABLE4a2052> 1.5VARIABLE4a2053- 1.3VARIABLE4a2054 1.5VARIABLE4a2055: 1.3VARIABLE4a2056= 1.6VARIABLE4a2057b 1.0VARIABLE4a2058eY 1.3VARIABLE4a2059i} 1.4VARIABLE4a2060 1.0VARIABLE4a2061j 1.5VARIABLE4a2062nh 1.0VARIABLE4a2063r 1.4VARIABLE4a2064 1.3VARIABLE4a2065~ 1.4VARIABLE4a2066v 1.5VARIABLE4a2067z 1.0VARIABLE4a20688 1.5VARIABLE4a2069\ 1.1VARIABLE4a2070H 1.8VARIABLE4a2071 1.3VARIABLE4a2072q 1.5VARIABLE4a2073{ 1.2VARIABLE4a2074D 1.3VARIABLE4a2075H& 1.4VARIABLE4a2076 1.4VARIABLE4a2077} 1.7VARIABLE4a2078M 1.5VARIABLE4a2079P0 1.0VARIABLE4a2080Q 1.1VARIABLE4a2081{ 1.0VARIABLE4a2082Ty 1.2VARIABLE4a2083X 1.6VARIABLE4a2084 1.5VARIABLE4a2085- 1.6VARIABLE4a20863 1.5VARIABLE4a20877 1.3VARIABLE4a2088 1.3VARIABLE4a2089 1.3VARIABLE4a2090 Y 1.2VARIABLE4a2091 } 1.4VARIABLE4a2092  1.2VARIABLE4a2093  1.0VARIABLE4a2094 h 1.3VARIABLE4a2095  1.2VARIABLE4a2096  1.3VARIABLE4a2097  1.3VARIABLE4a2098  1.1VARIABLE4a2099  1.3VARIABLE4a2100 8 2.5VARIABLE4a3101 \ 1.9VARIABLE4a3102 H 2.1VARIABLE4a3103  1.8VARIABLE4a3104 q 2.2VARIABLE4a3105  2.1VARIABLE4a3106 1.7VARIABLE4a3107 & 1.8VARIABLE4a3108  1.8VARIABLE4a3109 ( 2.5VARIABLE4a3110  2.0VARIABLE4a3111 0 1.9VARIABLE4a3112 Q 2.1VARIABLE4a3113 *{ 2.0VARIABLE4a3114 -y 2.4VARIABLE4a3115 2.3VARIABLE4a3116 1.8VARIABLE4a3117 %- 2.2VARIABLE4a3118 ׫ 2.3VARIABLE4a3119  1.5VARIABLE4a3120  2.3VARIABLE4a3121 @ 2.0VARIABLE4a3122 "Y 2.0VARIABLE4a3123 c} 1.8VARIABLE4a3124 g 2.1VARIABLE4a3125 1.8VARIABLE4a3126 h 1.8VARIABLE4a3127 l 1.8VARIABLE4a3128 p 2.1VARIABLE4a3129 1.6VARIABLE4a3130  1.9VARIABLE4a3131 t 2.0VARIABLE4a3132 w8 2.2VARIABLE4a3133 \ 1.5VARIABLE4a3134 H 1.4VARIABLE4a3135 b 2.3VARIABLE4a3136 q 2.4VARIABLE4a3137 1.8VARIABLE4a3138 ؝ 1.8VARIABLE4a3139 B& 2.1VARIABLE4a3140 F 2.4VARIABLE4a3141  2.3VARIABLE4a3142  1.9VARIABLE4a3143 J0 2.3VARIABLE4a3144 NQ 2.5VARIABLE4a3145 R{ 2.3VARIABLE4a3146 y 1.9VARIABLE4a3147 2.0VARIABLE4a3148 V 2.3VARIABLE4a3149 Z- 1.8VARIABLE4a3150 髦KL02SMALLESTNUMBEROFCLASSESTOBETRIED,KL,INFORMAT 2X,I2 1 KU09LARGESTNUMBEROFCLASSESTOBETRIED,KU,INFORMAT 2X,I2 5/*