under the conditions 1
N PN
i=1 YiYT
i = I and PN
i=1 Yi = 0, which provide a
unique solution. Note that the weights are kept constant in order to preserve
the local neighborhood of each data point. The most straightforward method for
computing the M-dimensional coordinates is to nd the bottom M + 1 eigen-
vectors of the sparse matrix