compgrot             package:wavethresh             R Documentation

_C_o_m_p_u_t_e _e_m_p_i_r_i_c_a_l _s_h_i_f_t _f_o_r _t_i_m_e _o_r_d_e_r_e_d _n_o_n-_d_e_c_i_m_a_t_e_d _t_r_a_n_s_f_o_r_m_s.

_D_e_s_c_r_i_p_t_i_o_n:

     Computes the empirical shift required for time-ordered
     non-decimated transform coefficients to bring them into time
     order.

_U_s_a_g_e:

     compgrot(J)

_A_r_g_u_m_e_n_t_s:

       J: The 'number of levels' in the non-decimated transform where
          coefficients are to be time-aligned.

_D_e_t_a_i_l_s:

     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).

_V_a_l_u_e:

     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'.

_R_E_L_E_A_S_E:

     Version 3.6 Copyright Guy Nason 1997

_N_o_t_e:

     GROT was the shop started by Reginald Perrin. Unfortunately, GROT
     stands for "Guy ROTation".

_S_e_e _A_l_s_o:

     'wst', 'wst.object', 'wpst', 'wpst.object'.

_E_x_a_m_p_l_e_s:

     #
     # Let's see how the resolution levels have to be shifted 
     #
     #compgrot(4)
     #[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.

