C     MSA2                                                              MS2   1 
C     A GENERAL NONMETRIC TECHNIQUE FOR MULTIDIMENSIONAL SCALOGRAM ANAL-MS2   2 
C     YSIS WITH SPHERICAL BOUNDARIES.  THIS RESEARCH IN NONMETRIC METH- MS2   3 
C     ODS IS SUPPORTED IN PART BY A GRANT TO L. GUTTMAN AND J.C. LINGOESMS2   4 
C     FROM THE NATIONAL SCIENCE FOUNDATION (GS-929).  PROGRAMMED FOR THEMS2   5 
C     UNIVERSITY OF MICHIGAN IBM-7090 COMPUTER IN FORTRAN-II BY LINGOES MS2   6 
C     (11/15/66).                                                       MS2   7 
C                                                                       MS2   8 
C     DECK SET-UP FOR G-L(MSA-II) -                                     MS2   9 
C                                                                       MS2  10 
C        1.   SYSTEM ID CARD/S.                                         MS2  11 
C        2.   BINARY PROGRAM.                                           MS2  12 
C        3.   TITLE CARD (PUNCH A 1 IN COLUMN 1 AND ANY BCD TITLE IN COLMS2  13 
C             -UMNS 2 TO 72, WHICH WILL BE PRINTED OUT FOR EACH PAGE OF MS2  14 
C             OUTPUT).                                                  MS2  15 
C        4.   PARAMETER CARD, 6 4-COLUMN FIELDS CONTAINING THE FOLLOW-  MS2  16 
C             ING INFORMATION SERIATUM -                                MS2  17 
C             A)  NS=NUMBER OF SUBJECTS,I.E.,ROWS OF MATRIX,            MS2  18 
C             B)  NV=NUMBER OF VARIABLES,I.E.,COLUMNS OF MATRIX,        MS2  19 
C             C)  NC=TOTAL NUMBER OF CATEGORIES OVER ALL VARIABLES      MS2  20 
C                 NC+NS .LE. 100,                                       MS2  21 
C             D)  MIND=0 IF PROGRAM IS TO DETERMINE MINIMUM DIMENSIONAL-MS2  22 
C                 ITY OF PROBLEM, OTHERWISE MIND=MINIMUM NUMBER OF DI-  MS2  23 
C                 MENSIONS DESIRED .LE. MAXD (V.I.),                    MS2  24 
C             E)  MAXD=MAXIMUM NUMBER OF DIMENSIONS DESIRED .LE. 10,    MS2  25 
C             F)  IT=0 IF NUMBER OF ITERATIONS BETWEEN COMPUTATIONS OF  MS2  26 
C                 LOSS FUNCTION =5, OTHERWISE SET TO DESIRED NUMBER (THEMS2  27 
C                 LARGER THE PROBLEM THE SMALLER IT SHOULD BE).         MS2  28 
C        5.   FORMAT CARD (DESCRIBING IN I-NOTATION WHERE DATA APPEARS  MS2  29 
C             ON THE CARDS).                                            MS2  30 
C        6.   DATA (PUNCH THE CATEGORY NUMERIC CODE CORRESPONDING TO THEMS2  31 
C             CATEGORY IN WHICH THE S FALLS. A SET OF CARDS REPRESENTS  MS2  32 
C             ONE ROW AND FIELDS ON THE CARD/S REPRESENT COLUMNS OF     MS2  33 
C             DATA MATRIX. AT LEAST ONE SUBJECT MUST FALL IN EVERY CATE-MS2  34 
C             GORY. ONLY INTEGERS (.GE.1) ARE TO BE USED).              MS2  35 
C        7.   REPEAT ITEMS 3-6 FOR ADDITIONAL RUNS.                     MS2  36 
C                                                                       MS2  37 
C     *** REFERENCES - LINGOES, J.C.  AN IBM-7090 PROGRAM FOR GUTTMAN-  MS2  38 
C                        LINGOES MULTIDIMENSIONAL SCALOGRAM ANALYSIS    MS2  39 
C                        - II.  BEHAV.SCI.,1967,12,268-270.             MS2  40 
C                      LINGOES, J.C.  THE MULTIVARIATE ANALYSIS OF QUAL-MS2  41 
C                        ITATIVE DATA.  MULT. BEHAV. RES., 1968,        MS2  42 
C                                                                       MS2  43 
      DIMENSION C(100,100),RHO(100,100),X(100,10),STR(100),FMT(18),DIST(MS2  44 
     12500),INDI(2500),INDJ(2500),NAT(100,100),MAT(100,100),NERR(50)    MS2  45 
      COMMON C,RHO,X,STR,DIST,INDI,INDJ,NERR                            MS2  46 
      EQUIVALENCE (C,NAT), (RHO,MAT)                                    MS2  47 
