Description Usage Arguments Value See Also Examples
This function is used for thresholding coefficients by group (or block) according to the hard or soft thresholding rule.
1 | BlockThresh(wc, j0, hatsigma, L, qmf, thresh = "hard")
|
wc |
wavelet coefficients. |
j0 |
coarsest decomposition scale. |
hatsigma |
estimator of noise variance. |
L |
Block size (n mod L must be 0). |
qmf |
Orthonormal quadrature mirror filter. |
thresh |
'hard' or 'soft'. |
wcb
wavelet coefficient estimators.
invblock_partition
, invblock_partition
.
1 2 3 4 5 6 7 8 9 | n <- 64
x <- MakeSignal('Ramp', n)
sig <- 0.01
y <- x + rnorm(n, sd=sig)
j0 <- 1
qmf <- MakeONFilter('Daubechies',8)
wc <- FWT_PO(y, j0, qmf)
L <- 2
wcb <- BlockThresh(wc, j0, sig, L, qmf, "hard")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.