We introduce a new approach on shape preserving estimation of cumulative distribution functions and probability density functions using the wavelet methodology for multivariate de- pendent data. Our estimators preserve shape constraints such as monotonicity, positivity and integration to one,...