C                                                                       MS2  48 
C     SUBROUTINES NEEDED=EIGEN, PLOT, AND RADIX SORT                    MS2  49 
C                                                                       MS2  50 
C     TAPE ASSIGNMENTS                                                  MS2  51 
      ITAPE=5                                                           MS2  52 
      JTAPE=6                                                           MS2  53 
C     SUBROUTINE DIMENSIONING                                                   
      MD=100                                                                    
      ND=10                                                                     
      MD3=3*MD+1                                                                
      MD4=MD3+MD                                                                
      MD5=MD4+MD                                                                
      MD6=MD5+MD                                                                
      MD7=MD6+MD                                                                
C     INPUT TITLE,PARAMETERS,AND DATA                                   MS2  54 
1     READ (ITAPE,101)                                                  MS2  55 
      READ (ITAPE,102) NS,NV,NC,MIND,MAXD,IT                            MS2  56 
      READ (ITAPE,103) (FMT(J),J=1,18)                                  MS2  57 
      DO 2 I=1,NS                                                       MS2  60 
2     READ (ITAPE,FMT) (MAT(I,J),J=1,NV)                                MS2  61 
C     DETERMINE IF CAPACITY EXCEEDED                                    MS2  62 
      NC=0                                                              MS2  63 
      DO 5 J=1,NV                                                       MS2  64 
      NERR(J)=0                                                         MS2  65 
      DO 4 I=1,NS                                                       MS2  66 
      IF (NERR(J)-MAT(I,J)) 3,4,4                                       MS2  67 
3     NERR(J)=MAT(I,J)                                                  MS2  68 
4     CONTINUE                                                          MS2  69 
5     NC=NC+NERR(J)                                                     MS2  70 
      IF (NC+NS-MD) 7,7,6                                               MS2  71 
6     WRITE (JTAPE,104)                                                 MS2  72 
      GO TO 1                                                           MS2  73 
C     CONSTRUCT MATRIX FOR FINDING INITIAL CONFIGURATION                MS2  74 
7     DO 11 I=1,NS                                                      MS2  75 
      DO 8 J=1,NC                                                       MS2  76 
8     NAT(1,J)=0                                                        MS2  77 
      JJ=0                                                              MS2  78 
      DO 9 K=1,NV                                                       MS2  79 
      LL=MAT(I,K)+JJ                                                    MS2  80 
      NAT(1,LL)=1                                                       MS2  81 
9     JJ=JJ+NERR(K)                                                     MS2  82 
      DO 10 L=1,NC                                                      MS2  83 
10    MAT(I,L)=NAT(1,L)                                                 MS2  84 
11    CONTINUE                                                          MS2  85 
C     INITIALIZATION                                                    MS2  86 
      IF (IT) 13,12,13                                                  MS2  87 
12    IT=5                                                              MS2  88 
13    KK=NC+NS                                                          MS2  89 
      NSP1=NS+1                                                         MS2  90 
      NCNS=NC*NS                                                        MS2  91 
      NCNSP1=NCNS+1                                                     MS2  92 
      NSNV=NS*NV                                                        MS2  93 
      NSNVP1=NSNV+1                                                     MS2  94 
      FNSNV=NSNV                                                        MS2  95 
      FNC=NC                                                            MS2  96 
      FNS=NS                                                            MS2  97 
      FN=KK                                                             MS2  98 
      ASSIGN 35 TO N1                                                   MS2  99 
      LFACT=1                                                           MS2 100 
      ITER=100                                                          MS2 101 
      PRESTR=1.                                                         MS2 102 
      RHO(NSP1,1)=KK*(KK-1)                                             MS2 103 
      RHO(NSP1,2)=NSNVP1                                                MS2 104 
      RHO(NSP1,2)=1.-RHO(NSP1,2)/RHO(NSP1,1)                            MS2 105 
      RHO(NSP1,3)=NS*(NV+NC)+1                                          MS2 106 
      RHO(NSP1,3)=1.-RHO(NSP1,3)/RHO(NSP1,1)                            MS2 107 
      RHO(NSP1,4)=NCNSP1                                                MS2 108 
      RHO(NSP1,4)=.5-RHO(NSP1,4)/RHO(NSP1,1)                            MS2 109 
      JJ=NS-1                                                           MS2 110 
      DO 14 I=1,JJ                                                      MS2 111 
      IP1=I+1                                                           MS2 112 
      DO 14 J=IP1,NS                                                    MS2 113 
      C(I,J)=RHO(NSP1,4)                                                MS2 114 
