Stationary Wavelet with Double Generalised Rayleigh Distribution

Statistics are mathematical tools applying scientific investigations, such as engineering and medical and biological analyses. However, statistical methods are often improved. Nowadays, statisticians try to find an accurate way to solve a problem. One of these problems is estimation parameters, which can be expressed as an inverse problem when independent variables are highly correlated. (is paper’s significant goal is to interpret the parameter estimates of double generalized Rayleigh distribution in a regressionmodel using a wavelet basis. It is difficult to use the standard version of the regressionmethods in practical terms, which is obtained using the likelihood. Since a noise level usually makes the result of estimation unstable, multicollinearity leads to various estimates. (is kind of problem estimates that features of the truth are complicated. So it is reasonable to use a mixed method that combines a fully Bayesian approach and a wavelet basis. (e usual rule for wavelet approaches is to choose a wavelet basis, where it helps to compute the wavelet coefficients, and then, these coefficients are used to remove Gaussian noise. Recovering data is typically calculated by inverting the wavelet coefficients. Some wavelet bases have been considered, which provide a shift-invariant wavelet transform, simultaneously providing improvements in smoothness, in recovering, and in squared-error performance.(e proposedmethod uses combining a penalizedmaximum likelihood approach, a penalty term, and wavelet tools. In this paper, real data are involved and modeled using double generalized Rayleigh distributions, as they are used to estimate the wavelet coefficients of the sample using numerical tools. In practical applications, wavelet approaches are recommended. (ey reduce noise levels. (is process may be useful since the noise level is often corrupted in real data, as a significant cause of most numerical estimation problems. A simulation investigation is studied using the MCMC tool to estimate the underlying features as an essential task statistics.


Introduction
Parameters' estimation, to provide an interpreted model, is often the biggest challenge in statistics since data might contain noise, blur, or both. ese kinds of problems were found in science, geophysics, engineering, and medicine.
is kind of situation received much attention from researchers over the past decade. In practical applications, the biggest challenge in estimating the unknown parameters is that real data usually contain white noise. Hence, using a pretreatment may reduce noise, where it might provide a suitable fit. More precisely, it is used in the statistical approaches of data corrupted with white noise arising from the collocation of equipment. ere are two types of statistical tools that are usually involved in processing the data. e first one is data pretreatment, which is applied to reduce the independent variable's correlation or noise level. e second is model calibration, which can be related to using Bayesian and wavelet methods. Hence, the key issues can be presented as working with many unknown features compared to the number of observations and then an ill-posed or ill-conditioned order in the model; that is, the maximum likelihood estimation is unsuitable for estimating underlying parameters. e widespread problem is to study real data collected by magnetometer or voltage reading, which are usually highly correlated. is process is needed since the sample's measured spectral characteristics may have noise levels and blur. Statistically, several established methods can be applied, such as classical thresholding approaches. Early work for studying this procedure can be found in [1,2] who introduced a new tool for removing noise (see [3,4] for explicit motivation). Bayesian approaches were studied using different probability distributions in many fields over the last century. A common practice would be to perform exponential [5], and others were applied to various density distributions. For example, the authors in [6] studied the exponential distribution and estimated their parameters, whereas those in [7] employed the Weibull distribution to estimate the parameters using censored data. Also, the authors in [8] studied the Rayleigh distribution using consorted data. Hence, the idea of this article is to combine Bayesian and wavelet methods for estimating underlying parameters. Wavelets can be powerful mathematical tools applied to reduce the impact of multicollinearity problems. Wavelet basis can be explained as a special complicated level of the Fourier transform. However, the main reason for using wavelet approaches is that it is easy to choose between different wavelet bases. Many summaries were written about this topic by several authors. For example, Mallat [9] states that a probability density function of wavelet coefficients is notably peaked and centered around at zero. Also, the algorithm of discrete wavelet transform can be found in [10]. In wavelet, the stationary basis is recommended for the reconstruction (see [11] for more details). en, the wavelets have received many comments from scientists, while several authors analyzed some real-statistical applications (see [12] for a direct result). Different approaches to the use of the wavelet can be found in [13]. Considerable details about wavelets can found in [14]. e central concept of the Bayesian approach is using the construction of theory. However, when the rules are built carefully, the model provides a good fit afterward as the estimation process. ere are several papers on Bayesian methods (see [15] who studied Bayesian approaches in the wavelet domain). Wavelet via Bayesian approaches can be studied in many articles, such as in [16]. More details about the combination of Bayesian and wavelet can be studied in [17]. Besides, using the MCMC algorithm is extracting a sample at each run of simulation from the rule. e posterior rule is more complicated for an analytic solution. e easy type of MCMC is in [18], which can be implemented to extract notation. More details about the MCMC tool can be studied in [19][20][21]. Moreover, the estimation of the unknown parameters of the double generalized Rayleigh DGRay (c j , κ j , λ j ) distributions is proposed to provide a new tool, where J � 0, 1, . . . , j − 1 for some indexes J. In practical terms, this type of investigation is sometimes called the "level-dependent" models since the distribution parameters are estimated for each level j, especially when the measurable characteristics are assumed under two or more different conditions. For example, some wavelet coefficients have defects that are close to be around zero, whereas wavelet coefficients without defects may take a form far from zero. Consider the linear inverse problem defined by with observed measurement x n×1 � X i : i � 1, . . . , n , the vector of the unknown parameters θ n×1 � θ i : i � 1, . . . , n , and errors ε n×1 . Furthermore, ε ∼ N n (0, σ 2 I n ), the noise level is usually assumed to be independent and identically distributed normally random, and n � 2 J . Consider the unknown parameters Θ defined by where K is an orthonormal matrix containing the wavelet basis. Hence, the unknown parameters Θ can be defined by their discrete wavelet transform Θ C,D � θ C 0,0 , θ D j,l : j � 0, 1, . . . , J − 1, l � 0, 1, . . . , n − 1 , and the stationary transform is used in this article. So the number of wavelet coefficients and observations is equal. Also, the wavelet coefficients of the observed data x are defined by where x C,D is the set of the wavelet coefficients of x and Θ C,D ⊂ R is also the set of Θ, where ϱ ∼ N n (0, σ 2 I n ). Leveldependent models play a significant role in wavelet applications-this procedure allows us to investigate the value of unknown parameters at each resolution j of wavelet coefficients. ere are numerous methods for specifying values of unknown parameters of the double generalized Rayleigh distributions. Moreover, the MCMC algorithms are implemented to investigate the unknown parameters from complicated or nonstandard posterior distributions [22]. In statistics, there are many tools that can be applied to estimate parameters, such as EM and MCMC algorithms. In this article, two types of methods are supposed; the first one is the posterior mean (PM), and the second is maximum a posteriori (MAP). Figure 1 illustrates the shape of the double Rayleigh distribution for different values of c. It can be seen that as c ⟶ 0, the density double Rayleigh approaches infinity, and this type of distribution can be used to fit the density of the empirical wavelet coefficients. More precisely, the wavelet coefficients are nearby the zero, which is found using the double generalized Rayleigh distribution with c � 0 and 0 < κ ≤ 0.5. In the other words, the density double Rayleigh approaches infinity as x approaches zero when κ ∈ (0, 0.5] and c � 0, which is equivalent to the summary of Mallat. is article is structured as follows: introduction to the double generalized Rayleigh distribution is explained in Section 2. All technical arguments are referred to in Sections 3 and 4. Numerical work confirming their features and simulation study to investigate estimation properties is provided in Sections 5 and 6. Section 7 gives the result of the proposed rule to real data. e final summary and conclusions are presented in Section 8.

Double Generalized Rayleigh Distribution
e generalized Rayleigh DGRay (c j , κ j , λ j ) distribution was proposed by Aykroyd et al. as a generalized distribution. ey showed the properties of the model, such as cumulative and survivor functions. Also, Aykroyd et al. [23] showed that the generalized Rayleigh distribution works well to fit data. ey also used the Bayesian approaches to estimate unknown parameters of the generalized Rayleigh distribution. In this paper, a double generalized Rayleigh distribution will be used to model the wavelet coefficients, equivalent to the density of the wavelet coefficients. Let single wavelet coefficient θ D j,l at the level j be the probability density function (pdf ) given by where |.| is the absolute value. e cumulative distribution function (cdf ) is defined by e survivor function (sf ) is given by and the failure function (hrf ) is given by where c j > 0, κ j > 0 and λ j > 0. In some indexes, J � log 2 (n). e parameters λ j and κ j are shape Mathematical Problems in Engineering parameters, and c j is a location parameter. Setting c j � 0 and κ j � 1 in (4)-(6), the results of the standard of the double Rayleigh distribution with parameter λ j are obtained.

Bayesian Approach
In statistics, Bayesian tools play important roles, where the approach has two keys. e first one is the likelihood, concocted between observation and unknown parameters, say p(x|ζ), where ζ and x are sets of underlying parameters and observations, respectively. e second key is the prior distribution, say p(ζ), and then the combining posterior distribution. Assuming the link between the model of x and the unknown of wavelet coefficients (K T Θ D ), where σ 2 is the variance of data and can be assumed by using equation (2) and the marginal likelihood given by e result of the previous integration can be found in [24]. e equivalent likelihood is defined by In addition, the posterior distribution for Θ D given x is where n J−1 , . . . , n 0 are the size of the coefficients at each level J − 1, . . . , 0. Hence, the value of κ j is suggested as 0 < κ j ≤ 0.5. e main reason for choosing the double generalized Rayleigh is that as the value of c j ⟶ 0 and |θ D j,l | ⟶ 0, the proposed distribution approaches infinity, which is followed by the saying of Mallat about the interpretation of the wavelet coefficients distribution. Clearly, equation (12) can be used to estimate the unknown parameters Θ D given x, λ J−1 , . . . , λ 0 , κ J−1 , . . . , κ 0 , and then these unknown parameters can be employed to describe the reconstruction. Hence, the unknown parameters are made up of one set, say with . . , κ 0 , τ}, and then, the previous form (12) becomes Aykroyd et al. considered gamma prior density for λ and κ with hyperprior parameters (α 1 , β 1 ) and (α 2 , β 2 ). Also, gamma distribution is proposed for τ with hyperparameters (α 3 , β 3 ), with density function en, the posterior density of the single value of θ D j,l with parameters τ, λ j , and κ j at the level j, given the data x, is given by and the joint posterior density given data, x, can be written as e hyperprior parameters τ, κ � κ 1,J−1 , . . . , κ 1,0 , κ 2,J−1 , . . . , κ 2,0 , and β � β 1,J−1 , . . . , β 1,0 , β 2,J−1 , . . . , β 2,0 can be fixed, as follows: let the expectation and variance of ω j at resolution j, say t i,j and r i,j , where i � 1, 2, 3. By solving the following equations the corresponding hyperprior parameters can be defined as

Stationary Approaches
e vital task in the wavelet approaches is to choose a basis. For more details, the interpretation of the wavelet basis is to start with two functions. e first one is scaling or father function ϕ, where the main task of this function is to compute the scaling coefficients. e other is a wavelet or mother function ψ, where it can be used to calculate the wavelet coefficients. Several wavelet bases are now available with different degrees of smoothness. However, the Haar basis is a simple version of the wavelet transform. Moreover, there are several established wavelet families demonstrated (see [25][26][27][28][29] for details). Stationary wavelet transforms (SWTs) attracted much attention for many applications over the last few years. In particular, the classical stationary wavelet transform was introduced in [30], while the authors in [31,32] applied at that time as the maximal overlap for discrete wavelet.
In 1995, Nason extended the discrete wavelet and recalled it as the "stationary." In the same year, Ronald and David [33] proposed a new tool: stationary wavelet coefficients and is sometimes referred to as "cycle spinning." In general, the SWT can be described as "fills in the gaps" between the decimated wavelet coefficients; that is, there is no missing computation between two different values of wavelet coefficients. Nason stated that this leads to an overdetermined redundancy of the original data (see the below example for more explanation). e producer gives a shiftinvariant removing noisy tool, which simultaneously shows improvements in reconstruction quality (see Ronald and Mathematical Problems in Engineering David). For example, of the SWT, the Haar wavelet is applied to the data x � x 1 , x 2 , x 3 , x 4 . e first and second sets of the scaling and detail coefficients can be computed: where θ C 1,l , θ C 0,l , θ D 1,l , and θ D 0,l are the matrices of transform at level j � 0, 1 { }, respectively. Hence, the number N of vanishing moments decreases.
is implicates that the smoothness of the corresponding shape decreases. In this paper, Daubechies father function ϕ and mother ψ with N � 8 vanishing moments are used to provide a smooth reconstruction. e plotting procedure for the stationary wavelet transform is shown in Figure 2. It can be seen that each level j has the same number of wavelet coefficients. Figure 3 shows the scaling and wavelet functions for Daubechies with N � 8 vanishing moments. Table 1 shows the wavelet coefficients for Daubechies compact, phase N � 8. Here, we present the idea of Daubechies, omitting some technical details.

Numerical Methods
e goal of the Bayesian computation is to extract a posterior sample for some unknown parameters ζ. However, computational statistics can be explained as inverse problems. Some tools can be used to make the estimation more efficient. ey include the standard version of the MCMC algorithms, Metropolis-Hastings tools, to extract a random sample from the posterior rule p(ζ|x) in (16). e procedure of the technique for parameter estimation, through the MCMC approach, can be found in [34]; for more information, see [35,36] and more recent works such as [37]. Figure 4 shows the diagram of the procedure of the proposed methods, where the procedure starts with data, which is corrupted by noise. e data are transformed to wavelet coefficients, which are used to estimate the unknown wavelet coefficients using the suggested method, and then, the underlying signal is calculated by inverting the estimation of wavelet coefficients. e main idea of the MCMC algorithms is that the parameter can take at any valued point in the parameter space Ω, say ζ i is the value point. en, at each step, MCMC creates values, say ζ 1 i , ζ 2 i , . . . , ζ r i . Each single parameter updates separately in the order that the MCMC algorithms depend on a random walk. More precisely, the general framework of the tool is defined as follows: (i) Starting with an initial value for Θ C,D � 0 and for each level j � 0, 1, . . . , 2 J− 1 , that is, for parameters, . Hence, the current value of the prior parameters is proposed with a variance parameter for each resolution j, which is chosen to obtain an acceptable convergence rate.
(2) Compute the posterior distribution in (16). Hence, all parameters are generated from the Gaussian distribution, while the current amount of the parameter is the expectation of the normal distribution with updating  Mathematical Problems in Engineering variance, which is based on the acceptance rate. It is essential to realize to pick up a random value around the current value, that is, both low and high, with variances ς 2 ζ,j chosen to depend on an acceptance rate. More precisely, choosing any valued point ω in the parameter space is accepted. e authors in [38] stated determining value is between 20% and 30% for an acceptance rate. Hence, we considered the following gamma prior density for the variance of noise σ 2 , where the starting point is computed from the finest level of the wavelet coefficients, Θ D J−1, (see Nason). Once the sample is collected from the posterior rule, the posterior mean for ζ can be calculated by and also, the posterior variance can be calculated by where K and M are the number of run and burn-in, respectively. Hence, there is an enormous method to compute the estimate point and interval. For the MAP rule, the previous procedure is changed into a simulated annealing process of Geman and Geman; this process can answer more quickly than the posterior mean. More accurately, the MAP estimate is chosen as the final iteration θ MAP � θ K . In other words, sample mean and variance can not be computed. e maximum a posteriori estimator (MAP) is defined as where K indicates the final iteration of the run of the MCMC algorithms.

Simulation
e investigation of the proposed rule is considered. en, the results are compared to some established wavelet-based methods. e authors in [39] introduced four test signals: bumps, Doppler, heavisine, and blocks. Moreover, these functions were corrupted by the independent Gaussian noise N n (0, I n ς 2 θ ). Different sizes are studied to investigate the proposed method's performance, which is n � 64 and 128, where the four test functions were simulated. Also, various wavelet bases were used: Daubechies with N � 8 applied for the test functions heavisine, Doppler, and bumps, while Haar basis was used for blocks. e starting level was j 0 � 3, as recommended in [40]. e average mean squared-error where N and K are the numbers of the data and the runs of the MCMC algorithms. Moreover, the results from the proposed method are denoted by θ k,i , i � 1, . . . , N at k-th run of MCMC algorithms. e proposed estimators were compared to various methods, such as the Bayesian wavelet thresholding (BAYES.THR) method of Abramovich and Silverman, the ABWS rule of Chipman, Kolaczyk, and McCulloch, and the BAMS rule of Vidakovic and Ruggeri. Table 2 shows the results of the simulation when decimated and the nonstationary wavelet were used. It shows the result of AMSE; for our simulation, two bases are used. e first one is the basis with zero vanishing moments, and the other is the Daubechies' wavelets with N � 8 vanishing moments. e proposed technique always gives the best reconstructions. e main interest is to improve the result of the reconstruction. is can be seen when the size of the sample is large because extensive observations contain massive information about the feature of the signal. In general, the MAP method provides a fair resolution in the test functions. However, the worst of the results is better than the other of the competed wavelet rules. e biggest problem in the MAP estimate is that the confidence intervals can not be computed because the latest sample of posterior is picked.

Application to Medical Data
e suggested method is studied and investigated to a realworld inductance plethysmography data to evaluate the excellent performance of the proposed rules, compared to the state-of-the-art methods. e Department of Anaesthesia at the Bristol Royal Infirmary collected these observations. e number of observations is 2048, equally spaced points. Readers can obtain these data within Wave resh using data (Baby-ECG). Also, the structure of the sleep state can be downloaded using data (BabySS). Figure 5 shows the plots of BabyECG and sleep state. Hence, the aim of the investigation of the BabyECG was to specify the sleep state successfully from the observations. ese data were studied and investigated by other authors (for example, [41]). e reconstruction of the unbalanced Haar approach (red line) is illustrated in Figure 6. It is not accessible to describe every moment using the unbalanced Haar method or to talk in general about the sleep state for the babies. Figures 7  and 8 show the reconstructions of the underlying feature with the MAP method using the Haar wavelet basis and Daubechies wavelet with N � 8 vanishing moments. In our reconstruction, the value of the shape κ j is set within the interval (0, 0.4). Table 3 shows the results of the simulations using the MAP and PM estimators. As the level j decreases, the value of κ increases. In contrast, the value of the parameter λ is slightly changed. Heart rate (beats per minute)

Conclusion
In this article, we show various ways in which the Bayesian rules and wavelet methods were used successfully in the practical problem. Also, a procedure for estimating the scale parameters, k and λ, of double generalized Rayleigh was estimated based on the BabyECG sample. is approach was adopted from the wavelet method for the independent level j and Bayesian approaches. Prior probability distributions for the parameters were assumed to be gamma distribution. Bayesian estimates for the points were proposed in the cases of artificial samples under the squared-error loss. e simulation studies are showing that the proposed rules worked well, and the proposed Bayesian estimate performed better than the existing state-of-the-art methods based on signal functions by reducing the AMSE. We discussed the proposed method estimates to estimate the underlying parameters. Numerical results were obtained to compare the theoretical performance results. Some points are observed from numerical results, which are summarized as follows: (i )From the results in Table 2, the suggested method process provides better excellent results for artificial data. (ii) Estimation results under the PM and MAP methods provide better estimation than the other established wavelet denoising methods according to the MSE. (iii) e use of the suggested method allows to describe the main feature of the real data, especially when observations are large.
is paper has confirmed that the wavelet approach provides attractive alternatives to other established wavelet methods, especially when underlying signals are inhomogeneous.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e author declares that there are no conflicts of interest.