This is exactly the same function as the wst function from wavethresh. It is only included here so that cthresh can pick up the new filter.select function.
wst(data, filter.number=10, family="DaubLeAsymm", verbose=FALSE)
- A vector containing the data you wish to decompose. The length of this vector must be a power of 2.
- 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.
- specifies the family of wavelets that you want to use. The options are "DaubExPhase" and "DaubLeAsymm".
- Controls the printing of "informative" messages whilst the computations progress. Such messages are generally annoying so it is turned off by default.
The packet-ordered non-decimated wavelet transform is more properly known as the TI-transform described by Coifman and Donoho, 1995. A description of this implementation can be found in Nason and Silverman, 1995.
The coefficients produced by this transform are exactly the same as those produced by the
wd function with the
type="station" option except in that function the coefficients are time-ordered. In the
wst function the coefficients are produced by a wavelet packet like algorithm with a cyclic rotation step instead of processing with the father wavelet mirror filter at each level.
The coefficients produced by this function are useful in curve estimation problems in conjunction with the thresholding function
threshold.wst and either of the inversion functions
InvBasis.wst The coefficients produced by the
time-ordered non-decimated wavelet transform are more useful for time series applications: e.g. the evolutionary wavelet spectrum computation performed by
ewspec. Note that a time-ordered non-decimated wavelet transform object may be converted into a packet-ordered non-decimated wavelet transform object (and vice versa) by using the
An object of class:
wst. The help for the
wst describes the intricate structure of this class of object.
Version 3.5.3 Copyright Guy Nason 1995 TEXT
# # Let's look at the packet-ordered non-decimated wavelet transform # of the data we used to do the time-ordered non-decimated wavelet # transform exhibited in the help page for wd. # test.data <- example.1()$y # # Plot it to see what it looks like (piecewise polynomial) # ## Not run:ts.plot(test.data)## End(Not run) # # Now let's do the packet-ordered non-decimated wavelet transform. # TDwst <- wst(test.data) # # And let's plot it.... # ## Not run:plot(TDwst)## End(Not run) # # The coefficients in this plot at each resolution level are the same # as the ones in the non-decimated transform plot in the wd # help page except they are in a different order. For more information # about how the ordering works in each case see # Nason, Sapatinas and Sawczenko, 1998. # # Next examples # ------------ # The chirp signal is also another good examples to use. # # # Generate some test data # test.chirp <- simchirp()$y ## Not run:ts.plot(test.chirp, main="Simulated chirp signal")## End(Not run) # # Now let's do the packet-ordered non-decimated wavelet transform. # For a change let's use Daubechies extremal phase wavelet with 6 # vanishing moments (a totally arbitrary choice, please don't read # anything into it). # chirpwst <- wst(test.chirp, filter.number=6, family="DaubExPhase") ## Not run:plot(chirpwst, main="POND WT of Chirp signal")## End(Not run)
Documentation reproduced from package cthresh, version 4.5. License: GPL-2