A bayesian wavelet shrinkage rule under LINEX loss function

This work proposes a wavelet shrinkage rule under asymmetric LINEX loss function and a mixture of a point mass function at zero and the logistic distribution as prior distribution to the wavelet coefficients in a nonparametric regression model with gaussian error. Underestimation of a significant wavelet coefficient can lead to a bad detection of features of the unknown function such as peaks, discontinuities and oscillations. It can also occur under asymmetrically distributed wavelet coefficients. Thus the proposed rule is suitable when overestimation and underestimation have asymmetric losses. Statistical properties of the rule such as squared bias, variance, frequentist and bayesian risks are obtained. Simulation studies are conducted to evaluate the performance of the rule against standard methods and an application in a real dataset involving infrared spectra is provided.


Introduction
Wavelet shrinkage methods are important tools to be applied in nonparametric regression models under wavelet basis expansion of the unknown function.In the wavelet domain, coefficients of the representation of a function are typically sparse and the few nonzero coefficients are localized at important positions of the function to be recovered, such as peaks, discontinuities and oscillations.However the observed wavelet coefficients, usually called empirical coefficients which are obtained by the application of a discrete wavelet transformation on the data, are not sparse in practice due the presence of noise.In this sense, wavelet shrinkage rules act by denoising the empirical coefficients and estimating the true wavelet coefficients of the representation of the unknown function.The concept of wavelet shrinkage was introduced in the seminal works by Donoho (1993Donoho ( , 1995)), Donoho andJohnstone (1994, 1995) and Donoho et al. (1995).
There are several wavelet shrinkage rules available in the literature, most of them are thresholding.A thresholding rule shrinks a small empirical coefficient that is less than a threshold value to exactly zero.The hard and soft thresholding rules of Donoho and Johnstone (1994) are the most famous and applied in the literature and policies to obtain the threshold value were proposed since then, such as the universal threshold by Donoho and Johnstone (1994), the Stein unbiased risk estimator of Donoho and Johnstone (1995) and the cross validation threshold by Nason (1996), among others.Bayesian wavelet shrinkage rules have also been proposed along the last decades.In fact, bayesian procedures allow the incorporation of prior information about the wavelet coefficients in the estimation process such as sparseness level and boundedness (if they are limited) through the elicitation of a prior probabilistic distribution to them.Chipman et al. (1997) proposed a mixture of two zero mean normal distributions as prior, one of them with very small variance, in the same idea of a spike and slab prior.Mixtures of a point mass function at zero and a symmetric and unimodal distribution centered at zero were proposed by Vidakovic and Ruggeri (2001) with the Laplace distribution, the Bickel distribution by Angelini and Vidakovic (2004), the double Weibull distribution by Reményi and Vidakovic (2015), the logistic distribution by Sousa (2020) and the beta distribution by Sousa et al. (2020).
Although the bayesian shrinkage rules mentioned above and others are well succeeded in many real data applications, they were built under symmetric loss functions, in particular, the quadratic loss function, which implies that overestimation and underestimation with the same magnitudes of a given wavelet coefficient have the same loss.Moreover, when the quadratic loss function is considered, mathematical advantages are gained such as differentiation and easier analytical manipulation, in addition to the fact that the Bayes rule is the posterior expected mean of the wavelet coefficient.However, when overestimation and underestimation should have different losses, the proposed shrinkage rules are not suitable and it is necessary the development of a shrinkage rule that is built under an asymmetric loss function.In this sense, this paper proposes a bayesian shrinkage rule under a point mass function at zero and the logistic distribution as prior to the wavelet coefficient and that is obtained under the so called linear-exponential (LINEX) loss function.This asymmetric loss function that was originally proposed by Varian (1975) penalizes exponentially the underestimation (overestimation) and linearly the overestimation (underestimation) according to convenient choices of its parameters.Under the wavelet estimation point of view, its application is interesting once underestimation of significant wavelet coefficients could not detect relevant features of the unknown function, such as peaks and/or discontinuities for example.Huang (2012) and Torehzadeh and Arashi (2014) proposed policies to obtain the threshold value to the soft thresholding rule under LINEX loss function with a bayesian interpretation of the estimation process, however the bayesian shrinkage rule to be proposed in this work allows prior information about the wavelet coefficients to be taken into account in the bayesian inference about them.This paper is organized as follows: the statistical model is defined in Section 2. The estimation procedure, including the proposed shrinkkage rule under LINEX loss, is described in Section 3. Section 4 is dedicated to simulation studies that were conducted to evaluate the performance of the proposed rule.An application in a real dataset about infrared spectra denoising is done in Section 5. Final considerations are in Section 6.

