For given counting number $n$, consider all permutations $\pi$ of {$1,\ldots,n$}, generate for every $\pi$ its Robinson-Schensted pair of standard tableaux $(P_\pi,Q_\pi)$ and average together all the $P_\pi$ (viewed as functions on $\Bbb Z_+ \times \Bbb Z_+$). Now take the limit (which exists) as $n\rightarrow\infty$ to define a function $F(m,n)$.
Clearly $F(1,1)=1$, curiously $F(1,2)=F(2,1)=e$ and according to a calculation I just made $F(2,2)=1+e^2(BesselI(0, 2)-BesselI(1, 2))\approx 6.090678724$. (Maple's convention here for Bessel functions.) [Added later:] Another calculation now shows that $F(1,3)=F(3,1)=F(2,2)-1$; nevertheless, as yet I lack a combinatorial explanation why these two presumably transcendental quantities happen to differ by 1.
My question: does my function $F$ already appear in the literature?
One can use the hook formula to express values of $F$ as infinite series. The series converge quite rapidly but grow increasingly unwieldy with $m$ and $n$. What is known (or what can one say) about closed form expressions or simplified series expressions for the values?
The values of $F(m,n)$ have probabilistic interpretations as expected stopping times for certain random processes, as follows. Generate a sequence $(a_i)$ of independent uniform samples from $[0,1]$ until the cardinality of the largest union of $m$ increasing subsequences of $(a_i)$ exceeds by $n$ the cardinality of the largest union of $m-1$ increasing subsequences.