14    C(J,I)=C(I,J)                                                     MS2 115 
      JJ=KK-1                                                           MS2 116 
      DO 15 I=NSP1,JJ                                                   MS2 117 
      IP1=I+1                                                           MS2 118 
      DO 15 J=IP1,KK                                                    MS2 119 
      C(I,J)=RHO(NSP1,4)                                                MS2 120 
15    C(J,I)=C(I,J)                                                     MS2 121 
      LL=0                                                              MS2 122 
      DO 18 I=1,NS                                                      MS2 123 
      DO 18 J=1,NC                                                      MS2 124 
      LL=LL+1                                                           MS2 125 
      JJ=NS+J                                                           MS2 126 
      IF (MAT(I,J)) 17,16,17                                            MS2 127 
16    C(I,JJ)=RHO(NSP1,3)                                               MS2 128 
      C(JJ,I)=C(I,JJ)                                                   MS2 129 
      INDI(LL)=0                                                        MS2 130 
      GO TO 18                                                          MS2 131 
17    C(I,JJ)=RHO(NSP1,2)                                               MS2 132 
      C(JJ,I)=C(I,JJ)                                                   MS2 133 
      INDI(LL)=1                                                        MS2 134 
18    CONTINUE                                                          MS2 135 
      DO 20 I=1,KK                                                      MS2 136 
      C(I,I)=FN                                                         MS2 137 
      DO 20 J=1,KK                                                      MS2 138 
      IF (I-J) 19,20,19                                                 MS2 139 
19    C(I,I)=C(I,I)-C(I,J)                                              MS2 140 
20    CONTINUE                                                          MS2 141 
C     CALL ON HOUSEHOLDER SUBROUTINE                                    MS2 142 
      JJ=MAXD+1                                                         MS2 143 
      CALL EIGEN (C,RHO,KK,STR,JJ,MD,X(1,1),X(1,4),X(1,5),X(1,6),X(1,7),MS2 144 
     1X(1,8))                                                           MS2 145 
      DO 21 J=2,JJ                                                      MS2 146 
      STR(J-1)=STR(J)                                                   MS2 147 
      DO 21 I=1,KK                                                      MS2 148 
21    X(I,J-1)=RHO(I,J)                                                 MS2 149 
C     DETERMINE IF UPPER-BOUND IS TO BE ESTIMATED                       MS2 150 
      IF (MIND) 29,22,29                                                MS2 151 
22    A0=FN/2.                                                          MS2 152 
      JJ=0                                                              MS2 153 
      DO 24 J=1,MAXD                                                    MS2 154 
      IF (STR(J)-A0) 25,23,23                                           MS2 155 
23    JJ=JJ+1                                                           MS2 156 
24    CONTINUE                                                          MS2 157 
25    IF (JJ) 27,26,27                                                  MS2 158 
26    MAXD=1                                                            MS2 159 
      GO TO 28                                                          MS2 160 
27    MAXD=JJ                                                           MS2 161 
28    MIND=MAXD                                                         MS2 162 
      LFACT=0                                                           MS2 163 
      ITER=25                                                           MS2 164 
29    M=MIND-1                                                          MS2 165 
C     LEAST-SQUARES-RANK-IMAGES ITERATIONS                              MS2 166 
30    NITER=0                                                           MS2 167 
      M=M+1                                                             MS2 168 
31    NITER=NITER+1                                                     MS2 169 
      ITCT=0                                                            MS2 170 
      M1=2                                                              MS2 171 
32    JJ=0                                                              MS2 172 
C     COMPUTE DISTANCES                                                 MS2 173 
      DO 34 I=1,NS                                                      MS2 174 
      DO 34 J=NSP1,KK                                                   MS2 175 
      JJ=JJ+1                                                           MS2 176 
      DIST(JJ)=0.                                                       MS2 177 
      DO 33 K=1,M                                                       MS2 178 
