Brushlet Packet Transform

DESCRIPTION:

Computes the brushlet packet transform. This function is defunct.

USAGE:

bpt(x, crystal.names=NULL, basis=NULL, partition=NULL, cost.fun="entropy",  
    n.levels=6, taper="poly2", boundary="periodic", 
    n.taper=NULL, scale=NULL, thresh=NULL, p=2, prob=.5) 
ibpt(x) 

REQUIRED ARGUMENTS:

x
a time series or a numeric vector, should demean first, or, for ibpt only, an object of class bpt.

OPTIONAL ARGUMENTS:

crystal.names
a character vector of crystal names of the form "bj.i" where j is an integer giving the level and i is an integer from 0 to 2^j-1 giving the block number within level j. If crystal.names is supplied, then basis is ignored.
basis
a character string specifying how the cosine packet basis is selected: "best.basis", "best.level", "block.bpt", or "gbpt". See the corresponding function for details.
partition
a vector specifying the desired partitioning of the data if the "gbpt" basis is used. See the help for function "gbpt" for details on how to specify the partition.
cost.fun
character string indicating which cost functional to use: "entropy", "threshold", "sure", and " lp" are available. Only used for "best.basis". See pcosts for details.
n.levels
a non-negative integer specifying the blocking factor: x is divided into 2^n.levels blocks each of length length(x)/2^n.levels. For "best.basis", n.levels gives the blocking factor for the finest level. If n.levels is bigger than ml, where ml is the maximum possible level, computed from the max.level function, then n.levels is set to ml and a warning message is given.
taper
a character string, indicating the taper function: "boxcar", "poly1", "poly2", "poly3", "poly4", "poly5", or "trig". See the function bp.table for details.
boundary
a character string, available boundary rules: "reflect", "periodic" and "zero". See the function bp.table for details.
n.taper
a non-negative integer. The length of the taper will be 2*n.taper. By default, n.taper is set to length(x)/2^(n.level+1), which is the maximum possible length at the finest blocking level.
scale
scaling factor for computing the costs. The default depends on cost.fun. See the function pcosts for details.
p
a number in the interval (0,2] giving the degree of the lpnorm when cost.fun is "lp". See the function pcosts for details.
thresh
a non-negative number giving the threshold for when cost.fun is "threshold" or "sure". See the function pcosts for details.
prob
a number in the interval (0,1) used to compute the threshold for when cost.fun is "threshold". See the function pcosts for details.

VALUE:

an object of class bpt, inheriting from the class crystal.vector.

ibpt(x) performs the inverse transform from the bpt object x, as the generic function reconstruct does for bpt objects.

DETAILS:

The algorithms for the taper functions are given in the S+WAVELETS User's Manual, in the section "Cosine Packet Algorithms". They are discussed in greater depth in Wickerhauser (1994). The default optional arguments n.levels, taper, boundary can be reset using function wavelet.options, see wavelet.options for details.

REFERENCES:

Wickerhauser, M. V. (1994). Adapted Wavelet Analysis from Theory to Software. A. K. Peters Ltd, Wellesley, MA.

Meyer, F.G. and Coifman, R.R. (1997), Brushlets: A tool for directional image analysis and image compression Applied and Computational Harmonic Analysis, Academic Press Publishers.

SEE ALSO:

, , , , , , , .

EXAMPLES:

x <- make.signal("linchirp") 
bbx <- bpt(x, basis="best.basis") 
eda.plot(bbx)