The local purity of large many-body quantum systems can be studied by following a statistical mechanical approach based on a random matrix model. Restricting the analysis to the case of global pure states, this method proved to be successful and a full characterization of the statistical properties of the local purity was obtained by computing the partition function of the problem. Here we generalize these techniques to the case of global mixed states. Since the computation of the partition function is far more challenging than in the pure case, we focus on the computation of the first moments of the local purity. Finally, we establish a connection with the theory of twirling maps in quantum channels.