33    DIST(JJ)=DIST(JJ)+(X(I,K)-X(J,K))**2                              MS2 179 
34    DIST(JJ)=SQRT(DIST(JJ))                                           MS2 180 
      GO TO N1, (35,36)                                                 MS2 181 
C     *** A CALL ON A FAST SORT SUBROUTINE SHOULD BE INSERTED HERE (E.G.MS2 182 
C     A RADIX SORT CODED IN ASSEMBLY LANGUAGE).                         MS2 183 
35    CALL SORT (-1,NCNS,DIST,0,INDJ)                                   MS2 184 
      ASSIGN 36 TO N1                                                   MS2 185 
      GO TO 37                                                          MS2 186 
C     *** CALL SORT                                                     MS2 187 
36    CALL SORT (1,0,DIST,0,INDJ)                                       MS2 188 
C     COMPUTE ERRORS                                                    MS2 189 
37    A=0.                                                              MS2 190 
      A1=0.                                                             MS2 191 
      LL=0                                                              MS2 192 
      DO 39 J=1,NSNV                                                    MS2 193 
      A=A+DIST(J)**2                                                    MS2 194 
      JJ=INDJ(J)                                                        MS2 195 
      IF (INDI(JJ)) 39,38,39                                            MS2 196 
38    LL=LL+1                                                           MS2 197 
      MAT(LL,1)=JJ                                                      MS2 198 
      A1=A1+DIST(JJ)                                                    MS2 199 
39    CONTINUE                                                          MS2 200 
      MM=LL                                                             MS2 201 
      A0=0.                                                             MS2 202 
      DO 41 J=NSNVP1,NCNS                                               MS2 203 
      A=A+DIST(J)**2                                                    MS2 204 
      JJ=INDJ(J)                                                        MS2 205 
      IF (INDI(JJ)) 40,41,40                                            MS2 206 
40    MM=MM+1                                                           MS2 207 
      MAT(MM,1)=JJ                                                      MS2 208 
      A0=A0+DIST(JJ)                                                    MS2 209 
41    CONTINUE                                                          MS2 210 
      IF (LL) 42,69,42                                                  MS2 211 
42    NN=LL+LL                                                          MS2 212 
      GL=FLOAT(LL)                                                      MS2 213 
      A1=A1/GL                                                          MS2 214 
      A0=A0/GL                                                          MS2 215 
      ITCT=ITCT+1                                                       MS2 216 
      IF (ITCT-IT) 45,45,43                                             MS2 217 
43    M1=1                                                              MS2 218 
      STR(NITER)=0.                                                     MS2 219 
C     COMPUTE LOSS FUNCTION                                             MS2 220 
      DO 44 J=1,LL                                                      MS2 221 
      MM=MAT(J,1)                                                       MS2 222 
      JJ=J+LL                                                           MS2 223 
      II=MAT(JJ,1)                                                      MS2 224 
44    STR(NITER)=STR(NITER)+(DIST(MM)-DIST(II))**2                      MS2 225 
      STR(NITER)=STR(NITER)/A                                           MS2 226 
      IF (NITER-6) 45,64,64                                             MS2 227 
C     COMPUTE C-MATRIX                                                  MS2 228 
45    DO 46 I=1,NS                                                      MS2 229 
      DO 46 J=NSP1,KK                                                   MS2 230 
46    C(J,I)=1.                                                         MS2 231 
      DO 52 K=1,NN                                                      MS2 232 
      MM=MAT(K,1)                                                       MS2 233 
      I=1                                                               MS2 234 
      JJ=MM-NC                                                          MS2 235 
47    IF (JJ) 49,49,48                                                  MS2 236 
48    I=I+1                                                             MS2 237 
      JJ=JJ-NC                                                          MS2 238 
      GO TO 47                                                          MS2 239 
49    J=KK+JJ                                                           MS2 240 
      IF (K-LL) 50,50,51                                                MS2 241 
50    A=A0                                                              MS2 242 
      GO TO 52                                                          MS2 243 
51    A=A1                                                              MS2 244 
52    C(J,I)=A/DIST(MM)                                                 MS2 245 
      DO 53 I=NSP1,KK                                                   MS2 246 
      C(I,I)=FNC                                                        MS2 247 
      DO 53 J=1,NS                                                      MS2 248 
