propest.wav {binhf} R Documentation

## Proportion estimation function

### Description

Proportion estimation procedure for simulations.

### Usage

propest.wav(proportion = P2, binsize=1,length = 256, times = 100, meth = "u", van = 6,
fam = "DaubLeAsymm", min.level = 3)


### Arguments

 proportion A Bernoulli proportion/binomial mean function. Examples are P2, P4 and sinlog. binsize The binomial size corresponding to the mean function proportion. length Length of vector to be produced. Must be a power of two. times The number of times to sample the proportion. meth A wavelet thresholding method. Possible values are "u" for universal thresholding, or "c" for cross-validation. van the number of vanishing moments of the wavelet used in the wavelet denoiser. fam the wavelet family used in the wavelet denoiser. Possible values are "DaubLeAsymm" and "DaubExPhase". min.level the primary resolution level for the wavelet transform denoiser.

### Details

This function creates a regularly-spaced vector on the unit interval of length length, and uses these values to create corresponding values using the proportion function. These values are then used as binomial probabilities to sample "observed" binomial random variables. The observation vector is then denoised using a wavelet transform defined by the arguments van, fam, min.level with thresholding method meth. This denoising is done for both Anscombe and the Haar-Fisz method for binomial random variables. The procedure is repeated times times, and the resulting proportion estimates averaged.

### Value

 x regular grid on which the proportion function is evaluated. y vector corresponding to x of proportion function values. b matrix of dimensions timesxlength of sampled binomial variables. e matrix of dimensions timesxlength of estimated values of the proportion function, for the binomial Haar-Fisz transform. ea matrix of dimensions timesxlength of estimated values of the proportion function, for Anscombe's transform. meanfhat averaged proportion estimate for the binomial Haar-Fisz transform. meanfhata averaged proportion estimate for Anscombe's transform. amse average mean square error for the binomial Haar-Fisz transform. amsea average mean square error for Anscombe's transform.

### Author(s)

Matt Nunes (m.nunes@ucl.ac.uk)

### Examples

## Not run:
sim<-propest.wav(proportion = P2, binsize=1,length = 256, times = 1000, meth = "u",
van = 6, fam = "DaubLeAsymm", min.level = 4)

plot(sim$x,sim$y,type="l",xlab="",ylab="Binomial mean function")

##^^ shows original proportion to estimate.

lines(sim$x,sim$meanfhat,col=2)
lines(sim$x,sim$meanfhata,col=3)

##^^shows the estimates of the proportion from the two transforms.

## End(Not run)


[Package binhf version 1.0-3 Index]