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:
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).
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).
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.
| Package: |
| deamer |
| Type: |
| Package |
| Version: |
| 1.0 |
| Date: |
| 2012-06-29 |
| License: |
| GPL |
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/
deamer-class,
deamerKE,
deamerSE,
deamerRO