GetRSSWST             package:wavethresh             R Documentation

_C_o_m_p_u_t_e_s _e_s_t_i_m_a_t_e _o_f _e_r_r_o_r _f_o_r _f_u_n_c_t_i_o_n _e_s_t_i_m_a_t_e.

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

     Computes estimate of error for function estimate. Given noisy data
     and a threshold value this function uses Nason's 1996 two-fold
     cross-validation algorithm, but using packet ordered non-decimated
     wavelet transforms to compute two estimates of an underlying
     ``true'' function and uses them to compute an estimate of the
     error in estimating the truth.

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

     GetRSSWST(ndata, threshold, levels, family = "DaubLeAsymm", 
             filter.number = 10, type = "soft", norm = l2norm, verbose = 0, 
             InverseType = "average")

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

   ndata: the noisy data. This is a vector containing the signal plus
          noise. The length of this vector should be a power of two.

threshold: the value of the threshold that you wish to compute the
          error of the estimate at

  levels: the levels over which you wish the threshold value to be
          computed (the threshold that is used in computing the
          estimate and error in the estimate). See the explanation for
          this argument in the 'threshold.wst' function. 

  family: specifies the family of wavelets that you want to use. The
          options are "DaubExPhase" and "DaubLeAsymm".

filter.number: This selects the smoothness of wavelet that you want to
          use in the decomposition. By default this is 10, the
          Daubechies least-asymmetric orthonormal compactly supported
          wavelet with 10 vanishing moments.

    type: whether to use hard or soft thresholding. See the explanation
          for this argument in the 'threshold.wst' function.

    norm: which measure of distance to judge the dissimilarity between
          the estimates. The functions 'l2norm' and 'linfnorm' are
          suitable examples.

 verbose: If 'TRUE' then informative messages are printed during the
          progression of the function, otherwise they are not.

InverseType: The possible options are "average" or "minent". The former
          uses basis averaging to form estimates of the unknown
          function. The "minent" function selects a basis using the
          Coifman and Wickerhauser, 1992 algorithm to select a basis to
          invert.

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

     This function implements the component of the cross-validation
     method detailed by Nason, 1996 for computing an estimate of the
     error between an estimate and the ``truth''. The difference here
     is that it uses the packet ordered non-decimated wavelet transform
     rather than the standard Mallat 'wd' discrete wavelet transform.
     As such it is an example of the translation-invariant denoising of
     Coifman and Donoho, 1995 but uses cross-validation to choose the
     threshold rather than SUREshrink. 

     Note that the procedure outlined above can use 'AvBasis' basis
     averaging or basis selection and inversion using the Coifman and
     Wickerhauser, 1992 best-basis algorithm

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

     A real number which is estimate of the error between estimate and
     truth at the given threshold.

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

     Version 3.6 Copyright Guy Nason 1995

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

     'linfnorm', 'linfnorm', 'wstCV', 'wstCVl'.

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

     #
     # This function performs the error estimation step for the
     # this-is-escaped-codenormal-bracket52bracket-normal function and so is not intended for
     # user use. 
     #

