C     CS-I                                                              CS1   1 
C     GIVEN ANY TWO OF THE FOLLOWING 4 MATRICES DETERMINE THE NONMETRIC CS1   2 
C     FIT AS A FUNCTION OF IMPOSING THE RANK ORDERS OF ONE MATRIX ON THECS1   3 
C     OTHER MATRIX - 1) AN NXN MATRIX OF COEFFICIENTS REPRESENTING THE  CS1   4 
C     SIMILARITIES AMONG N OBJECTS OR THEIR DISSIMILARITIES BASED UPON  CS1   5 
C     ONE SET OF OBSERVATIONS, 2) AN NXN MATRIX OF COEFFICIENTS REPRE-  CS1   6 
C     SENTING THE INTERRELATIONS AMONG THE SAME N OBJECTS BASED UPON AN-CS1   7 
C     OTHER SET OF OBSERVATIONS, 3) AN NXM MATRIX OF VALUES REPRESENT-  CS1   8 
C     ING THE ORTHOGONAL COORDINATES RESULTING FROM THE ANALYSIS OF THE CS1   9 
C     FIRST MATRIX WHICH GENERATES A MATRIX OF DISTANCES, AND 4) AN NXK CS1  10 
C     COORDINATE MATRIX REPRESENTING A SOLUTION FOR THE SECOND MATRIX   CS1  11 
C     WHICH GENERATES A DISTANCE MATRIX.  THE NORMALIZED PHI (THE COEF- CS1  12 
C     FICIENT OF MONOTONICITY) WHICH IS OUTPUT FOR EACH OBJECT AND OVER CS1  13 
C     ALL OBJECTS FOR ANY PAIR OF THESE MATRICES IS AN INDEX OF CONFIG- CS1  14 
C     URATIONAL OR PATTERN SIMILARITY WHICH ANSWERS SLIGHTLY DIFFERENT  CS1  15 
C     QUESTIONS DEPENDING UPON WHICH TWO MATRICES ARE INVOLVED IN THE   CS1  16 
C     COMPARISON YIELDING SIX POSSIBLE INTERPRETATIONS,I.E., 1) HOW WELLCS1  17 
C     DO THE RELATIVE DISTANCES GENERATED FROM THE THIRD MATRIX MAP INTOCS1  18 
C     THE COEFFICIENTS OF THE FIRST MATRIX, 2) SIMILARLY FOR THE FOURTH CS1  19 
C     AND THE SECOND MATRICES,MUTATIS MUTANDIS, 3) HOW WELL DOES A SOL- CS1  20 
C     UTION OF THE FIRST MATRIX (REPRESENTED BY THE THIRD MATRIX) REPRE-CS1  21 
C     SENT A REASONABLE SOLUTION OF THE SECOND MATRIX, 4) M.M. FOR THE  CS1  22 
C     FOURTH AND SECOND MATRICES, 5) HOW WELL DOES THE CONFIGURATION    CS1  23 
C     GENERATED FROM ONE SOLUTION MATCH THAT FROM ANOTHER SOLUTION (THE CS1  24 
C     THIRD AND FOURTH MATRICES), AND 6) HOW SIMILAR ARE THE PATTERNS OFCS1  25 
C     TWO MATRICES OF SIMILARITIES OR DISSIMILARITIES. THE FIRST TWO    CS1  26 
C     QUESTIONS RELATE TO THE CONVENTIONAL GOODNESS OF FIT CRITERION    CS1  27 
C     THAT IS STANDARD OUTPUT IN THE G-L PROGRAMS, SO THESE COMPARISONS CS1  28 
C     WOULD NOT USUALLY BE MADE.  IF, HOWEVER, PHI WAS NOT AVAILABLE,   CS1  29 
C     E.G., FOR THE INPUT AND OUTPUT OF A FACTOR ANALYSIS, ONE MIGHT    CS1  30 
C     WISH TO DETERMINE THIS INDEX OF PATTERN SIMILARITY.  THE SECOND   CS1  31 
C     TWO QUESTIONS, ON THE OTHER HAND, RELATE TO HOW WELL COULD ONE    CS1  32 
C     SOLUTION OF A GIVEN MATRIX BE ALSO CONSIDERED A SOLUTION TO YET   CS1  33 
C     ANOTHER MATRIX BASED UPON DIFFERENT OBSERVATIONS OF THE SAME N    CS1  34 
C     OBJECTS.  INTERPRETATIONS 3 AND 4 ARE NOT NECESSARILY SYMMETRIC,  CS1  35 
C     SINCE M MAY NOT EQUAL K OR EVEN IF M=K ONE FIT MAY BE BETTER THAN CS1  36 
C     THE OTHER FOR ITS CORRESPONDING MATRIX.  REGARDLESS OF HOW WELL A CS1  37 
C     GIVEN SOLUTION FITS ITS CORRESPONDING MATRIX WE MAY WELL BE INTER-CS1  38 
C     ESTED IN HOW WELL TWO DIFFERENT SOLUTIONS CAN BE CONSIDERED MONO- CS1  39 
C     TONIC TRANSFORMATIONS OF EACH OTHER.  AND,FINALLY,REGARDLESS OF   CS1  40 
C     WHAT KIND OF SOLUTION COULD BE OBTAINED FOR TWO MATRICES SEPARATE-CS1  41 
C     LY WE MAY WANT AN ANSWER TO THE QUESTION OF HOW CLOSE IS THE CON- CS1  42 
C     FIGURATION OF ONE MATRIX OF RELATIONS TO THAT OF ANOTHER GIVEN THECS1  43 
C     SAME OBJECTS.  ANOTHER POSSIBLE TEST THAT MIGHT BE OF INTEREST    CS1  44 
C     WOULD BE THAT OF TESTING WHETHER A GIVEN PATTERN IS LIKE SOME HY- CS1  45 
C     POTHESIZED PATTERN,E.G.,A SIMPLEX OR A CIRCUMPLEX.                CS1  46 
C                                                                       CS1  47 
C     THIS PROGRAM FOR PATTERN MATCHING WAS WRITTEN FOR THE UNIVERSITY  CS1  48 
C     OF MICHIGAN IBM-7090 COMPUTER IN FORTRAN II.  THE GUTTMAN-LINGOES CS1  49 
C     NONMETRIC SERIES RESEARCH IS SUPPORTED IN PART BY A GRANT FROM THECS1  50 
C     NATIONAL SCIENCE FOUNDATION (GS-929).  DATE = 2/15/67.            CS1  51 
C                                                                       CS1  52 
C     DECK SET-UP FOR G-L(CS-I) -                                       CS1  53 
C                                                                       CS1  54 
C        1.  SYSTEM ID CARD/S.                                          CS1  55 
C        2.  BINARY PROGRAM.                                            CS1  56 
C        3.  TITLE CARD (PUNCH A 1 IN COLUMN 1 AND ANY BCD TITLE IN COL-CS1  57 
C            UMNS 2-72 FOR OUTPUT LABELING).                            CS1  58 
C        4.  PARAMETER CARD (8 4-COLUMN FIELDS CONTAINING THE FOLLOWING CS1  59 
C            INFORMATION SERIATUM -                                     CS1  60 
C            A) N = NUMBER OF OBJECTS (.LE. 70),                        CS1  61 
C            B) M = NUMBER OF DIMENSIONS IN FIRST COORDINATE MATRIX (IF CS1  62 
C               NOT TO BE INPUT, LEAVE FIELD BLANK) .LE. 35,            CS1  63 
C            C) K = NUMBER OF DIMENSIONS IN SECOND COORDINATE MATRIX (IFCS1  64 
C               NOT TO BE INPUT, LEAVE FIELD BLANK) .LE. 35,            CS1  65 
C            D) NC = NUMBER OF COMPARISONS TO BE MADE .LE. 6,           CS1  66 
C            E) ISIMA=1 IF SIMILARITY DATA FOR MATRIX 1, OTHERWISE ZERO,CS1  67 
C            F) ISIMB=1 IF SIMILARITY DATA FOR MATRIX 2, OTHERWISE ZERO,CS1  68 
C            G) ISIMC=1 IF MATRIX 3 REPRESENTS A SCALAR SOLUTION,E.G.,  CS1  69 
C               FROM A FACTOR ANALYSIS, AND 0 IF IT REPRESENTS A SOLU-  CS1  70 
C               TION FOR EUCLIDEAN DISTANCES AND -1 IF IT REPRESENTS A  CS1  71 
C               SOLUTION FOR CITY-BLOCK DISTANCES,                      CS1  72 
C            H) ISIMD= SAME AS G) FOR MATRIX 4, I.E., 1,0, OR -1.       CS1  73 
C        5.  PARAMETER CARD (IF NC = 6, LEAVE BLANK, OTHERWISE PUNCH IN CS1  74 
C            SUCCESSIVE 4-COLUMN FIELDS THE PAIR NUMBERS CORRESPONDING  CS1  75 
C            TO THE COMPARISONS TO BE MADE,E.G.,IF NC=3 AND MATRIX 1 IS CS1  76 
C            TO BE COMPARED WITH 2, 3 WITH 4, AND 2 WITH 3, THEN PUNCH  CS1  77 
C            1 2 3 4 2 3, WHICH WOULD REQUIRE 4 MATRICES TO BE INPUT.   CS1  78 
C            THE ORDERING WITHIN PAIRS IS MEANINGFUL,I.E.,THE ORDER TO  CS1  79 
C            BE PRESERVED REFERS TO THAT OF THE FIRST PAIR MEMBER, BUT  CS1  80 
C            THE PAIRS THEMSELVES CAN BE IN ANY ORDER,E.G., 3 4 2 3 1 2 CS1  81 
C            WOULD RESULT IN THE SAME COMPARISONS.  HOWEVER, 2 1 4 3 3 2CS1  82 
C            WOULD NOT.  WHEN NC = 6, THE COMPARISON PAIRS WILL BE 1 2, CS1  83 
C            1 3,1 4,2 3,2 4, AND 3 4, IN THAT ORDER.  THE SECOND MEMBERCS1  84 
C            OF ANY PAIR MUST HAVE A HIGHER SUBSCRIPT THAN THE FIRST,   CS1  85 
C            E.G.,2 1 IS NOT LEGITIMATE.                                CS1  86 
C        6.  FORMAT CARD FOR MATRICES 1 AND 2 (DESCRIBING IN F-NOTATION CS1  87 
C            HOW DATA HAS BEEN PUNCHED).  IF BOTH NOT INPUT,LEAVE BLANK.CS1  88 
C        7.  FORMAT CARD FOR MATRICES 3 AND 4.  IF NEITHER IS TO BE IN- CS1  89 
C            PUT,LEAVE BLANK.                                           CS1  90 
C        8.  DATA (PUNCH THE UPPER-HALF OFF-DIAGONAL MATRIX FOR MATRICESCS1  91 
C            1 AND 2 AS IN SSA-I.  THE COORDINATE MATRICES ARE PUNCHED  CS1  92 
C            SUCH THAT EACH CARD (OR SET OF CARDS IN CASE 1 CARD CANNOT CS1  93 
C            ACCOMODATE THE INFORMATION) REPRESENTS ONE ROW OF THE COOR-CS1  94 
C            DINATE MATRIX FOR AN OBJECT AND THE FIELDS CORRESPOND TO   CS1  95 
C            THE DIMENSIONS.  THE RELATION MATRICES (IF ANY) MUST PRE-  CS1  96 
C            CEDE THE COORDINATE MATRICES (IF ANY) AND THE ORDER OF THE CS1  97 
C            MATRICES SHOULD BE ISOMORPHIC TO THE COMPARISON NUMBERS.   CS1  98 
C            DATA SHOULD BE PUNCHED IN COLS. 1-72, RESERVING 73-80 FOR  CS1  99 
C            OPTIONAL ID INFORMATION.                                   CS1 100 
C        9.  FOR ADDITIONAL RUNS REPEAT ITEMS 3-8.                      CS1 101 
C                                                                       CS1 102 
C     *** REFERENCES - LINGOES, J.C.  AN IBM-7090 PROGRAM FOR GUTTMAN-  CS1 103 
C                        LINGOES CONFIGURATIONAL SIMILARITY - I.  BEHAV.CS1 104 
C                        SCI., 1967,12,502-503.                         CS1 105 
C                                                                       CS1 106 
C     *** OBSERVE STARRED COMMENT CARDS FOR CONVERSIONS TO OTHER SYSTEMSCS1 107 
C                                                                       CS1 108 
      DIMENSION FMTA(18), FMTB(18), CORD(70,35), DIST(70,70), RANK(2415)CS1 109 
     1, INDI(2415), INDJ(2415), NOTIES(1208), ITIE(4900), M1(12), L1(4),CS1 110 
     2 DSTAR(70,70)                                                     CS1 111 
      COMMON CORD,DIST,RANK,INDI,INDJ,NOTIES,DSTAR                      CS1 112 
      EQUIVALENCE (DIST,ITIE)                                           CS1 113 
C                                                                       CS1 114 
C     SUBROUTINES NEEDED=RADIX SORT                                     CS1 115 
C                                                                       CS1 116 
      DEFINE FILE 13(3000,255,L,KT),14(3000,255,L,MT)                   CS1 117 
C     *** TAPE ASSIGNMENT -                                             CS1 118 
      ITAPE=5                                                           CS1 119 
      JTAPE=6                                                           CS1 120 
1     KT=1                                                              CS1 121 
      MT=1                                                              CS1 122 
C     READ TITLE, PARAMETERS, AND FORMATS                               CS1 123 
      READ (ITAPE,59)                                                   CS1 124 
      READ (ITAPE,60) N,M,K,NC,ISIMA,ISIMB,ISIMC,ISIMD                  CS1 125 
      ML=NC*2                                                           CS1 126 
      READ (ITAPE,60) (M1(J),J=1,ML)                                    CS1 127 
      READ (ITAPE,61) (FMTA(J),J=1,18)                                  CS1 128 
      READ (ITAPE,61) (FMTB(J),J=1,18)                                  CS1 129 
C     INITIALIZE                                                        CS1 130 
      NM1=N-1                                                           CS1 131 
      NEL=(N*NM1)/2                                                     CS1 132 
      NC=NC-6                                                           CS1 133 
      NN=ISIMA                                                          CS1 134 
      DO 6 L=1,2                                                        CS1 135 
      IF (NC) 2,4,2                                                     CS1 136 
2     DO 3 J=1,ML                                                       CS1 137 
      IF (M1(J)-L) 3,4,3                                                CS1 138 
3     CONTINUE                                                          CS1 139 
      L1(L)=0                                                           CS1 140 
      GO TO 6                                                           CS1 141 
4     L1(L)=1                                                           CS1 142 
C     INPUT DATA                                                        CS1 143 
      JJ=0                                                              CS1 144 
      KK=N                                                              CS1 145 
      DO 5 I=1,NM1                                                      CS1 146 
      II=JJ+1                                                           CS1 147 
      KK=KK-1                                                           CS1 148 
      JJ=JJ+KK                                                          CS1 149 
5     READ (ITAPE,FMTA) (RANK(LL),LL=II,JJ)                             CS1 150 
      WRITE (13'KT) (RANK(MM),MM=1,NEL)                                 CS1 151 
C     *** CALL ON FAST SORT ROUTINE, PREFERABLY A MACHINE CODED RADIX   CS1 152 
C         SORT WHICH STORES THE ORDERED SUBSCRIPTS AND LEAVES THE VALUESCS1 153 
C         TO BE ORDERED ALONE ***                                       CS1 154 
      CALL SORT (-1,NEL,RANK,NN,INDI)                                   CS1 155 
      WRITE (13'KT) (INDI(IP1),IP1=1,NEL)                               CS1 156 
6     NN=ISIMB                                                          CS1 157 
      NN=ISIMC                                                          CS1 158 
      MM=M                                                              CS1 159 
      DO 21 L=1,2                                                       CS1 160 
      IF (NC) 7,9,7                                                     CS1 161 
7     DO 8 J=1,ML                                                       CS1 162 
      IF (M1(J)-L-2) 8,9,8                                              CS1 163 
8     CONTINUE                                                          CS1 164 
      L1(L+2)=0                                                         CS1 165 
      GO TO 20                                                          CS1 166 
9     L1(L+2)=1                                                         CS1 167 
      DO 10 I=1,N                                                       CS1 168 
10    READ (ITAPE,FMTB) (CORD(I,JJ),JJ=1,MM)                            CS1 169 
      KK=0                                                              CS1 170 
      IF (NN) 11,16,13                                                  CS1 171 
11    NN=0                                                              CS1 172 
      DO 12 I=1,NM1                                                     CS1 173 
      IP1=I+1                                                           CS1 174 
      DO 12 J=IP1,N                                                     CS1 175 
      KK=KK+1                                                           CS1 176 
      RANK(KK)=0.                                                       CS1 177 
      DO 12 LL=1,MM                                                     CS1 178 
12    RANK(KK)=RANK(KK)+ABS(CORD(I,LL)-CORD(J,LL))                      CS1 179 
      GO TO 19                                                          CS1 180 
13    DO 15 I=1,NM1                                                     CS1 181 
      IP1=I+1                                                           CS1 182 
      DO 15 J=IP1,N                                                     CS1 183 
      KK=KK+1                                                           CS1 184 
      RANK(KK)=0.                                                       CS1 185 
      DO 14 LL=1,MM                                                     CS1 186 
14    RANK(KK)=RANK(KK)+CORD(I,LL)*CORD(J,LL)                           CS1 187 
15    CONTINUE                                                          CS1 188 
      GO TO 19                                                          CS1 189 
16    DO 18 I=1,NM1                                                     CS1 190 
      IP1=I+1                                                           CS1 191 
      DO 18 J=IP1,N                                                     CS1 192 
      KK=KK+1                                                           CS1 193 
      RANK(KK)=0.                                                       CS1 194 
      DO 17 LL=1,MM                                                     CS1 195 
17    RANK(KK)=RANK(KK)+(CORD(I,LL)-CORD(J,LL))**2                      CS1 196 
18    RANK(KK)=SQRT(RANK(KK))                                           CS1 197 
19    WRITE (13'KT) (RANK(LL),LL=1,NEL)                                 CS1 198 
C     *** CALL SORT SUBROUTINE ***                                      CS1 199 
      CALL SORT (-1,NEL,RANK,NN,INDI)                                   CS1 200 
      WRITE (13'KT) (INDI(II),II=1,NEL)                                 CS1 201 
20    NN=ISIMD                                                          CS1 202 
21    MM=K                                                              CS1 203 
      KT=1                                                              CS1 204 
C     DETERMINE IF TIES EXIST IN ANY BUT THE LAST MATRIX                CS1 205 
      DO 32 J=1,3                                                       CS1 206 
      IF (L1(J)) 22,32,22                                               CS1 207 
22    READ (13'KT) (RANK(L),L=1,NEL)                                    CS1 208 
      READ (13'KT) (INDI(I),I=1,NEL)                                    CS1 209 
      JJ=1                                                              CS1 210 
      II=-1                                                             CS1 211 
      KK=JJ                                                             CS1 212 
      NT=0                                                              CS1 213 
      NOTIES(1)=0                                                       CS1 214 
23    LL=1                                                              CS1 215 
      MM=INDI(JJ)                                                       CS1 216 
24    KK=KK+1                                                           CS1 217 
      NN=INDI(KK)                                                       CS1 218 
      IF (RANK(MM)-RANK(NN)) 26,25,26                                   CS1 219 
25    LL=LL+1                                                           CS1 220 
      IF (KK-NEL) 24,27,27                                              CS1 221 
26    IF (LL-2) 28,27,27                                                CS1 222 
27    II=II+2                                                           CS1 223 
      NT=NT+2                                                           CS1 224 
      NOTIES(II)=LL                                                     CS1 225 
      NOTIES(II+1)=JJ                                                   CS1 226 
28    JJ=KK                                                             CS1 227 
      IF (JJ-NEL) 23,29,29                                              CS1 228 
29    IF (NOTIES(1)) 31,30,31                                           CS1 229 
30    L1(J)=-1                                                          CS1 230 
      GO TO 32                                                          CS1 231 
31    L1(J)=NT                                                          CS1 232 
      WRITE (14'MT) (NOTIES(LL),LL=1,NT)                                CS1 233 
32    CONTINUE                                                          CS1 234 
      KT=1                                                              CS1 235 
      MT=1                                                              CS1 236 
C     INPUT MATRICES AND CALCULATE COEFFICIENT OF CONFIGURATIONAL       CS1 237 
C     SIMILARITY                                                        CS1 238 
      JJ=1                                                              CS1 239 
      DO 58 I=1,3                                                       CS1 240 
      NT=L1(I)                                                          CS1 241 
      IF (NT) 34,58,33                                                  CS1 242 
33    READ (14'MT) (NOTIES(L),L=1,NT)                                   CS1 243 
34    READ (13'KT) (RANK(II),II=1,NEL)                                  CS1 244 
      READ (13'KT) (INDI(KK),KK=1,NEL)                                  CS1 245 
      IP1=I+1                                                           CS1 246 
      DO 56 J=IP1,4                                                     CS1 247 
      IF (L1(J)) 35,56,35                                               CS1 248 
35    READ (13'KT) (RANK(LL),LL=1,NEL)                                  CS1 249 
      READ (13'KT) (INDJ(MM),MM=1,NEL)                                  CS1 250 
      IF (NC) 36,39,36                                                  CS1 251 
36    DO 38 K=2,ML,2                                                    CS1 252 
      IF (M1(K)-J) 38,37,38                                             CS1 253 
37    IF (M1(K-1)-I) 38,39,38                                           CS1 254 
38    CONTINUE                                                          CS1 255 
      GO TO 56                                                          CS1 256 
C     CHECK FOR TIES                                                    CS1 257 
39    IF (NT) 48,48,40                                                  CS1 258 
40    DO 47 M=1,NT,2                                                    CS1 259 
      KK=NOTIES(M)                                                      CS1 260 
      LL=NOTIES(M+1)                                                    CS1 261 
      II=LL                                                             CS1 262 
      DO 41 L=1,KK                                                      CS1 263 
      MM=INDI(II)                                                       CS1 264 
      CORD(L,1)=RANK(MM)                                                CS1 265 
      ITIE(L)=MM                                                        CS1 266 
41    II=II+1                                                           CS1 267 
      MM=KK-1                                                           CS1 268 
42    NN=0                                                              CS1 269 
      DO 44 IJ=1,MM                                                     CS1 270 
      IP1=IJ+1                                                          CS1 271 
      DO 44 L=IP1,KK                                                    CS1 272 
      IF (CORD(L,1)-CORD(IJ,1)) 43,44,44                                CS1 273 
43    BIG=CORD(IJ,1)                                                    CS1 274 
      CORD(IJ,1)=CORD(L,1)                                              CS1 275 
      CORD(L,1)=BIG                                                     CS1 276 
      NN=1                                                              CS1 277 
      II=ITIE(IJ)                                                       CS1 278 
      ITIE(IJ)=ITIE(L)                                                  CS1 279 
      ITIE(L)=II                                                        CS1 280 
44    CONTINUE                                                          CS1 281 
      IF (NN) 42,45,42                                                  CS1 282 
45    II=LL                                                             CS1 283 
      DO 46 L=1,KK                                                      CS1 284 
      INDI(II)=ITIE(L)                                                  CS1 285 
46    II=II+1                                                           CS1 286 
47    CONTINUE                                                          CS1 287 
C     COMPUTE RANK-IMAGES                                               CS1 288 
48    DO 49 K=1,NEL                                                     CS1 289 
      KK=INDI(K)                                                        CS1 290 
      LL=INDJ(K)                                                        CS1 291 
49    CORD(KK,1)=RANK(LL)                                               CS1 292 
C     COMPUTE COEFFICIENT OF MONOTONICITY                               CS1 293 
      CAL=0.                                                            CS1 294 
      BIG=0.                                                            CS1 295 
      DO 50 K=1,NEL                                                     CS1 296 
      BIG=BIG+RANK(K)**2                                                CS1 297 
50    CAL=CAL+(RANK(K)-CORD(K,1))**2                                    CS1 298 
      CAL=SQRT(1.-(1.-(CAL/(2.*BIG)))**2)                               CS1 299 
      II=0                                                              CS1 300 
      DO 51 K=1,NM1                                                     CS1 301 
      IP1=K+1                                                           CS1 302 
      DO 51 L=IP1,N                                                     CS1 303 
      II=II+1                                                           CS1 304 
      DIST(K,L)=RANK(II)                                                CS1 305 
      DIST(L,K)=DIST(K,L)                                               CS1 306 
      DSTAR(K,L)=CORD(II,1)                                             CS1 307 
51    DSTAR(L,K)=DSTAR(K,L)                                             CS1 308 
C     COMPUTE INDIVIDUAL PHI'S                                          CS1 309 
      BIG=4.*BIG                                                        CS1 310 
      DO 54 K=1,N                                                       CS1 311 
      CORD(K,1)=0.                                                      CS1 312 
      DO 53 L=1,N                                                       CS1 313 
      IF (K-L) 52,53,52                                                 CS1 314 
52    CORD(K,1)=CORD(K,1)+(DIST(K,L)-DSTAR(K,L))**2                     CS1 315 
53    CONTINUE                                                          CS1 316 
      CORD(K,1)=CORD(K,1)/BIG                                           CS1 317 
54    CONTINUE                                                          CS1 318 
      WRITE (JTAPE,59)                                                  CS1 319 
      WRITE (JTAPE,62) I,J                                              CS1 320 
      DO 55 K=1,N                                                       CS1 321 
55    WRITE (JTAPE,63) K,CORD(K,1)                                      CS1 322 
      WRITE (JTAPE,64) CAL                                              CS1 323 
56    CONTINUE                                                          CS1 324 
      KT=1                                                              CS1 325 
      DO 57 K=1,JJ                                                      CS1 326 
      READ (13'KT) (RANK(L),L=1,NEL)                                    CS1 327 
57    READ (13'KT) (INDI(M),M=1,NEL)                                    CS1 328 
      JJ=JJ+1                                                           CS1 329 
58    CONTINUE                                                          CS1 330 
      GO TO 1                                                           CS1 331 
C     *** FORMAT STATEMENTS ***                                         CS1 332 
C                                                                       CS1 333 
59    FORMAT (72H                                                       CS1 334 
     1                 )                                                CS1 335 
60    FORMAT (18I4)                                                     CS1 336 
61    FORMAT (18A4)                                                     CS1 337 
62    FORMAT (55H0NORMALIZED PHI CONTRIBUTIONS FOR EACH ITEM OF MATRICESCS1 338 
     1,I2,4H AND,I2/9H0VARIABLE,2X6HN. PHI/1H ,17(1H-))                 CS1 339 
63    FORMAT (I6,F11.5)                                                 CS1 340 
64    FORMAT (44H0GUTTMAN-LINGOES COEFFICIENT OF ALIENATION =,E13.6)    CS1 341 
      END                                                               CS1 342-
1PATTERN SIMILARITY FOR TBLS. 5.3,9.13,LINGOES                                  
   8   3   3   3   1   1   1   1                                                
   1   3   1   4   3   4                                                        
(7F3.3)                                                                         
(3F4.3)                                                                         
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
 691 614                                                                        
 591 740                                                                        
 581 704                                                                        
 598 652                                                                        
 694     623                                                                    
 611     560                                                                    
 562     453                                                                    
 596     473                                                                    
 860-334-135                                                                    
 845-416 134                                                                    
 824-394 088                                                                    
 848-322-122                                                                    
 727 551-073                                                                    
 627 505-036                                                                    
 568 523-009                                                                    
 608 366 189                                                                    
1PATTERN SIMILARITY FOR TBLS. 5.3,9.14,LINGOES                                  
   8   2   3   3   1   1   1   1                                                
   1   3   1   4   3   4                                                        
(7F3.3)                                                                         
(3F4.3)                                                                         
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
 858-328                                                                        
 849-414                                                                        
 810-412                                                                        
 825-339                                                                        
 747 561                                                                        
 637 507                                                                        
 561 488                                                                        
 619 371                                                                        
 860-334-135                                                                    
 845-416 134                                                                    
 824-394 088                                                                    
 848-322-122                                                                    
 727 551-073                                                                    
 627 505-036                                                                    
 568 523-009                                                                    
 608 366 189                                                                    
1PATTERN SIMILARITY FOR TBLS. 5.3,10.11,LINGOES                                 
   8   2   3   3   1   1   1   1                                                
   1   3   1   4   3   4                                                        
(7F3.3)                                                                         
(3F4.3)                                                                         
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
 830-396                                                                        
 818-469                                                                        
 777-470                                                                        
 798-401                                                                        
 786 500                                                                        
 672 458                                                                        
 594 444                                                                        
 647 333                                                                        
 860-334-135                                                                    
 845-416 134                                                                    
 824-394 088                                                                    
 848-322-122                                                                    
 727 551-073                                                                    
 627 505-036                                                                    
 568 523-009                                                                    
 608 366 189                                                                    
1PATTERN SIMILARITY FOR TBLS. 5.3,17.4,LINGOES                                  
   8   3   3   3   1   1   1   1                                                
   1   3   1   4   3   4                                                        
(7F3.3)                                                                         
(3F4.3)                                                                         
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
 874-258-102                                                                    
 874-370 081                                                                    
 838-358-013                                                                    
 849-281-093                                                                    
 710 625-141                                                                    
 596 529-024                                                                    
 534 544-028                                                                    
 612 423 493                                                                    
 860-334-135                                                                    
 845-416 134                                                                    
 824-394 088                                                                    
 848-322-122                                                                    
 727 551-073                                                                    
 627 505-036                                                                    
 568 523-009                                                                    
 608 366 189                                                                    
1PATTERN SIMILARITY FOR TBLS. 5.3,KAISER AND CAFFREY,LINGOES                    
   8   2   3   3   1   1   1   1                                                
   1   3   1   4   3   4                                                        
(7F3.3)                                                                         
(3F4.3)                                                                         
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
 814 420                                                                        
 803 496                                                                        
 758 494                                                                        
 787 432                                                                        
 805-482                                                                        
 680-420                                                                        
 620-444                                                                        
 645-290                                                                        
 860-334-135                                                                    
 845-416 134                                                                    
 824-394 088                                                                    
 848-322-122                                                                    
 727 551-073                                                                    
 627 505-036                                                                    
 568 523-009                                                                    
 608 366 189                                                                    
1PATTERN SIMILARITY FOR TBLS. 5.3,LINGOES,KAISER AND CAFFREY                    
   8   3   2   6   1   1   1   1                                                
   1   2   1   4   2   3   1   3   3   4   2   4                                
(7F3.3)                                                                         
(3F4.3)                                                                         
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
846805859473398301382                                                           
881826376326277415                                                              
801380319237345                                                                 
436329327365                                                                    
762730629                                                                       
583577                                                                          
539                                                                             
 860-334-135                                                                    
 845-416 134                                                                    
 824-394 088                                                                    
 848-322-122                                                                    
 727 551-073                                                                    
 627 505-036                                                                    
 568 523-009                                                                    
 608 366 189                                                                    
 814 420                                                                        
 803 496                                                                        
 758 494                                                                        
 787 432                                                                        
 805-482                                                                        
 680-420                                                                        
 620-444                                                                        
 645-290                                                                        
1PATTERN SIMILARITY FOR GUILFORD DATA - SSA-I.                                  
   8   0   3   1   0   1   0   0                                                
   2   4                                                                        
(2X35F2.2)                                                                      
(3F8.3)                                                                         
0067401912252639                                                                
00502620282638                                                                  
005239311824                                                                    
0055492522                                                                      
00462914                                                                        
004238                                                                          
0040                                                                            
-1364769 -556475  321444                                                        
-1007889 -712033  148124                                                        
  183891-1175082  -88725                                                        
 1204698 -439840 -338363                                                        
 1543582  -60803  300403                                                        
  659985  734183 -219934                                                        
 -161158 1356641  627813                                                        
-1058339  853423 -750765                                                        