53    C(I,I)=C(I,I)+C(I,J)                                              MS2 249 
      DO 54 J=1,NS                                                      MS2 250 
      C(J,J)=FNS                                                        MS2 251 
      DO 54 I=NSP1,KK                                                   MS2 252 
54    C(J,J)=C(J,J)+C(I,J)                                              MS2 253 
      DO 55 I=NSP1,KK                                                   MS2 254 
      DO 55 J=1,NS                                                      MS2 255 
55    C(I,J)=1.-C(I,J)                                                  MS2 256 
C     APPLY CORRECTIONS                                                 MS2 257 
      DO 63 K=1,M                                                       MS2 258 
      DO 62 I=1,KK                                                      MS2 259 
      RHO(I,1)=X(I,K)*C(I,I)                                            MS2 260 
      IF (I-NS) 56,56,59                                                MS2 261 
56    DO 58 J=NSP1,KK                                                   MS2 262 
      IF (C(J,I)) 57,58,57                                              MS2 263 
57    RHO(I,1)=RHO(I,1)+X(J,K)*C(J,I)                                   MS2 264 
58    CONTINUE                                                          MS2 265 
      GO TO 62                                                          MS2 266 
59    DO 61 L=1,NS                                                      MS2 267 
      IF (C(I,L)) 60,61,60                                              MS2 268 
60    RHO(I,1)=RHO(I,1)+X(L,K)*C(I,L)                                   MS2 269 
61    CONTINUE                                                          MS2 270 
62    RHO(I,1)=RHO(I,1)/FN                                              MS2 271 
      DO 63 N=1,KK                                                      MS2 272 
63    X(N,K)=RHO(N,1)                                                   MS2 273 
      GO TO (31,32), M1                                                 MS2 274 
C     DETERMINE IF FINISHED                                             MS2 275 
64    IF (STR(NITER-5)-STR(NITER)-1E-06) 66,66,65                       MS2 276 
65    IF (NITER-ITER) 45,66,45                                          MS2 277 
66    IF (LFACT) 67,70,67                                               MS2 278 
67    IF (PRESTR-STR(NITER)) 1,1,68                                     MS2 279 
68    PRESTR=STR(NITER)                                                 MS2 280 
      GO TO 70                                                          MS2 281 
C     CALCULATE STATISTICS AND PRINT TITLE                              MS2 282 
69    R=1.                                                              MS2 283 
      CAL=0.                                                            MS2 284 
      GO TO 71                                                          MS2 285 
70    R=1.-GL/FNSNV                                                     MS2 286 
      CAL=SQRT(1.-(1.-STR(NITER))**2)                                   MS2 287 
71    WRITE (JTAPE,101)                                                 MS2 288 
      WRITE (JTAPE,105) M,(MM,MM=1,10)                                  MS2 289 
      A=DIST(INDJ(NSNV))                                                MS2 290 
      IF (M-1) 72,76,72                                                 MS2 291 
C     PRINCIPAL AXIS ROTATION OF FULL SPACE                             MS2 292 
72    DO 74 I=1,KK                                                      MS2 293 
      DO 74 J=I,KK                                                      MS2 294 
      C(I,J)=0.                                                         MS2 295 
      DO 73 K=1,M                                                       MS2 296 
73    C(I,J)=C(I,J)+X(I,K)*X(J,K)                                       MS2 297 
74    C(J,I)=C(I,J)                                                     MS2 298 
C     CALL ON HOUSEHOLDER SUBROUTINE                                    MS2 299 
      CALL EIGEN (C,RHO,KK,STR,M,MD,DIST(1),DIST(MD3),DIST(MD4),DIST(MD5MS2 300 
     1),DIST(MD6),DIST(MD7))                                            MS2 301 
      DO 75 K=1,M                                                       MS2 302 
      DO 75 I=1,KK                                                      MS2 303 
75    X(I,K)=RHO(I,K)                                                   MS2 304 
      GO TO 78                                                          MS2 305 
76    DO 77 I=1,KK                                                      MS2 306 
77    RHO(I,1)=X(I,1)                                                   MS2 307 
C     SCALE ROTATED VECTORS AND PRINT                                   MS2 308 
78    A0=0.                                                             MS2 309 
      DO 80 K=1,M                                                       MS2 310 
      C(1,K)=RHO(1,K)                                                   MS2 311 
      C(2,K)=RHO(1,K)                                                   MS2 312 
      DO 79 I=2,KK                                                      MS2 313 
      C(1,K)=AMIN1(C(1,K),RHO(I,K))                                     MS2 314 
