convert.wst            package:wavethresh            R Documentation

_C_o_n_v_e_r_t _a _n_o_n-_d_e_c_i_m_a_t_e_d _w_s_t _o_b_j_e_c_t _i_n_t_o _a _w_d _o_b_j_e_c_t.

_D_e_s_c_r_i_p_t_i_o_n:

     Convert a packed-ordered non-decimated wavelet transform object
     into a time-ordered non-decimated wavelet transform object.

_U_s_a_g_e:

     convert.wst(wst, ...) 

_A_r_g_u_m_e_n_t_s:

     wst: The 'wst' class object that you wish to convert.

     ...: any other arguments

_D_e_t_a_i_l_s:

     In WaveThresh3 a non-decimated wavelet transform can be ordered in
     two different ways: as a time-ordered or packet-ordered
     representation. The coefficients in the two objects are _exactly
     the same_ it is just their internal representation and ordering
     which is different. The two different representations are useful
     in different situations. The packet-ordering is useful for curve
     estimation applications and the time-ordering is useful for time
     series applications. 

     See Nason, Sapatinas and Sawczenko, 1998 for further details on
     ordering and weaving. 

     Note that the input object must be of the non-decimated type. In
     other words the type component of the input object must be
     "'station'".  Once the input object has been converted the output
     can be used with any of the functions suitable for the
     'wd.object'.

     The actual weaving permutation for shuffling coefficients from one
     representation to another is achieved by the 'getarrvec' function.

_V_a_l_u_e:

     An object of class 'wd' containing exactly the same information as
     the input object but ordered differently as a packet-ordered
     object.

_R_E_L_E_A_S_E:

     Version 3.6 Copyright Guy Nason 1997

_S_e_e _A_l_s_o:

     'convert', 'getarrvec', 'wd', 'wd.object', 'wst', 'wst.object'.

_E_x_a_m_p_l_e_s:

     #
     # Generate a sequence of 32 random normals (say) and take their
     # this-is-escaped-code{
     #
     #myrand <- wst(rnorm(32))
     #
     # Print out the result (to verify the class and type of the object)
     #
     #myrand
     #Class 'wst' : Stationary Wavelet Transform Object:
     #       ~~~  : List with 8 components with names
     #              wp Carray nlevels filter date 
     #
     #$WP and $Carray are the coefficient matrices

     #
     #Created on : Tue Sep 29 12:29:45 1998 
     #
     #summary(.):
     #----------
     #Levels:  5 
     #Length of original:  32 
     #Filter was:  Daub cmpct on least asymm N=10 
     #Boundary handling:  periodic 

     #Date:  Tue Sep 29 12:29:45 1998 
     #
     # Yep, the myrand object is of class: this-is-escaped-codenormal-bracket41bracket-normal.
     #
     # Now let's convert it to class this-is-escaped-codenormal-bracket43bracket-normal. The object
     # gets returned and, as usual in S, is printed.
     #
     #convert(myrand)
     #Class 'wd' : Discrete Wavelet Transform Object:
     #       ~~  : List with 8 components with names
     #              C D nlevels fl.dbase filter type bc date 
     #
     #$ C and $ D are LONG coefficient vectors !
     #
     #Created on : Tue Sep 29 12:29:45 1998 
     #Type of decomposition:  station 
     #
     #summary(.):
     #----------
     #Levels:  5 
     #Length of original:  32 
     #Filter was:  Daub cmpct on least asymm N=10 
     #Boundary handling:  periodic 
     #Transform type:  station 
     #Date:  Tue Sep 29 12:29:45 1998 
     #
     # The returned object is of class this-is-escaped-codenormal-bracket45bracket-normal with a
     # type of "station".
     # I.e. it has been converted successfully.

