compgrot | R Documentation |
Computes the empirical shift required for time-ordered non-decimated transform coefficients to bring them into time order.
compgrot(J, filter.number, family)
J |
The |
filter.number |
The wavelet filter number to be used, see |
family |
The wavelet family, see |
Time-ordered non-decimated transform coefficients when raw are not in exact time alignment due to the phase of the underlying wavelet. This function returns the shifts that are necessary to apply to each resolution level in the transform to bring back each set of time-ordered coefficients into time alignment. Note that the shifts returned are approximate shifts which work for any Daubechies wavelet. More accurate shifts can be computed using detailed knowledge of the particular wavelet used.
Each shift is "to the left". I.e. higher indexed coefficients should take the place of lower-indexed coefficients. Periodic boundaries are assumed.
This realignment is mentioned in Walden and Contreras Cristan, (1997) and Nason, Sapatinas and Sawczenko, (1998).
A vector containing the shifts that need to be applied to each scale level to return them to the correct time alignment.
There are J
entries in the vector. The first entry corresponds to the shift required for the finest level coefficients (i.e. level J-1
) and the last entry corresponds to the coarsest level (i.e. level 0). Entry j
corresponds to the shift required for scale level J-j
.
Version 3.6 Copyright Guy Nason 1997
GROT was the shop started by Reginald Perrin. Unfortunately, GROT stands for "Guy ROTation".
G P Nason
wst
, wst.object
, wpst
, wpst.object
.
#
# Let's see how the resolution levels have to be shifted
#
compgrot(4, filter.number=10, family="DaubExPhase")
#[1] 2 6 15 31
#
# In other words. Scale level 3 needs to be shifted two units.
# Scale level 2 needs to be shifted 6 units
# Scale level 1 needs to be shifted 15 units
# Scale level 0 needs to be shifted 31 units.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.