Statistical Model
We start with the unidimensional nonparametric regression problem, i.e, one observes n = 2 J points (x 1 , y 1 ), • • • , (x n , y n ), J ∈ N, from the model where are independent zero mean normal random errors with unknown common variance σ 2 , σ > 0. The goal is to estimate the function f from the sample {(x i , y i )} n i=1 without any assumption regarding the functional structure of f , except that it is squared integrable.
The classical approach to estimate f under the nonparametric point of view is to expand it in terms of some functional basis, such as polynomials, Fourier, splines and wavelets for example, see Takezawa (2005).
In this work we represent f in model (2.1) in terms of wavelet basis, where {ψ j,k (x) = 2 j/2 ψ(2 j x − k), j, k ∈ Z} is an orthonormal wavelet basis for L 2 (R) constructed by dilations j and translations k of a function ψ called wavelet or mother wavelet and θ j,k are wavelet coefficients that describe features of f at spatial location 2 −j k and scale 2 j or resolution level j.Note that the problem of estimating f becomes a problem of estimating the wavelet coefficients in (2.2).In fact, wavelet is a function that satisfies the admissibility condition |Ψ(ω)| 2 /|ω| < ∞, where Ψ(•) is the Fourier transformation of ψ.
The admissibility condition implies that ψ = 0, which motivates the name wavelet.There are several wavelet functions available in the literature, see Vidakovic (1999).For the development of this work, we considered the Daubechies' compactly supported orthogonal wavelet, that does not have closed form but it is shown in Figure 2.1.We can rewrite the model (2.1) in vector notation, where We apply a wavelet discrete transform (DWT) on the original data to take them to the wavelet domain.Although DWT is usually performed by fast algorithms such as the pyramidal algorithm, it is possible to represent it by a transformation matrix W with dimension n × n which is applied on both sides of (2.3).Once DWT is linear, we obtain the following model in wavelet domain where The estimation of the vector of wavelet coefficients θ will be done coefficient by coefficient under a bayesian framework.In this sense, a prior distribution is assigned to a single wavelet coefficient θ that incorporates our knowledge about its sparsity and symmetry around zero.Then, we consider the prior distribution π(•; α, τ ) based on a mixture of a point mass function at zero δ 0 (•) and a symmetric around zero logistic distribution g(•; τ ) proposed by Sousa (2020), where α ∈ (0, 1), τ > 0 and I R (•) is the indicator function on the real set R. Thus, the prior distribution has two hyperparameters to be elicited, α and τ .As discussed in Sousa (2020), these hyperparameters control the severity of the shrinkage imposed by the bayesian rule on the empirical coefficients.In fact, the empirical coefficients are shrunk more for higher values of α and smaller values of τ since the prior distribution becomes more concentrated around zero.In the same work, the author suggests values of τ ≤ 10 and the elicitation of α according to Angelini and Vidakovic (2004), where J 0 ≤ j ≤ J − 1, J 0 is the primary resolution level, J is the number of resolution levels, J = log 2 (n) and γ > 0. They also suggest that in the absence of additional information, γ = 2 can be adopted.Finally, it is necessary to estimate σ for a complete (hyper)parameter elicitation.According to Donoho and Johnstone (1994), based on the fact that much of the noise information present in the data can be obtained on the finer resolution scale, they proposed 3 Estimation Procedure The estimation of the wavelet coefficients vector θ is done by the application of the bayesian wavelet shrinkage rule δ * (•) associated to the models (2.4) and (2.5) on the empirical wavelet coefficients vector d, i.e, for a single wavelet coefficient θ, its estimator θ is The bayesian rule δ * (•) is obtained by minimizing the posterior expected loss according to a loss function specification L(δ, θ), which is done depending on the estimation problem and its costs.For symmetric losses as example, the well known quadratic loss function L(δ, θ) = (δ − θ) 2 is generally applied and its associated Bayes rule is the posterior expected mean δ * (d) = E π (θ|d), see Robert (2007) for a complete development of bayesian frameworks.In this work however, we consider wavelet coefficients estimation problem under asymmetric losses, in the sense that underestimation and overestimation of θ lead to different costs.Thus, the LINEX loss function is a very suitable function to be considered in this context.
The next subsections describe the LINEX loss function, the associated Bayes rule and its statistical properties.

