scdensity: A program for self-consistent density estimation
AbstractEstimating the density of a distribution from a finite number of data points is an important tool in the statistician's and data analyst's toolbox. In their recent paper, Bernacchia and Pigolotti (JRSS-B, 2011) introduce a new non-parametric method for the density estimation of univariate distributions. Whereas conventional methods, like plotting histograms or kernel density estimates, rely on the need to make arbitrary choices beforehand (e.g., choosing a smoothing parameter), Bernacchia and Pigolotti's approach does not rely on any a priori assumptions but estimates the density in a `self-consistent' way by iteratively finding an optimal shape of the kernel. The method of self-consistent density estimation is implemented in Stata as an ado-file (-scdensity-), with its main engine written in Mata. The underlying theory and main features of this program will be discussed. In addition, results of Monte Carlo simulations will be presented comparing the performance of the self-consistent density estimate with various kernel estimates and Maximum Likelihood fits. Finally, the potential usefulness of the self-consistent estimator in other contexts such as non-parametric regression modeling will be evaluated.
Download InfoIf you experience problems downloading a file, check if you have the proper application to view it first. In case of further problems read the IDEAS help page. Note that these files are not on the IDEAS site. Please be patient as the files may be large.
Bibliographic InfoPaper provided by Stata Users Group in its series SAN12 Stata Conference with number 14.
Date of creation: 01 Aug 2012
Date of revision:
This paper has been announced in the following NEP Reports:
- NEP-ALL-2012-08-23 (All new papers)
You can help add them by filling out this form.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Christopher F Baum).
If references are entirely missing, you can add them using this form.