

putD(wavethresh)                             R Documentation

_P_u_t _W_a_v_e_l_e_t _C_o_e_f_f_i_c_i_e_n_t_s _I_n_t_o _W_a_v_e_l_e_t _S_t_r_u_c_t_u_r_e

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

     Makes a copy of the wd object, replaces some wavelet
     coefficients into the copy, and then returns the copy.

_U_s_a_g_e_:

     putD(wd, level, v, boundary=FALSE)

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

      wd: object of class `wd' that is to be copied and have
          wavelet coefficients replaced.

   level: integer; the level at which the replacement is to
          take place.

       v: the replacement coefficients which should be of
          the correct length.

boundary: logical; if `FALSE' then only the "real" coeffi-
          cients are replaced (and it is easy to predict the
          required length of `v', just the correct power of
          2).

          If boundary is `TRUE', then you can replace the
          boundary coefficients as well (but it is hard to
          predict the required length of `v', and the infor-
          mation has to be extracted from the `first.last'
          database component of `wd').

          `boundary' has no meaning if `wd' was obtained
          with the periodic boundary handling method (`bc').

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

     The function `accessD'obtains the wavelet coefficients
     for a particular level.  This function, `putD' replaces
     coefficients at a particular level and returns a modi-
     fied wd object reflecting the change.

     As opposed to the utility of `putC', the `putD' func-
     tion is actually quite useful.  It is fun to replace
     coefficients, since then you can dream up your own
     functions, get pictures of the wavelets etc. etc.

_V_a_l_u_e_:

     A `wd' class object containing the replaced coeffi-
     cients, see `wd.object'.

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

     Release 2.2 Copyright Guy Nason 1993

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

     `wd.object', `accessD', `draw'

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

     example(wd)

     ## Set all the wavelet coefficients to zero
     for(i in 0:(wds$nlevels - 1))
       wds <- putC(wds, level=i, v=rep(0,2^i))

     ## and now see what you've done
     plot(wds)
     matplot(x, cbind(wr(wd(y)), wr(wds)), type = 'l')