LINEX Loss Function
We consider the unidimensional LINEX loss function by Varian (1975) given by for a = 0 and b > 0. The sign of a indicates the direction of the exponential and linear losses.In fact, for a > 0, the loss is exponential when δ > θ and linear when δ < θ, i.e, the loss is greater in overestimation.On the other hand, for a < 0, the loss is exponential when δ < θ and linear when δ > θ, i.e, the loss function penalizes more the underestimation.The parameter b controls the general magnitudes of the losses.

Bayesian Shrinkage Rule
The bayesian rule δ * is obtained by minimizing the posterior expected loss E π [L(δ, θ)|d].According to Zellner (1986), under the LINEX loss function we have and the bayesian rule is where log is the natural logarithm.Note that the bayesian rule under LINEX loss function does not depend on the parameter b.Finally, it is easy to show that under the model, one has that where φ(•) is the standard normal density function and the integrals involved in (3.7) can be obtained numerically.is not a shrinkage rule in the sense of wavelet shrinkage for all values of a, although the shape of the rules mimics wavelet shrinkage rules.For example, when a = 4, small values of the empirical coefficient d are not shrunk toward zero by the bayesian rule.Further, if d = 0 then δ * (0) = −0.768.The main idea of a wavelet shrinkage rule is to shrink small empirical coefficients toward zero, even in an asymmetrically way, as it is the case under asymmetric loss functions or priors for the wavelet coefficients θ.
Hopefully, we can obtain shrinkage rules for convenient choices of a, as it is shown in Figure 3.2 (b), for a ∈ {0.5, 1, 1.5, 2, 2.5}.Note that the rules shrink small empirical coefficients toward zero in an asymmetrically manner.Moreover, for empirical coefficients far from zero, the rule acts with different severity according to the sign of the coefficient.Since the LINEX loss penalizes exponentially the overestimation when a > 0, the associated shrinkage rules shrink severely empirical coefficients greater than zero to avoid overestimation of θ, given the prior belief that most of the θ's are zero, and the severity increases for higher values of a.
Although it is not shown in the paper, the behaviour of the rule is similar when a < 0, but in the opposite sense to avoid underestimation.

Statistical Properties of the Shrinkage Rule
The squared bias Bias 2 δ (θ) = [E(δ(d)) − θ] 2 of the shrinkage rules (3.6) for σ = 1, α = 0.9, τ = 1 and a ∈ {0.5, 1, 1.5, 2, 2.5} were obtained numerically and are shown in Figure 3.3 (a).In fact, the rules have small bias when θ is close to zero and the minimum squared bias is achieved when θ = 0. Since LINEX loss is exponential for δ > θ under positive values of a, i.e, the estimation procedure penalizes more overestimation, the squared bias is higher when θ is greater than zero.The shrinkage rule is built under the prior information that most of the wavelet coefficients are zero or close to zero.Then, for coefficients greater than zero, the shrinkage rule tends to underestimate them due exponential loss imposed by the LINEX loss function against overestimation.Further, the bias increases when a increases, once the loss increases for higher a values.On the other hand, the bias is small when θ is less than zero because the loss is only linear for δ < θ and then, even under the prior belief of sparsity, the estimation procedure allows a good recover of wavelet coefficients less than zero.Analogous interpretation could be made for the frequentist risks R δ (θ) = E[L(δ, θ)] of the rules (3.6) and displayed in Figure 3.3 (c).

The variances Var
of the shrinkage rules are provided in Figure 3.3 (b).Unlike the squared bias and frequentist risks, the variances have symmetric shapes and are small when θ is close to zero.They increase as |θ| increases and achieve two maxima.Finally, the variances stabilize close to one for higher values of |θ|.Further, the obtained maxima are higher as the LINEX parameter a decreases.
The general behaviours of the squared bias, variance and frequentist risk are quite similar.They are small when θ is close to zero, increase until two local maxima and then stabilize for sufficiently high values of |θ|.
The Bayes risks r δ = E π [R δ (θ)] of the shrinkage rules are available in Table 3.1.We note that the Bayes risk increases as the LINEX parameter a increases.For example, for a = 0.5, we obtain r δ = 0.1521 and for a = 2.5 we have r δ = 0.4644, which is almost three times the risk for a = 0.5.Since shrinkage rules under higher values of a shrink in a more severe way empirical coefficients greater than zero (once the loss is exponential), we have that the Bayes risks associated to more severe shrinkage rules are greater than less severe ones.