79    C(2,K)=AMAX1(C(2,K),RHO(I,K))                                     MS2 315 
80    A0=AMAX1(A0,(C(2,K)-C(1,K)))                                      MS2 316 
      A0=2./A0                                                          MS2 317 
      DO 81 K=1,M                                                       MS2 318 
      DO 81 I=1,KK                                                      MS2 319 
81    RHO(I,K)=A0*(RHO(I,K)-C(1,K))-1.                                  MS2 320 
      DO 82 I=1,NS                                                      MS2 321 
82    WRITE (JTAPE,106) I,(RHO(I,K),K=1,M)                              MS2 322 
      WRITE (JTAPE,110)                                                 MS2 323 
      DO 83 K=NSP1,KK                                                   MS2 324 
83    WRITE (JTAPE,106) K,(RHO(K,J),J=1,M)                              MS2 325 
      A=A0*A                                                            MS2 326 
      WRITE (JTAPE,107) R,CAL,A,NITER                                   MS2 327 
C     REPRODUCE SCORE MATRIX AND PRINT                                  MS2 328 
      WRITE (JTAPE,108) (MM,MM=1,30)                                    MS2 329 
      IF (M-1) 84,87,84                                                 MS2 330 
C     COMPUTE DISTANCES                                                 MS2 331 
84    JJ=0                                                              MS2 332 
      DO 86 I=1,NS                                                      MS2 333 
      DO 86 J=NSP1,KK                                                   MS2 334 
      JJ=JJ+1                                                           MS2 335 
      DIST(JJ)=0.                                                       MS2 336 
      DO 85 K=1,M                                                       MS2 337 
85    DIST(JJ)=DIST(JJ)+(X(I,K)-X(J,K))**2                              MS2 338 
86    DIST(JJ)=SQRT(DIST(JJ))                                           MS2 339 
87    JJ=0                                                              MS2 340 
      DO 91 I=1,NS                                                      MS2 341 
      DO 90 J=1,NV                                                      MS2 342 
      II=JJ+1                                                           MS2 343 
      JJ=JJ+NERR(J)                                                     MS2 344 
      A0=10000.                                                         MS2 345 
      LL=0                                                              MS2 346 
      DO 89 K=II,JJ                                                     MS2 347 
      LL=LL+1                                                           MS2 348 
      IF (DIST(K)-A0) 88,89,89                                          MS2 349 
88    A0=DIST(K)                                                        MS2 350 
      MM=LL                                                             MS2 351 
89    CONTINUE                                                          MS2 352 
90    NAT(J,1)=MM                                                       MS2 353 
91    WRITE (JTAPE,109) I,(NAT(L,1),L=1,NV)                             MS2 354 
      IF (M-1) 92,99,92                                                 MS2 355 
C     CALL PLOT SUBROUTINE                                              MS2 356 
92    CALL PLOT (RHO,C(1,1),C(1,2),C(1,3),C(1,4),C(1,5),KK,M,MD,ND)     MS2 357 
      IF (LFACT) 99,93,99                                               MS2 358 
C     IF GOING DOWN - DROP VECTORS                                      MS2 359 
93    JJ=1                                                              MS2 360 
      DO 95 J=2,M                                                       MS2 361 
      A1=STR(J)/STR(1)                                                  MS2 362 
      IF (A1-.3) 96,96,94                                               MS2 363 
94    JJ=JJ+1                                                           MS2 364 
95    CONTINUE                                                          MS2 365 
96    IF (M-JJ) 97,1,97                                                 MS2 366 
97    M=JJ-1                                                            MS2 367 
      MAXD=JJ                                                           MS2 368 
      IF (ITER-100) 98,30,98                                            MS2 369 
98    ITER=ITER+25                                                      MS2 370 
      GO TO 30                                                          MS2 371 
C     IF FINISHED - DO NEXT PROBLEM, OTHERWISE REPEAT WITH NEW M        MS2 372 
99    IF (M-MAXD) 100,1,100                                             MS2 373 
100   IF (1E-04-CAL) 30,1,1                                             MS2 374 
C     *** FORMAT STATEMENTS ***                                         MS2 375 
C                                                                       MS2 376 
101   FORMAT (72H                                                       MS2 377 
     1                 )                                                MS2 378 
