Scaling positive random matrices: concentration and asymptotic convergence

It is well known that any positive matrix can be scaled to have prescribed row and column sums by multiplying its rows and columns by certain positive scaling factors. This procedure is known as matrix scaling, and has found numerous applications in operations research, economics, image processing, and machine learning. In this work, we establish the stability of matrix scaling to random bounded perturbations. Specifically, letting à ϵ R^M×N be a positive and bounded random matrix whose entries assume a certain type of independence, we provide a concentration inequality for the scaling factors of à around those of A = E[Ã]. This result is employed to study the convergence rate of the scaling factors of à to those of A, as well as the concentration of the scaled version of à around the scaled version of A in operator norm, as M, N → ∞. We demonstrate our results in several simulations.