Published online by Cambridge University Press: 01 January 2025
As a consequence of the complexity of the iterative optimum seeking procedures used by practical nonmetric multidimensional scaling algorithms, many of their computational properties have not been well understood. In particular, the following questions are of interest: (a) from the user’s point of view, are there significant differences between alternative available programs, (b) are suboptimal solutions frequently encountered, and how does this depend on the characteristics of the algorithm? Using Monte Carlo techniques to generate dissimilarity matrices with known underlying configurations, Kruskal’s M-D-SCAL, Guttman-Lingoes’ SSA-I, and Young-Torgerson’s TORSCA-9 programs were compared. It was found: (a) that differences between the solutions obtained by the algorithms were typically so small as to be of little practical importance, (b) deviant solutions were occasionally produced by each of the algorithms, but most often by M-D-SCAL, and furthermore, most frequently in one dimension. The likelihood of being trapped in a nonoptimal position is reduced by a good choice of initial configuration, and also possibly by constructing the iterative process to favor the possibility of stepping over small local depressions when far from the location of the optimum. The conclusions of this study are based on a total of 2160 scaling solutions.
This paper is based on a study which formed part of a thesis submitted to the University of Toronto in conformity with the requirements for the degree of Doctor of Philosophy. I am indebted to a large number of people for advice and criticism; these include J. C. Ogilvie, C. D. Creelman, R. S. Lockhart, and S. Nishisato. I am also grateful for a number of useful criticisms and suggestions by the referees. This research was supported by a predoctoral fellowship and by grant APA-176 to J. C. Ogilvie from the National Research Council of Canada.