waved (version 1.3)

maxithresh: Returns maxiset threshold

Description

Returns maxiset threshold for specified resolution levels.

Usage

maxithresh(data, g, L = 2, F = (log2(length(data)) - 1), eta = sqrt(6))

Value

Returns maxiset threshold for coarse resolution level equal to L, and wavelet resolution levels L,...,F.

Arguments

data

data to be processed

g

sample of g

L

Lowest resolution level

F

Finest resolution level

eta

constant describing tail of noise distribution

Author

Marc Raimondo

References

Johnstone, I., Kerkyacharian, G., Picard, D. and Raimondo, M. (2004), `Wavelet deconvolution in a periodic setting', Journal of the Royal Statistical Society, Series B 66(3),~547--573. with discussion pp.627-652.

Raimondo, M. and Stewart, M. (2006), `The WaveD Transform in R', preprint, School and Mathematics and Statistics, University of Sydney.

See Also

WaveD

Examples

library(waved)
data=waved.example(TRUE,FALSE)
maxithresh(data$lidar.noisy,data$g,L=3,F=7)