R: Compute mean residual sum of squares for odd prediction of...

rsswav {wavethresh} | R Documentation

## Compute mean residual sum of squares for odd prediction of even ordinates and vice versa

### Description

Compute mean of residual sum of squares (RSS) for odd prediction of even ordinates and vice versa using wavelet shrinkage with a specified threshold. This is a subsidiary routine of the `WaveletCV` cross validation function. A version implemented in C exists called `Crsswav`.

### Usage

    rsswav(noisy, value = 1, filter.number = 10, family = "DaubLeAsymm",
        thresh.type = "hard", ll = 3)

### Arguments

`noisy` | A vector of dyadic (power of two) length that contains the noisy data that you wish to compute the averaged RSS for.
`value` | The specified threshold.
`filter.number` | This selects the smoothness of wavelet that you want to perform wavelet shrinkage by cross-validation.
`family` | specifies the family of wavelets that you want to use. The options are "DaubExPhase" and "DaubLeAsymm".
`thresh.type` | this option specifies the thresholding type which can be "hard" or "soft".
`ll` | The primary resolution that you wish to assume. No wavelet coefficients that are on coarser scales than ll will be thresholded.

### Details

Note: a faster C based implementation of this function called `Crsswav` is available. It takes the same arguments and returns the same values.

Two-fold cross validation can be computed for a wd object using the "cv" policy option in `threshold.wd`. As part of this procedure for each threshold value that the CV optimisation algorithm selects a RSS value must be computed (the CV optimisation algorithm seeks to minimize this RSS value).

The RSS value computed is this. First, the even and odd indexed values are separated. The even values are used to construct an estimate of the odd true values using wavelet shrinkage with the given threshold. The sum of squares between the estimate and the noisy odds is computed. An equivalent calculation is performed by swapping the odds and evens. The two RSS values are then averaged and the average returned. This algorithm is described more fully in Nason, (1996).

### Value

A list with the following components

`ssq` | The RSS value that was computed
`df` | The dof value computed on the thresholded wavelet transform of the data with the given threshold and thresholding options. (Although this is not really used for anything).
`value` | The value argument that was specified.
`type` | the `thresh.type` argument that was specified.
`lev` | The vector `ll:(nlevelsWT(noisy)-1)` (i.e. the levels that were thresholded).

### Author(s)

G P Nason

### See Also

`Crsswav`,`threshold.wd`, `WaveletCV`

[Package wavethresh version 4.6.8 Index]