deamer-package: Non-parametric deconvolution density estimation of variables...

In deamer: Deconvolution density estimation with adaptive methods for a variable prone to measurement error

Description

The deamer package provides routines for non-parametric estimation of the density f of a variable with additive noise e. The general error model is y=x+e, where y is the noisy observation, x is the unobserved variable and e is a measurement error.

Technically, the estimation is performed using penalized deconvolution contrasts and data-driven adaptive model-selection. The estimator is projected on an orthonormal sinus cardinal basis using Fast Fourrier Transform for efficiency.
This technical framework is implemented for three situations depending on the available information or data regarding the noise e. It is assumed in all cases that the noise is homoscedastic and that its charateristic function never vanishes.
Each situation is encapsulated in a specific function:

1. deamerKE for estimation with known error density. The density of the error is assumed Gaussian or Laplace with a known mean and standard deviation (Gaussian) or scale parameter (Laplace).

2. deamerSE for estimation with an auxiliary sample of errors. This situation arises when the density of e is unknown, but an auxiliary sample of independent and identically distributed pure errors is available. Examples of such situations are found in engineering for example, when errors can be freely generated by a controlled system (like a measuring device).

3. deamerRO for estimation with an auxiliary sample of replicate observations (see deamerRO for a formal definition of replicate observations). Here, the density of the noise is specifically assumed symmetric around zero. This situation is likely to arise in biological and medical research, when pure errors cannot be observed but replicate (or repeated) noisy observations can be achieved in a sample of individuals.

Each of these functions will produce an object of class deamer for which generic methods exist. Alternatively, estimation can also be conducted using the default function deamer for users familiar with all three situations by specifying an argument 'method' and the appropriate arguments (see deamer-class and the example within).
It is worth mentionning that unlike any other deconvolution procedure, deamer does not require an estimation of a "bandwidth" parameter prior to density estimation, thus making the usage much easier. Furthermore, deamerKE and deamerSE directly handle non-centered noise.
However, none of the deamer functions are implemented for heteroscedastic errors.

Details

Package:	deamer
Type:	Package
Version:	1.0
Date:	2012-06-29
License:	GPL

Author(s)

Julien Stirnemann, Adeline Samson, Fabienne Comte with contribution from Claire Lacour.

Maintainer: Julien Stirnemann <[e-mail]>

References

Stirnemann JJ, Comte F, Samson A. Density estimation of a biomedical variable subject to measurement error using an auxiliary set of replicate observations. Statistics in medicine. 2012 May 17 [Epub ahead of print]

Comte F, Lacour C. Data-driven density estimation in the presence of additive noise with unknown distribution. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2011 Sep 1;73(4):601-27.

Comte F, Rozenholc Y, Taupin M-L. Penalized Contrast Estimator for Adaptive Density Deconvolution. The Canadian Journal of Statistics / La Revue Canadienne de Statistique. 2006; 34(3):431-52.

Comte F, Samson A, Stirnemann J. Deconvolution estimation of onset of pregnancy with replicate observations [Internet]. 2011 [cited 2011 Oct 25]. Available from: http://hal.archives-ouvertes.fr/hal-00588235_v2/

See Also

deamer-class, deamerKE, deamerSE, deamerRO

deamer documentation built on May 2, 2019, 12:36 p.m.