Heavisine function

Description

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

Usage

1
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(i/n), i=1,…,n and its variability ||f|| = ∑_{i=1}^n (f_i - mena(f))^2 / (n-1) where mena(f) = ∑_{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

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