 Solution

For centred data, classical multidimensional scaling from Euclidean distances constructs the observation Gram matrixwhere $D_{ij}^{(2)}=\Vert x_i-x_j{\Vert }^2$. By part (i), its nonzero eigenvalues agree with those of $X^{T}X$, the matrix used for principal component analysis before dividing by a covariance normalization.

If $u_i=X{v}_i/\sqrt{{\lambda }_i}$ is the unit observation-space eigenvector, the corresponding classical multidimensional scaling coordinate column isThat is exactly the column of principal component scores. Consequently classical multidimensional scaling of Euclidean pairwise distances recovers the same reduced coordinates as principal component analysis, up to signs and arbitrary orthogonal transformations within equal-eigenvalue subspaces. Distances alone cannot fix translation or orientation. The equivalence presumes centred data and Euclidean distances; general dissimilarities need not yield a positive semidefinite Gram matrix or a corresponding data covariance matrix.