Comparison with Quadratic Loss Function
The quadratic loss function L(δ, θ) = (δ − θ) 2 is perhaps the most considered loss function in statistical decision problems due to several reasons.Under a mathematical point of view, it is more tractably analytically than other loss functions and differentiable.Further, in statistical inference, when it is dealt with an unbiased estimator, the quadratic loss function becomes the variance of the estimator, with several nice properties and a rich theory, see Schervish (1995).
One important feature of the quadratic loss function is the symmetry around θ, which means that underestimation and overestimation with the same magnitudes of a parameter θ by an estimator δ imply in the same polynomial (quadratic) loss.On the other hand, the LINEX loss function penalizes exponentially one side and linearly the other side depending on its parameters, i.e, the losses are assymmetric.In the bayesian framework, the bayesian rule under the quadratic loss function is given by the posterior expected value of θ, i.e, δ * (d) = E π (θ|d).In this case, Sousa (2020) showed that under the prior model (2.5), the shrinkage rule is given by .
The shrinkage rules under LINEX loss for a = ±1 and under quadratic loss and prior model (2.5) for α = 0.9 and τ = σ = 1 are shown in Figure 3.4(b).In fact, the shrinkage rule under quadratic loss acts symmetrically around zero and shrinks the empirical coefficients with intermediate severity in both sides.The shrinkage rule under LINEX loss with a = −1 shrinks empirical coefficients greater than zero with strong severity and the rule with a = 1 is severe for coefficients less than zero.
In general, the choice of the loss function has a great impact on the behaviour of the associated shrinkage rule.Symmetric loss functions like the quadratic assigns equals losses for underestimation and overestimation and imply in a symmetric shrinkage rule.Asymmetric loss functions like LINEX are suitable when the loss of underestimation is different to overestimation.In this case, the associated shrinkage rule is then asymmetric around zero.

Simulation Studies
Simulation studies were conducted to evaluate the performance of the shrinkage rule (3.6) and to compare it against standard and recent proposed shrinkage and thresholding rules.Wavelet coefficients were generated from a mixture of a point mass function at zero and the beta distribution with right asymmetries given by coefficient of kurtosis equals to 0.48 (scenario 1) and 1.73 (scenario 2) and two sample sizes, n = 512 and 2048.Then, empirical wavelet coefficients were obtained from model (2.4) with iid zero mean gaussian random errors and common variance σ 2 selected according to two signal to noise ratios (SNR), SNR = 3 and 9.For each replication r , we applied the proposed shrinkage rule under LINEX loss and the comparison rules to the generated empirical coefficients and calculated the mean squared error (MSE) and the mean absolute error (MAE) Table 4.1 shows the results of the scenario 1 under asymmetric distributed wavelet coefficients with kurtosis equals to 0.48.The proposed shrinkage rule had great performances for both values of sample sizes and SNR.
Actually, the SNR had the main impact on the performances of the rules and the sample size did not have significant influence on them.For SNR = 3, the shrinkage rule under asymmetric beta prior was the best one in terms of AMSE and AMAE.The proposed rule and the shrinkage rule under symmetric beta prior also worked well and had performances close to the asymmetric beta rule one.SURE and the thresholding rule by Torehzadeh and Arashi had intermediate performances and Universal and CV thresholding rules did not worked well.For example, the AMSE and AMAE of the CV thresholding rule were respectively 10.24 and 3.38 times the AMSE and AMAE of the asymmetric beta shrinkage rule.For SNR = 9, the proposed rule was the best for both performance measures, with AMSE = 0.0004 and AMAE = 0.0121 when n = 512 and AMSE = 0.0003 and AMAE = 0.0113 when n = 2048.The shrinkage rule under asymmetric beta prior and SURE thresholding rule also had good behaviour, but with performances far for the proposed rule ones.For example, the SURE thresholding rule (the second better performance) had AMSE and AMAE respectively almost 3 and 2 times of the proposed rule ones when n = 2048.On the other hand, the thresholding rule by The results of the simulation study under scenario 2 is available in Table 4.2.In fact, the proposed shrinkage rule had the best performance for the both measures in the four combinations of sample sizes and SNR.
It suggests that the rule works better for denoising higher asymmetric distributed wavelet coefficients.The symmetric and asymmetric beta shrinkage rules also had good performances and so the SURE thresholding Thus, the simulation studies showed that the proposed shrinkage rule can be considered under asymmetric distributed wavelet coefficients context as an alternative to the shrinkage rule under asymmetric beta prior.
Their performances were closed with each other and for higher kurtosis value, the proposed rule overcome the asymmetric beta shrinkage rule for both measures and all the combinations of sample size and SNR.

