cthresh               package:cthresh               R Documentation

_E_s_t_i_m_a_t_e _r_e_a_l _s_i_g_n_a_l _u_s_i_n_g _c_o_m_p_l_e_x-_v_a_l_u_e_d _w_a_v_e_l_e_t_s

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

     Implements the multiwavelet style and empirical Bayes shrinkage
     procedures described in Barber & Nason (2004)

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

     cthresh(data, j0 = 3, dwwt, dev = madmad, rule = "hard", 
         filter.number = 3.1, family = "LinaMayrand", plotfn = FALSE,  
         TI = FALSE, details = FALSE, policy = "mws", code = "NAG", tol = 0.01)

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

    data: The data to be analysed. This should be real-valued and of
          length a power of two.

      j0: Primary resolution level; no thresholding is done below this
          level.

    dwwt: description to come

     dev: A function to be used to estimate the noise level of the
          data. The function supplied must return a value of spread on
          the variance scale (i.e. not standard deviation) such as the
          var() function. A popular, useful and robust alternative is
          the madmad function.

    rule: The type of thresholding done. If policy = "mws", available
          rules are "hard" or "soft"; if policy = "ebayes", then rule
          can be "hard", "soft" or "mean".

filter.number, family: These parameters specify the wavelet used. See
          'filter.select' for details. 

          Also, if filter.number = 5, estimation is done with all the
          complex-valued wavelets with 5 vanishing moments and the
          results averaged. If filter.number = 0, then he averaging is
          over all available complex-valued wavelets.

  plotfn: If 'plotfn = true', then a plot of the noisy data and
          estimated signal are produced.

      TI: If TI = T, then the non-decimated transform is used. See the
          help pages for wd and wst for more on the non-decimated
          transform.

 details: If 'details = FALSE' (the default), only the estimate of the
          underlying signal is returned. If 'details = TRUE', many
          other details are also returned.

  policy: Controls the type of thresholding done. Available policies
          are multiwavelet style (policy = "mws") and empirical Bayes
          (policy = "ebayes").

    code: Tells cthresh whether external C or NAG code is available to
          help with the calculations.

     tol: A tolerance parameter used in searching for prior parameters
          if the empirical Bayes policy is used.

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

     If a real-valued signal is decomposed using a complex-valued
     wavelet (like the Lina-Mayrand wavelets supplied by
     filter.select), then the wavelet coefficients are also
     complex-valued. Wavelet shrinkage can still be used to estimate
     the signal, by asking the question "which coefficients are small
     (and represent noise) and which are large (and represent signal)?"
     Two methods of determining which coefficients are small and which
     are large are proposed by Barber & Nason (2004). One is
     "multiwavelet style" thresholding (similar to that in Downie &
     Silverman (1998) where the coefficients are treated like the
     coefficients of a multiwavelet. Here, the "size" of the wavelet
     coefficient is determined as modulus of a standardised version of
     the coefficient. The standardisation is by the square root of the
     covariance matrix of the coefficient. A Bayesian method is to
     place a mixture prior on each coefficient. The prior has two
     components: a bivariate normal and a point mass at (0,0). The
     parameters are determined by an empirical Bayes argument and then
     the prior is updated by the data.

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

     Either a vector containing the estimated signal (if details =
     FALSE), or a list with the following components:  

    data: The original data as supplied to cthresh.

 data.wd: The wavelet decomposition of the data.

  thr.wd: The thresholded version of data.wd.

estimate: The estimate of the underlying signal.

   Sigma: The covariance matrices induced by the wavelet transform. See
          'make.dwwt' for more details.

   sigsq: The estimate of the variance of the noise which corrupted the
          data.

    rule: Which thresholding rule was used

  EBpars: The empirical Bayes parameters found by the function
          find.parameters. Only present if the "ebayes" policy was
          used.

 wavelet: A list with components filter.number and family which, when
          supplied to 'link{filter.select}', determine the wavelet used
          to decompose the data.

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

     Part of the CThresh addon to WaveThresh. Copyright Stuart Barber
     and Guy Nason 2004.

_N_o_t_e:

     The estimates returned by cthresh have an imaginary component. In
     practice, this component is usually negligible.

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

     'filter.select', 'find.parameters', 'link{make.dwwt}',
     'test.data', and the undocumented functions in CThresh.

