In a previous paper we introduced a general class of shapepreserving wavelet approximating operators (approximators) which transform cumulative distribution functions (cdf) and densities into functions of the same type. Empirical versions of these operators are used in this paper to introduce, in an unified way, shape- preserving wavelet estimators of cdf and densities, with a priori prescribed smoothness properties. We evaluate their risk for a variety of loss functions and analyze their asymptotic behavior. We study the convergence rates depending on minimal additional assumptions about the cdf/ density. These assumptions are in terms of the function belonging to certain homogeneous Besov or Triebel- Lizorkin spaces and others. As a main evaluation tool the integral p-modulus of smoothness is used