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.