Numerical simulations indicate that cosmological halos display power-law radial profiles of pseudo phase-space density (PPSD), Q ≡ ρ/σ 3 , where ρ is the mass density and σ is the velocity dispersion. We tested these predictions for Q(r) using the parameters derived from the Markov Chain Monte Carlo (MCMC) analysis performed with the MAMPOSSt mass-orbit modeling code on the observed kinematics of a velocity dispersion based stack (σ v ) of 54 nearby regular clusters of galaxies from the WINGS data set. In the definition of PPSD, the density is either in total mass ρ (Q ρ ) or in galaxy number density ν (Q ν ) of three morphological classes of galaxies (ellipticals, lenticulars, and spirals), while the velocity dispersion (obtained by inversion of the Jeans equation using the MCMC parameters) is either the total (Q ρ and Q ν ) or its radial component (Q r,ρ and Q r,ν ). We find that the PPSD profiles are indeed power-law relations for nearly all MCMC parameters. The logarithmic slopes of our observed Q ρ (r) and Q r,ρ (r) for ellipticals and spirals are in excellent agreement with the predictions for particles in simulations, but slightly shallower for S0s. For Q ν (r) and Q r,ν (r), only the ellipticals have a PPSD slope matching that of particles in simulations, while the slope for spirals is much shallower, similar to that of subhalos. However, for cluster stacks based on the richness or gas temperature, the fraction of power-law PPSDs is lower (esp. Q ν ) and the Q ρ slopes are shallower, except for S0s. The observed PPSD profiles, defined using ρ rather than ν, appear to be a fundamental property of galaxy clusters. They would be imprinted during an early phase of violent relaxation for dark matter and ellipticals, and later for spirals as they move toward dynamical equilibrium in the cluster gravitational potential, while S0s are either intermediate (richness and temperature-based stacks) or a mixed class (σ v stack).