102   FORMAT (6I4)                                                      MS2 379 
103   FORMAT (18A4)                                                     MS2 380 
104   FORMAT (55H0PROGRAM CAPACITY EXCEEDED.  NC + NS .G. 100. CHECK NC.MS2 381 
     1)                                                                 MS2 382 
105   FORMAT (32H0G-L(MSA-II) COORDINATES FOR M =,I3,1H./10H0DIMENSION,1MS2 383 
     10I10/1H ,130(1H-)/9H0 PERSON )                                    MS2 384 
106   FORMAT (I5,5X2P10F10.3)                                           MS2 385 
107   FORMAT (33H0COEFFICIENT OF REPRODUCIBILITY =F6.3/28H COEFFICIENT OMS2 386 
     1F ALIENATION =E12.6/22H RADIUS OF INCLUSION =2PF8.3,4H FOR,I4,12H MS2 387 
     2ITERATIONS.)                                                      MS2 388 
108   FORMAT (1H1,51X30HMSA-II REPRODUCED SCORE MATRIX/58X17HV A R I A BMS2 389 
     1 L E S/8H SUBJECT,2X30I4/1H ,131(1H-))                            MS2 390 
109   FORMAT (I5,5X30I4)                                                MS2 391 
110   FORMAT (9H0CATEGORY)                                              MS2 392 
      END                                                               MS2 393-
1MSA-II OF BELL-SIRJAMAKI STUDY OF TYPES OF GROUPS.                             
   7   5   1   2   3                                                            
(7I1)                                                                           
11122                                                                           
22222                                                                           
11211                                                                           
42423                                                                           
44423                                                                           
33312                                                                           
23322                                                                           
1RADEX - ROTARY AND EXPANSIVE DICHOTOMOUS ITEMS                                 
  36   8  16   3   3                                                            
(10I1)                                                                          
22111122                                                                        
21111122                                                                        
11111122                                                                        
11111222                                                                        
11112222                                                                        
11122222                                                                        
11222222                                                                        
12222222                                                                        
22222222                                                                        
22222122                                                                        
22221122                                                                        
22211122                                                                        
22111121                                                                        
21111121                                                                        
11111121                                                                        
11111221                                                                        
11112221                                                                        
11122221                                                                        
11222221                                                                        
12222221                                                                        
22222221                                                                        
22222121                                                                        
22221121                                                                        
22211121                                                                        
22111111                                                                        
21111111                                                                        
11111111                                                                        
11111211                                                                        
11112211                                                                        
11122211                                                                        
11222211                                                                        
12222211                                                                        
22222211                                                                        
22222111                                                                        
22221111                                                                        
22211111                                                                        
1MSA-II OF HALL AND LINDZEY THEORISTS.                                          
  17  18  54   4   4   3                                                        
(18I1)                                                                          
333213333212321321                                                      FREUD   
323333333332223333                                                      MURPHY  
112113312111111211                                                      EYSENCK 
121113313313111311                                                      SHELDON 
311223211313231313                                                      ALLPORT 
311121211321331311                                                      GOLDSTEI
123231133111112232                                                      D + M   
223233322112132313                                                      CATTELL 
321112333233233231                                                      ADLER   
321121112322332111                                                      ROGERS  
332122122222233131                                                      HORNEY  
322122222222223131                                                      FROMM   
332113321312221312                                                      JUNG    
322322223332333232                                                      SULLIVAN
321123213321331211                                                      ANGYAL  
312122111133323123                                                      LEWIN   
332113233332322333                                                      MURRAY  
1MSA-II OF PERSONALITY THEORISTS (IPSATIZED SCORES-HALL+LINDZEY,TBL. 1).        
  17  18  54   4   4   3                                                        
(18I1)                                                                          
233212333112311321                                                              
213332232231112223                                                              
113223323121112322                                                              
122223323323111312                                                              
211223211313231313                                                              
312231221321331312                                                              
123331133122112232                                                              
123232321112132313                                                              
221111332133223231                                                              
321231112332332112                                                              
332131122122233131                                                              
323131232132213131                                                              
332113321312221312                                                              
211321113231323132                                                              
321223213321331211                                                              
312222111133323123                                                              
232112132231311233                                                              
