heav {CVThresh}

Heavisine function

Description

This function generates Heavisine function values for n equally spaced points in [0,1].

Usage

heav(norx=1024)

Arguments

norx

the number of data or x values in [0, 1]

Details

Heavisine function is introduced by Donoho and Johnstone (1994) and is useful test function evaluating a wavelet shrinkage method.

Value

Heavisine function values f(\frac{i}{n}), i=1,\ldots,n and its variability ||f|| = \frac{\sum_{i=1}^n (f_i - \bar f)^2}{n-1} where \bar f = \frac{\sum_{i=1}^n f_i}{n}.

References

Donoho, D.L. and Johnstone, I.M. (1994) Ideal spatial adaptation by wavelet shrinkage. Biometrika, 81, 425–455.

See Also

dopp, ppoly, fg1.

Examples

testheav <- heav(1024)
plot(testheav$x, testheav$meanf, xlab="", ylab="", 
     main="Plot of Heavisine function", type="l")