accessD.wd: Get detail (mother wavelet) coefficients from wavelet object...

In wavethresh: Wavelets Statistics and Transforms

Get detail (mother wavelet) coefficients from wavelet object (wd).

Description

This function extracts and returns a vector of mother wavelet coefficients, corresponding to a particular resolution level, from a wd wavelet decomposition object.

The pyramid of coefficients in a wavelet decomposition (returned from the wd function, say) are packed into a single vector in WaveThresh.

Usage

## S3 method for class 'wd'
accessD(wd, level, boundary=FALSE, aspect="Identity", ...)

Arguments

wd	Wavelet decomposition object from which you wish to extract the mother wavelet coefficients.
level	The resolution level at which you wish to extract coefficients.
boundary	some methods of wavelet transform computation handle the boundaries by keeping some extra bookkeeping coefficients at either end of a resolution level. If this argument is TRUE then these bookkeeping coefficients are returned when the mother wavelets are returned. Otherwise, if FALSE, these coefficients are not returned.
aspect	The aspect argument permits the user to supply a function to modify the returned coefficients. The function is applied to the vector of coefficients before it is returned. This can be useful, say, with the complex DWT where you could supply aspect="Mod" if you wanted to return the modulus of the coefficients at a given resolution level. The default argument, "Identity", ensures that the coefficients are not modified before returning.
...	any other arguments

Details

The need for this function is a consequence of the pyramidal structure of Mallat's algorithm and the memory efficiency gain achieved by storing the pyramid as a linear vector. AccessD obtains information about where the smoothed data appears from the fl.dbase component of an wd object, in particular the array

fl.dbase$first.last.d

which gives a complete specification of index numbers and offsets for

wd.object$D.

Note that this function is a method for the generic function accessD.

Note also that this function only retrieves information from wd class objects. To insert coefficients into wd objects you have to use the putD function (or more precisely, the putD.wd method).

Value

A vector containing the mother wavelet coefficients at the required resolution level (the coefficients might have been modified depending on the value of the aspect argument).

RELEASE

Version 3.5.3 Copyright Guy Nason 1994

Author(s)

G P Nason

References

Mallat, S. G. (1989) A theory for multiresolution signal decomposition: the wavelet representation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 11, 674–693.

Nason, G. P. and Silverman, B. W. (1994). The discrete wavelet transform in S. Journal of Computational and Graphical Statistics, 3, 163–191

See Also

wr, wd, accessD, filter.select, threshold

Examples

#
# Get the 4th resolution level of wavelet coefficients. 
#
dat <- rnorm(128)
accessD(wd(dat), level=4)

wavethresh documentation built on Sept. 11, 2024, 9:33 p.m.