Application in Real Dataset
Serenoa repens or Saw Palmetto is a palm native from North America whose extracted essential oil is applied in several medical treatments, such as Benign Prostatic Hyperplasia (BPH) in men, disorder of urinary system, sexual impotence and others.When the essential oil is adulterated by the addition of substances such as olive oil, the contaminated oil losses its chemical and biological properties and becomes inactive against BPH for example.Chemists usually apply infrared (IR) spectroscopy on samples of Saw Palmetto essential oil to detect peaks of absorbance and classify a given sample in pure or adulterated.  1868.Since the goal is to obtain the peaks of absorbance, it is necessary to filter the data and remove noise that can lead to a misclassification of a peak.In this sense, the application of a wavelet shrinkage rule is suitable to denoise this data.The first step is to obtain the DWT of the data.This procedure can be applied in the other samples for better identification of peaks and noise removal.

Final Considerations
The present work proposed a wavelet shrinkage rule under LINEX loss function and a mixture of a point mass function at zero and the logistic distribution as prior distribution to the wavelet coefficients.
The application of this rule is suitable when overestimation and underestimation of the coefficients have asymmetric losses.In the wavelet domain, underestimated significant coefficients can lead to a bad detection of important features to be recovered of the unknown function such as peaks, discontinuities and oscillations.
In this sense, to assign higher losses to underestimation should be interesting.Further, the simulation studies suggested a good performance of the proposed shrinkage rule in terms of averaged mean squared error and averaged mean absolute error measures to estimate asymmetrically distributed wavelet coefficients against standard shrinkage and thresholding rules.
The magnitudes of the losses and the severity that the empirical coefficients are shrunk by the rule are well controlled according to convenient choices of the parameters of the LINEX function and the hyperparameters of the prior distribution, which is important to have flexibility in real datasets analysis.The impact of the wavelet basis in the estimation process and the proposal of other asymmetric loss functions are suggested as future works.
is the sparse vector of unknown wavelet coefficients of f and ε = W e = [ε 1 , • • • , ε n ] is the vector of random errors.Thus we can see the empirical wavelet coefficients d as noised versions of the unknown wavelet coefficients θ.Further, the random errors in the wavelet domain ε 1 , • • • , ε n remain independent zero mean normal with common variance σ 2 since the DWT is orthogonal.We will drop the subindices and consider the model d = θ + ε for a single wavelet coefficient along the text for simplicity.

. 1 )
In fact the estimator θ is a shrunk (denoised) version of the empirical coefficient d.Once we have the estimated wavelet coefficients vector θ = [ θ1 , • • • , θn ] , the function values f is then estimated by the inverse discrete wavelet transformation (IDWT) application that can be represented by the transpose matrix W ,

3. 4
(a) presents the LINEX loss functions for a = ±1 and the quadratic loss function when θ = 0.Under overestimation, i.e, when δ 0, the LINEX loss for a = −1 is considerably greater than the quadratic loss, once its loss increases exponentially.The LINEX loss for a = 1 increases only linearly and has the lowest loss.The same interpretation occurs for the underestimation case, but the roles of the LINEX rules change with each other.(a) Loss functions.(b) Shrinkage rules.

