Statistical convex-cocompactness for mapping class groups of non-orientable surfaces

We show that a finite volume deformation retract $\mathcal{T}_{\varepsilon_t}^{-}(\mathcal{N}_g)/\mathrm{MCG}(\mathcal{N}_g)$ of the moduli space $\mathcal{M}(\mathcal{N}_g)$ of non-orientable surfaces $\mathcal{N}_g$ behaves like the convex core of $\mathcal{M}(\mathcal{N}_g)$, despite not even being quasi-convex. We then show that geodesics in the convex core leave compact regions with exponentially low probabilities, showing that the action of $\mathrm{MCG}(\mathcal{N}_g)$ on $\mathcal{T}_{\varepsilon_t}^{-}(\mathcal{N}_g)$ is statistically convex-cocompact. Combined with results of Coulon and Yang, this shows that the growth rate of orbit points under the mapping class group action is purely exponential, pseudo-Anosov elements in mapping class groups of non-orientable surfaces are exponentially generic, and the action of mapping class group on the limit set in the horofunction boundary is ergodic with respect to the Patterson-Sullivan measure. A key step of our proof relies on complexity length, developed by Dowdall and Masur, which is an alternative notion of distance on Teichm\"uller space that accounts for geodesics that spend a considerable fraction of their time in the thin part.