The usual convergence proof of the SMACOF algorithm model for least squares multidimensional scaling critically depends on the assumption of nonnegativity of the quantities to be fitted, called the pseudodistances. When this assumption is violated, erratic convergence behavior is known to occur. Three types of circumstances in which some of the pseudodistances may become negative are outlined: nonmetric multidimensional scaling with normalization on the variance, metric multidimensional scaling including an additive constant, and multidimensional scaling under the city-block distance model. A generalization of the SMACOF method is proposed to resolve the difficulty that is based on the same rationale frequently involved in robust fitting with least absolute residuals.