For each scenario of
asymmetry, sample size and signal to noise ratio values, 200 replications were run and the averaged mean squared error (AMSE) and the averaged mean absolute error (AMAE) were calculated as performance measures of the rules, given respectively by We compared the proposed shrinkage rule (LINEX -Logistic) against the soft thresholding rules with the threshold value chosen according to the universal thresholding (Universal) proposed byDonoho and Johnstone (1994), cross validation (CV) byNason (1996) and the Stein unbiased risk estimator (SURE) of Donoho and Johnstone (1995), the bayesian shrinkage rules under symmetric beta prior (Beta symmetric) proposed by Sousa et al. (2020) and asymmetric beta prior (Beta asymmetric) by Sousa (2021) and the soft thresholding rule under LINEX loss proposed by Torehzadeh and Arashi (2014).
rule.It should be noted that the thresholding rule under LINEX loss by Torehzadeh and Arashi did not work well in these contexts.Further, its performance decreased for SNR = 9 under both sample size values.The standard thresholding rules under Universal and CV policies had intermediate performances relative to the comparison rules, but they were overcome by the proposed rule.For example, under n = 512 and SNR = 3, the AMSE and AMAE of CV rule were about 10.34 and 3.34 times of the proposed rule ones respectively.Finally, Figure 4.1 presents the boxplots of the MSEs and MAEs of the rules in the scenario 1 (a) and (b) and scenario 2 (c) and (d), for n = 512 and SNR = 3.The boxplots of scenario 1 show the well behaviour of the shrinkage rules under beta prior and the proposed rule in both measures.Universal and CV had high MSEs and MAEs relative to the others.On the other hand, the boxplots of scenario 2 show emphasizes the atypical behaviour of the thresholding rule of Torehzadeh and Arashi, with MSEs and MAEs considerably higher than the others ones.Similar behaviours were obtained for the other values of sample size and SNR.

Figure 5 .
Figure 5.1(a) shows the IR spectra of a sample of pure Saw Palmetto oil in n = 2 10 = 1024 points.This dataset is available in the R package ChemoSpec by Hanson (2023).Actually the complete dataset consists of IR spectra of 14 samples of pure and adulterated samples in 1868 points.The particular data considered here is the IR spectra of the second oil sample and the considered sample size is the closest dyadic number that is smaller than 1868, i.e, n = 1024 = 2 log 2 (1868) .Since the goal is to obtain the peaks of absorbance, it

Figure 5 . 1 :
Figure 5.1: Infrared spectra of Saw Palmetto essential oil (a) and its empirical wavelet coefficients by resolution level after a DWT application (b).

Figure 5 .
1(b)  presents the magnitudes of the empirical wavelet coefficients of the dataset by resolution level.The considered wavelet basis was Daubechies with 10 null moments.In fact large empirical coefficients in the coarser resolution levels are associated to the positions of the absorbance peaks.On the other hand, the empirical coefficients in the finest resolution level are typically related to noise.The wavelet shrinkage rule acts mainly in these coefficients by reducing their magnitudes.The proposed shrinkage rule (3.6) was applied in the empirical coefficients and the denoised version of the IR spectra was recover by the IDWT application (3.2).The chosen parameter of the LINEX loss function was a = 1, the scale hyperparameter of the logistic density in the prior (2.5) was τ = 5, the hyperparameter α was chosen to be level dependent according to (2.6) and the estimated standard deviation was σ = 0.0004 according to (2.7).The obtained denoised IR spectra is showed in Figure5.2(a).It is possible to observe the action of denoising in the wavenumber interval [500, 700] and [800, 1024], where local noise peaks were smoothed.Figure 5.2(b) presents a plot of the empirical coefficients in the interval [−0.10, 0.10] versus estimated wavelet coefficients.This plot is interesting to see the range of shrunk coefficients, which occurred mainly in the interval of empirical coefficients [−0.05, 0.05].Once the signal-to-noise ratio of the data is high, the shrinkage rule was active only in small empirical coefficients typically localized in the finer resolution levels.

Figure 5 . 2 :
Figure 5.2: Denoised version of the infrared spectra of Saw Palmetto essential oil by the application of the shrinkage rule under LINEX loss function (a) and empirical coefficients in [−0.10, 0.10] versus estimated (shrunk) wavelet coefficients (b).

Table 4 .
Torehzadeh and Arashi did not have good performance, which suggests that it works better for low signal to 1: AMSE and AMAE of the shrinkage and thresholding rules for generated wavelet coefficients under asymmetric distribution with kurtosis equals to 0.48 -scenario 1.

Table 4 .
2: AMSE and AMAE of the shrinkage and thresholding rules for generated wavelet coefficients under asymmetric distribution with kurtosis equals to 1.73 -scenario 2.