2
$\begingroup$

Given a non-negative real matrix $A \in R_+^{m \times n}$, how do i convert the matrix to doubly stochastic matrix i.e, each row summed to 1 and each column sums to one. In math terms,
Row Sum : $\sum_{j=1}^n A_{ij}= 1, \forall i= 1 \cdots m$
Column Sum : $\sum_{i=1}^m A_{ij}= 1, \forall j= 1 \cdots n$

Is the conversion is possible? if not, can we find a nearest matrix which is doubly stochastic matrix?

  • 0
    You got it, Learner.2012-08-08

1 Answers 1

5

There is a paper by Richard Sinkhorn: A relationship between arbitrary positive matrices and doubly stochastic matrices, The Annals of Mathematical statistics, 35 (1964), 876–879.

There he proves the following

Theorem. If $A$ is a square matrix with strictly positive entries then there are a unique doubly stochastic matrix $T_A$ and diagonal matrices $D_1$, $D_2$ such that $T_A=D_1AD_2$. The matrices $D_1$ and $D_2$ are themselves unique up to a scalar factor.

  • 0
    Thanks for your reference. Anything similar on rectangular matrix?2012-08-07