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 (A) over tilde is an element of R-MxN 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 (A) over tilde around those of A = E[(A) over tilde]. This result is employed to study the convergence rate of the scaling factors of (A) over tilde to those of A, as well as the concentration of the scaled version of (A) over tilde around the scaled version of A in operator norm, as M, N -> infinity. We demonstrate our results in several simulations.