In this paper we develop wavelet methods for detecting and estimating jumps and cusps in the mean function of a non-parametric regression model. An important characteristic of the model considered here is that it allows for conditional heteroscedastic variance, a feature frequently encountered...