12
$\begingroup$

In this question Anthony Quas asks about the expected absolute value of the determinant of an $n\times n$ row stochastic matrix $A$, where the rows are independently selected from the uniform distribution on the unit $(n-1)$-dimensional simplex $x_1+\cdots+x_n=1$, $x_i\geq 0$. I can show by a messy computation that the integral over all such matrices of $(\det A)^2$ is $1/(n+1)!^{n-1}$. Is there a noncomputational reason for such a simple value?

$\endgroup$
  • $\begingroup$ Interesting. Is this the $(n-1)$st power of the factorial? This looks a lot smaller than I was expecting; or maybe I just misunderstood the random matrix result. $\endgroup$ – Anthony Quas Dec 3 '16 at 15:53
  • 1
    $\begingroup$ This is $(n+1)!$ to the $(n-1)$st power. However, this is not the average value of $(\det A)^2$. To get the average value one must divide by the volume, which is $1/(n-1)!^n$. $\endgroup$ – Richard Stanley Dec 4 '16 at 1:11
  • $\begingroup$ Great - so this gives the sort of size that I was expecting. Thanks for taking an interest in my original question. $\endgroup$ – Anthony Quas Dec 4 '16 at 3:12
  • 1
    $\begingroup$ @RichardStanley: that's not an answer to the "reason" question, but I think I see how to make a perhaps nicer computation: by computing the expectation of the $det^2$ for the i.i.d. exponential variables case (see mathoverflow.net/a/13040/31371, though I would do the handling of $Fix$ summation part differently). Then, divide it by $(n(n+1))^n$ that corresponds to the normalization (expectation of $(a_{i1}+\dots+a_{in})^2$ per row), that sends it to the simplex. It provides $\frac{(n+1)!}{(n(n+1))^n}= \frac{(n-1)!}{(n(n+1))^{n-1}}$ for the expectation; perhaps this way is less messy? $\endgroup$ – Victor Kleptsyn Dec 5 '16 at 15:40
  • 1
    $\begingroup$ @Hans: The technique is that of Exercise 5.64 in my book Enumerative Combinatorics, vol. 2. $\endgroup$ – Richard Stanley Jan 30 '18 at 23:31
7
+50
$\begingroup$

Hidden in the comment by Victor Kleptsyn is a really nice argument. Since nobody upvoted that comment yet (I just did) it's probably worth expanding it.

From a probabilistic approach it's more natural to show the equivalent statement that the expectation of $(\det A)^2$ equals $\frac{(n+1)!}{(n(n+1))^n}$. The trick is to realize the uniform measure on the simplex as the distribution of $$ (X_1,\dots,X_n) := \frac{(Y_1,\dots,Y_n)}{Y_1 + \cdots +Y_n} $$ where $(Y_i)$ are i.i.d. exponential random variables (say, of parameter 1). All you need to know is that $\mathbf{E} [Y_i] =1$ and $\mathbf{E} [Y_i^2] =2$. Moreover in this representation $(X_i)$ are independent from $S:=Y_1+\cdots+Y_n$. It follows that we have the following identity in distribution $$ B = D A $$ where $B$ is a matrix with i.i.d. exponential entries and $D$ a diagonal matrix whose entries are independent copies of $S$, with $A$ and $D$ moreover independent. Therefore $$ \mathbf{E} [\det B^2] = \mathbf{E} [\det D^2] \cdot \mathbf{E} [\det A^2] .$$ It is very easy to check that $\mathbf{E} [\det D^2] = (n(n+1))^n$, and $\mathbf{E} [\det B^2]$ can be expanded as $$ \sum_{\sigma,\tau \in \mathfrak{S}_n} \varepsilon(\sigma) \varepsilon(\tau) 2^{Fix(\sigma \tau^{-1})} = n! \sum_{\pi \in \mathfrak{S}_n} \varepsilon(\pi) 2^{Fix(\pi)} = n! \det \begin{pmatrix} 2 & 1 & \cdots & 1 \\ 1 & 2 & \ddots & \vdots \\ \vdots & \ddots & \ddots & \vdots \\ 1 & \cdots & \cdots & 2 \end{pmatrix} = (n+1)!$$ (Here $Fix$ denotes the number of fixed points.)

| cite | improve this answer | |
$\endgroup$
1
$\begingroup$

Not a full answer.

The determinant of an $n\times n$ matrix $A$ is the volume of the parallelotope spanned by its rows. That volume is equal to $n!$ times the volume of the $n$-simplex whose vertices are at the origin and at the vectors which are the rows of $A$.

When those vectors are on the standard $(n-1)$-simplex, we can compute the volume of the $n$-simplex by the formula for the volume of a pyramid: the volume of the base times the height, divided by $n$.

In this case, the base is the $(n-1)$-simplex formed by the rows of $A$, and the height is $\frac{1}{\sqrt{n}}$.

Thus the determinant of $A$ is $\frac{n! V}{n\sqrt{n}}$, where $V$ is the volume of the $(n-1)$-simplex spanned by the rows of $A$.

The square of the determinant is $(\det A)^2 = \frac{n!^2 V^2}{n^3}$.

Therefore the average value over all such $A$ of the square of the determinant is $\frac{n!^2}{n^3}E(V^2)$.

"All" that remains is to compute the expected value $E(V^2)$ of the squared volume of a subsimplex of a standard simplex. Note that the analogous question for the expected volume $E(V)$ is apparently quite difficult, see http://mathworld.wolfram.com/SimplexSimplexPicking.html

| cite | improve this answer | |
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.