WaveThresh
Help
putpacket.wst
Put a packet of coefficients into a packet ordered non-decimated
wavelet object (wst).
DESCRIPTION
This function inserts a packet of coefficients
into a packet-ordered non-decimated wavelet object
(wst) object. The
wst objects are computed by the
wst function amongst others.
USAGE
putpacket.wst(wst, level, index, packet)
REQUIRED ARGUMENTS
- wst
- Packet-ordered non-decimated wavelet object into which you wish to
insert the packet.
- level
- The resolution level of the coefficients that you wish to insert.
- index
- The index number within the resolution level of the packet of
coefficients that you wish to insert.
- packet
- A vector of coefficients that you wish to insert into the
wst object. The length that the packet has to
may be
determined by extracting the same packet of coefficients using the
getpacket.wst function and using the S-Plus
length function to determine the length!
OPTIONAL ARGUMENTS
None.
VALUE
An object of class wst object containing the
packet ordered non-decimated wavelet coefficients that have been
modified: i.e. with packet inserted.
DETAILS
This function actually calls the putpacket.wp
to do the insertion.
In the future this function will be extended to insert father
wavelet coefficients as well.
RELEASE
Version 3.9 Copyright Guy Nason 1998
SEE ALSO
getpacket.wst,
putpacket,
putpacket.wp,
wst,
wst object.
EXAMPLES
#
# Take the packet-ordered non-decimated transform of some random data
#
MyWST <- wst(rnorm(1:512))
#
# The above data set was 2^9 in length. Therefore there are
# coefficients at resolution levels 0, 1, 2, ..., and 8.
#
# The high resolution coefficients are at level 8.
# There should be 256 coefficients at level 8 in index location 0 and 1.
#
length(getpacket(MyWST, level=8, index=0))
# [1] 256
length(getpacket(MyWST, level=8, index=1))
# [1] 256
#
# There should be 4 coefficients at resolution level 2
#
getpacket(MyWST, level=2, index=0)
# [1] -0.92103095 0.70125471 0.07361174 -0.43467375
#
# O.k. Let's insert the packet containing the numbers 19,42,21,32
#
NewMyWST <- putpacket(MyWST, level=2, index=0, packet=c(19,42,31,32))
#
# Let's check that it put the numbers in correctly by reaccessing that
# packet...
#
getpacket(NewMyWST, level=2, index=0)
# [1] 19 42 31 32
#
# Yep. It inserted the packet correctly.