We analyze the properties of entangled random pure states of a quantum system partitioned into two smaller subsystems of dimensions N and M. Framing the problem in terms of random matrices with a fixed-trace constraint, we establish, for arbitrary N≤M, a general relation between the n-point densities and the cross moments of the eigenvalues of the reduced density matrix, i.e., the so-called Schmidt eigenvalues, and the analogous functionals of the eigenvalues of the Wishart-Laguerre ensemble of the random matrix theory. This allows us to derive explicit expressions for two-level densities, and also an exact expression for the variance of von Neumann entropy at finite N,M. Then, we focus on the moments E{Ka} of the Schmidt number K, the reciprocal of the purity. This is a random variable supported on [1,N], which quantifies the number of degrees of freedom effectively contributing to the entanglement. We derive a wealth of analytical results for E{Ka} for N=2 and 3 and arbitrary M, and also for square N=M systems by spotting for the latter a connection with the probability P(xGUEmin≥2N−−−√ξ) that the smallest eigenvalue xGUEmin of an N×N matrix belonging to the Gaussian unitary ensemble is larger than 2N−−−√ξ. As a by-product, we present an exact asymptotic expansion for P(xGUEmin≥2N−−−√ξ) for finite N as ξ→∞. Our results are corroborated by numerical simulations whenever possible, with excellent agreement.